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ABSTRACT 


We  derive  the  closed  form  expression  for  the  probability  of  bit  error  in  direct  detection, 
dense  wavelength  division  multiplexed  (WDM)  fiber  optic  systems  employing  OOK 
as  the  modulation  technique  and  single-cavity  Fabry- Perot  (FP)  filters  in  the  receiver 
as  demultiplexers.  The  expression  is  derived  in  the  time  domain  using  the  impulse 
response  of  the  single-cavity  Fabry-Perot  filter  and  the  complex  baseband  equivalent 
of  the  received  dense  WDM  optical  signal.  The  two  are  convolved  to  produce  the  FP 
filtered  output  signal  s{t).  We  then  integrate  over  one  bit  period,  where  1Z 

is  the  responsivity  (A/W)  of  the  photodetector  following  the  FP  filter  in  the  receiver 
structure.  This  integral  is  the  deterministic  X  of  the  decision  variable  Y  where 
Y  =  X  +  N.  N  is  the  postdetection  thermal  noise  (amplifier  generated),  a  zero  mean 
Gaussian  random  variable  with  variance  NqT  where  No  is  the  noise  current  spectral 
density  (A^/Hz).  Both  X  and  N  are  combined  into  an  expression  for  probability  of 
bit  error.  A  limited  case  of  the  complete  model  is  assumed,  and  probability  of  bit 
error  graphs  are  generated. 
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I.  INTRODUCTION 


Direct  detection  optical  frequency  division  multiple  access  (FDMA)  networks 
are  increasingly  becoming  an  attractive  alternative  to  coherent  optical  FDMA  net¬ 
works  [1].  One  of  the  primary  reasons  is  that  noncoherent  systems  do  not  require 
expensive  synchronization  circuitry  for  proper  operation.  Also,  present  optical  filter 
technology  allows  designers  to  closely  pack  the  channels  in  frequency,  resulting  in 
dense  wavelength  division  multiplexed  (WDM)  systems  that  can  provide  aggregate 
bit  rates  of  many  terabits  per  second  (IT  b/sec  =  10^^  bits/sec)  [2].  The  possible 
uses  of  dense  WDM  systems  are  many  (local  area  networks,  undersea  surveillance, 
etc.),  but  one  can  easily  see  the  enormous  economic  benefit  and  importance  of  being 
able  to  transmit  aggregate  bit  rates  of  terabits  per  second  on  a  single  fiber  without 
the  economic  burden  of  expensive  synchronization  circuitry  imposed  by  a  coherent 
system.  In  this  thesis,  we  derive  the  complete  closed  form  expression  in  the  time 
domain  for  the  probability  of  bit  error  of  dense  WDM  systems  employing  OOK  as 
the  modulation  technique  and  single-cavity  Fabry-Perot  (FP)  filters  as  channel  de¬ 
multiplexers.  In  our  derivation,  we  make  no  simplifying  mathematical  assumptions 
and  use  the  impulse  response  of  the  single-cavity  Fabry-Perot  filter,  which  is  an  in¬ 
finite  sum  of  delayed  impulses  whose  intensities  decrease  geometrically.  The  results 
of  our  work  are  presented  in  Chapters  II,  III,  and  IV.  A  detailed  derivation  of  the 
deterministic  signal  component  of  the  decision  variable  appearing  at  the  output  of 
the  integrator  of  the  Channel  0  (Channel  of  Interest)  receiver  appears  in  Appendix  A. 
See  Fig.  1  for  a  sketch  of  this  receiver.  Appendix  B  shows  how  the  complete  model 
is  reduced  for  a  limited  case  and  how  probability  of  bit  error  calculations  are  made 
for  this  limited  case.  Appendix  C  shows  the  strategy  and  programs  used  to  obtain 
the  probability  of  bit  error  for  four  values  of  free  spectral  range-bit  period  product 
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and  a  given  range  of  signal-to-noise  ratios.  It  is  interesting  to  note  at  this  point 
that,  although  the  model  presented  in  this  thesis  is  strictly  derived,  complete,  and 
done  without  approximation,  its  major  weakness  is  that  it  is  computationally  very 
intensive.  In  fact,  it  takes  several  months  of  computer  time  on  several  SPARC-10 
workstations  working  simultaneously  to  generate  a  single  graph  of  the  probability  of 
bit  error  for  this  system.  At  the  end  of  the  thesis  we  conclude  that,  although  our 
model  is  mathematically  correct,  a  discrete  time  approach  to  the  problem  is  probably 
a  more  efficient  investigative  tool  in  analyzing  the  performance  of  dense  WDM  fiber 
optic  networks. 


Photodetector  Output  Current  =  7J|s(t)p 


X  =  'Jlj^\s{t)\'^dt  N  =  j\{t)dt 

Figure  1:  Channel  0  OOK  receiver  structure. 


DATA 

- ► 
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II.  ANALYSIS 


An  OOK,  dense  WDM  system  utilizing  single-cavity  Fabry-Perot  (FP)  optical 
filters  as  demultiplexers  consists  of  Af  -|- 1  transmitters  (fixed  wavelength  lasers)  con¬ 
nected  over  a  fiber  link  to  M  +  l  receivers.  Each  receiver  contains  a  frequency  selective 
Fabry-Perot  filter  to  demultiplex  one  of  the  M  +  l  channels.  In  our  derivation  of  the 
closed  form  expression  for  the  probability  of  bit  error  for  this  system,  M  is  the  even 
number  of  adjacent  channels  symmetrically  placed  in  frequency  around  the  Channel 
of  Interest  (Channel  0)  whose  carrier  is  transmitted  on  an  arbitrary  wavelength  \a- 
The  Channel  0  receiver  is  shown  in  Fig.  1.  The  dense  WDM  optical  signal  is  received 
by  the  Fabry-Perot  filter  tuned  to  Aa,  which  allows  the  Channel  0  data  signal  to  pass 
and  rejects  signals  in  the  adjacent  channels.  A  photodetector  then  converts  the  fil¬ 
tered  optical  signal  to  a  current.  The  photodetector  has  an  arbitrary  responsivity  1Z 
(A/W).  The  current  is  then  amplified  by  a  low-noise  amplifier  which  adds  a  postde¬ 
tection  thermal  noise  n{t)  with  two-sided  current  spectral  density  Nq  (A^/Hz).  Both 
the  signal  and  noise  currents  are  now  passed  to  an  integrator  (low  pass  filter).  The 
output  of  the  integrator  is  the  decision  variable  F,  which  is  compared  to  a  threshold 
Vt  to  determine  whether  an  “0”  or  “1”  was  sent. 

A.  THE  DENSE  WDM  SIGNAL  AT  THE  INPUT  OF  THE  FABRY- 
PEROT  FILTER:  r{t) 

Again  we  note  that  Channel  0  is  our  channel  of  interest.  We  can  write  the 
expression  for  the  Channel  0  data  signal  in  the  complex  baseband  as 

bo{t)=  E  bo,iPT{t-iT)  (1) 

where  T  is  the  data  bit  period  (s),  and  6o,,-  is  the  bit  in  Channel  0  during  the  time 
period  [iT,  (i  +  l)r].  Note  that  6o,t  G  {0, 1).  To  is  a  positive  integer  which  represents 
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the  number  of  bits  in  Channel  0  that  are  trailing  the  detected  bit  or  bit  of  interest 
^0,0-  Pt(0  rectangular  pulse  function  defined  as 


Pt(0  = 


1, 

0, 


0  <t  <T 

otherwise 


(2) 


Channel  k  is  any  arbitrary  adjacent  channel.  There  are  M  adjacent  channels  placed 
symmetrically  in  frequency  around  Channel  0.  M  is  an  even  integer.  The  indices  for 
k  are  as  follows  k  =  —M/2,  •  •  • ,  — 1, 1,  ■  ■  • ,  +M/2.  We  may  now  write  the  complex 
baseband  expression  for  the  k*'^  channel  data  signal 

bk{t)=  E  (3) 

l=-L 

where  Uk  is  the  radian  frequency  spacing  between  Channel  0  and  Channel  k  and 
uJk  =  We  have  already  defined  the  pulse  function  prit)  in  Eq.  (2)  above.  L  is 

a  positive  integer  which  represents  the  number  of  bits  in  Channel  k  which  trail  the 
bit  in  Channel  k  that  is  the  bit,  b^fi.  We  now  note  that  b^/  is  the  bit  in  Channel 
k  during  the  time  period  [£T,  {£  +  IjT]  and  that  6*,^  G  {0,  where  4>k  is  Channel 

k's  phase  offset  from  Channel  0.  4>k  is  assumed  to  be  a  uniformly  distributed  random 
variable  between  0  and  2x  ~  t/[0,27r]).  We  can  now  use  Eqs.  (1)  and  (3)  to  write 
the  complex  baseband  equivalent  dense  WDM  received  optical  signal  which  appears 
at  the  input  of  the  Channel  0  Fabry- Perot  filter  (see  Fig.  1) 

+A//2 

r{t)  =  Vp  bo{t)  -f  ^  y/P  bk{t)  (4) 

k=-M/2 

k^O 


and 

/  0  +M/2  0 

>■(0  =  ^(  E  V.PT(i  -iT)+  E  E  -  IT) 

k=-M/2 

k:^0 

where  P  is  the  received  optical  power. 


(5) 
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B.  THE  OPTICAL  OUTPUT  OF  THE  FABRY-PEROT  FILTER: 

s{t)  =  ro(0 

To  arrive  at  the  optical  output  of  the  Fabry-Perot  filter  s{t),  we  must  convolve 
the  input  signal  r{t)  with  the  impulse  response  h[t)  of  the  filter 

s{t)  =  ro{t)  =  r{t)  *  h{t)  =  h{t)  *  r{t)  (6) 


First,  however,  we  need  the  expression  for  h{t).  The  Fabry-Perot  filter  is  a  causal, 
linear,  time-invariant  (LTI)  system.  It  can  be  shown  [1]  that 


h{t)  =  {I  -  p)f^  p^6  ft- I 


(7) 


5=0 

where  8{t)  is  the  Dirac  Delta  Function  defined  in  two  parts 

8{t)  = 

and 


0, 

0, 


t  <  0 
t  >  0 


(8a) 


8{t)dt  =  1  (8b) 

Jo- 

Looking  at  Eq.  (7),  we  note  that  p  is  the  power  reflectivity  of  the  single-cavity  Fabry- 
Perot  filter  and  P  is  the  filter’s  free  spectral  range  (Hz).  Performing  the  convolution 
operation  yields  the  following  (see  Appendix  A  for  details) 


s(t)  =  roit)  =  \/P  {I  -  p) 


K 


^  p^bo^oPT  (  t  - 


Lj9=0 


—  1  oo  / 

+  P^KiPx  ( t 

t=-Lo  5=0  ^ 


/5 


«b(^) 


iT 


+M/2 


sisi(t) 
0  oo 


+  E  EE  it-^-iT' 

k=-M/2  ^  ^  ^ 

k^o 

' - V - 

MC/W 


(9) 
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where  5b (t)  is  the  desired  signal  (signal  of  interest),  sisi{t)  is  the  intersymbol  inter¬ 
ference  signal,  and  SAci{t)  is  the  adjacent  channel  interference  signal.  Then 

s{t)  —  ro{t)  =  K  (10) 

We  are  now  interested  in  writing  an  expression  for  s{t)  during  the  detection 
interval  0  <  t  <  T.  We  do  this  by  sketching  SB{t)  from  Eq.  (9)  and  then  writing 
another  expression  which  accounts  for  how  the  rectangular  pulses  behave  during  the 
detection  interval.  We  then  sketch  the  pulse  functions  for  the  first  three  values  of  i 
for  sisi{t)  in  Eq.  (9).  We  note  how  they  appear  during  the  detection  interval,  and 
recognizing  a  pattern,  we  write  another  expression  for  Sjsiit)  for  0  <  t  <  T.  A  similar 
process  is  followed  to  arrive  at  SACi{t)  for  0  <  t  <  T.  We  now  present  SB{t),  5/s/(t), 
and  SAci{t)  for  0  <  t  <  T  (see  Appendix  A  for  details) 


g-1 

■5b(0  =  P^bofiPT 

3=0 


rji  d 


for  0  <t  <T  (11) 


where  q  =  Tf{lf  13)  =  I3T  is  an  integer.  Recall  that  T  is  the  data  bit  period  (s),  and 
P  is  the  single-cavity  Fabry-Perot  filter’s  free  spectral  range  (Hz). 


1  ‘'A'  .1  (  fiT(t) 

Sisi{t)  =2^  2^  P^bo,iPT  — ,  ,  -x 

i=-Lo  [g  =  -(l+,),+  l  +  (1  +  0? 


-(t-l)g-l 

+  P^bo,iPT 


9  +  iq 

13 


9  +  iq 

(3 


,  for  0  <  t  <  r  (12) 


6 


/ 


+Af/2 

SACl{t)  =  ^ 

k=-M/2 

fc/0 


/ 


Lp=0 


V 


t-^ 


T- 


+  E 


e=-L 


i-eq)-i 

S=-(l+^)q4-l 


vv  /5; 

/?r(0 


v5  + (1  +  ^)9. 


-(^-i),-i 

g--iq 


( 


Pt 


V 


t  - 


'9_±^ 


\ 


/ 


,  for  0  <  i  <  r  (13) 


/ 


Recall  Eqs.  (9)  and  (10)  where 

s{t)  =  K{sB{t)  +  S[si{t)  +  SACl{i)) 


and 

k  =  Vp{i-p) 

We  can  see  that  s{t)  in  the  interval  0  <t  <T  can  now  be  calculated  by  substituting 
Eqs.  (11),  (12),  and  (13)  into  their  appropriate  positions  in  Eq.  (10). 

C.  COMPUTATION  OFX  =  n  /^ls(t)pdt,  THE  SIGNAL  COMPONENT 

Jo 

OF  THE  DECISION  VARIABLE 

Looking  at  Fig.  1,  we  see  that  the  filtered  optical  output  of  the  Fabry-Perot 
filter  s{t)  is  passed  to  a  photodetector.  The  output  current  of  the  photodetector 
is  7^|5(t)p  where  TZ  is  the  responsivity  of  the  photodetector  (A/W).  Thus,  we  can 
see  that  the  deterministic  signal  component  of  the  decision  variable  at  the  integrator 
output  is 

X  =  n  r  \s{t)fdt  (14) 

Jo 
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We  have  already  derived  SB{t),  sisi{t),  and  SACl{t)  for  0  <  f  <  T.  Directly 
substituting  these  three  signals  into  Eq.  (10)  yields  s{t)  for  0  <  <  <  T.  Now  our  true 
task  is  to  compute  |5(t)P  and  \s{t)\'^dt  as  TZ  is  an  arbitrary  constant. 

s{t)  =  '/P  (1  —  p)[Sjg(t)  +  sisi{t)  +  (15) 

^ ^  - - - ^ 

K  s^{t) 

Then 

|s(t)p  =  K^s\t)s'{ty  (16) 

and  dropping  the  {t)  notation  for  convenience  in  the  three  terms  of  s'{t)  and  noting 
that  complex  conjugation  is  a  linear  operator 


S'{t)*  —  Sg  +  S*jgj  T 


’ACI 


So 


S  {t)s  (t)  —  [sb  +  Sisi  +  SaCi){Sb  +  ■S/5/  +  Saci) 


Multiplication  yields 


S  —  SbS*b  +  SBS*jg]r  -)-  sbs*^ci 

+SlSlS*B  +  Sisrs}si  +  SISIS*ACI 
+SacIS*b  +  SaCIS}si  +  ^ACIS*aci 


Rearranging  and  simplifying  yields 


(17) 


(18) 


(19) 
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|5(i)|' =  A^V(ty(o*  = 


sbs*b  +  srsis*isi  +  ^aciSaci 


|sbP  |s/S/P  |s/4C/P 

+  SbS*jsj  +  SisiSg 


2Re[sBs}si]  =  2Re[sisiSB] 


+  ^bSaci  +  sacis*b 


2Re[s*BSACi]  =  2Re[sBS*Aci. 


+  ^ISlS*ACl  +  ^AClS*[si 


2Re[sj^jSAci,  —  ^R^I^isi^aci]  (20) 


We  will  integrate  each  of  these  terms  over  0  to  T  to  compute  Jq  \s{t)\‘^dt. 

Substituting  the  appropriate  definitions  into  the  six  terms  of  Eq.  (20)  yields 
21  terms  (“clusters  of  summations”)  which  are  integrated  from  0  to  T.  Also  note 
that  we  had  to  create  several  “gating”  functions  to  turn  off  integrals  when  the  pulse 
function  products  in  some  of  the  21  terms  fail  to  overlap  for  certain  summation 
indices.  Although  intuitively  obvious,  see  Appendix  A  [below  Eqs.  (104)  and  (112)] 
for  definitions  of  the  max  and  min  functions  which  appear  in  some  of  the  forthcoming 
equations. 

Recall  that  K  =  y/P  {1  —  p)  and  —  P{\  —  p)^.  Then  we  have 

r  \s{t)\^dt = K'^  r  s'{t)s'{tydt = r  is'(i)|2dt  (21) 

0  *'0  •/  0 

which  has  the  following  terms  [see  Eq.  (20)] 


/(•"  r \,sfdt  =  P(1  - pfbioY:  E  (t -  { 

g=0m=0  V  V 


max(p,m)' 


(22) 
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(23) 


where 


in  which 


[  \sisi\^dt  —  Aisi  +  Bisi  +  Cisi 

Jo 

-1  (-*9)-l  -1  {-rq)-l 

A,s,  =  p{i-prY:  EE  E  v’o 

i=z—Lo  g=-{l'^i)q-\-l  r=-Lo  (14-r)g+l 


-1  (~tg)-l  -1  -(r~l)q-l 

B,s,  =2P(l-rt\E  E  E  E  /+’“VA,G.(s,i,m,r) 

i=z~-Lo  — (l  +  2)g-f  1  T~~Lci  m—  —  rq 

in  which 

[  'e[(5' +  (1  +  0?)  “  for  m  +  +  (1  + 


Gi{g,i,m,r)  = 


otherwise 


-1  -(t-l)g-l  ^1  -(r-l)g-l 

C/s;  =  P(l-^)"  E  E  E  E  P’*’"b„M.r  T 

i^-^Lo  g^—iq  r=  —  Lo  m^^—rq 


ma,x{g  +  iq,m  +  rq) 


where 


f  \sACi\^dt  —  Aaci  +  Baci  +  Caci  +  Daci  +  Eaci  +  Faci  (24) 

J  0 


+M/2  g  — 1  -{'Mf2  9— 1 

A^ci  =  p(i-pr  E  E  E 

k=-Mf2  n=-M/2 


in  which 


^k,0^n,0  ■(  ^  ^  ' 

’  3\^k  -^n) 

\pi[(bOk-'^n)T-uj};(glP)+Wn{m/g)\  _  „j[(a>fc-a;„)(l//3) max(5,m)-u;j:(p//3)+i»;n(m//3)] 


,  ioi  k  ^ 


ioT  k  =  n 
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4-A//2  -1  {-iq)-l  +M/2  -1  (-rg)-l 

B,c,  =  p{i-P?  E  E  E  E  E  E 


A:=:-M/2 

A:^0 


e—Lg=-{l+e)q+l 

n^O 


r=—L  m=  — (lH-»')g+l 


in  which 


J[-‘^fc(s//3)+‘^n(m//3)] 


^ei[(‘^fc-‘*'n)(l//3)rain(s+(H-^)g,m+(l+r)9)]  _  ^ 


for  k  ^  n 


hM,T  ^  min(i?  +  (1  +  ^)9,  m  +  (1  +  r)q)  ^  , 


for  k  =  n 


+M/2  -1  -(^-1)?-1  +M/2  -1  -(r-l)9-l 

c,c,  =  p(i-pr  E  E  E  E  E  E 


in  which 


[(t^lt  -i^n  )T-w;t  (g//3)  +c.;„  (m//3)] 


^k,eK,r  ~( - \  1^ 


_  J  _gi[(c*'it-Wn)(l//3)max(5+^9,Tn+r9)-a2fc(5//3)+w„(m//3)]|  ^  £qj.  ^  ^ 


T  —  -n  max(5f  +  iq,m  +  rq)  ,  for  A:  = 
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-|-iVf/2  ^—1  -\-Mf2  -1  — (t-  — 1)7  — 1 

D,c;  =  Pii-p?  E  E  E  E  E  ■^'4 


k.-M/2  »=”  n=-M/2  '=-‘-  "=■'» 


in  which 


_  \^3[{^k-^n)T-ujk{g/P)+ujn(m/P)] 


oj[{^k-^n)(l/P)Tnax(g,m-{-Tq)-uJk(g/P)-\-ujn{m/P)] 


-\-ri9+rnl^*  A  _  \ p-j[(^k-i^n)T-ujk{g/P)+(^n(7n/P)] 

“r  ^k,0^rL,r  r  •/  M  ^ 

-^nj]  ^ 


—  p-jl{^k-<^n)il/P)  max(p,m+7'g)-t^lc (5//?)+<^n (m/^)] 


for  A:  7^  n 


T-  -m&x{g,m  +  rq) 


for  k  =  n 
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-l-M/2  q-1  +M/2  -1  (-rg)-l 

EAci^Pii-pf  E  E  E  E  E 

k=-M/2  n=-M/2  ^=-£'m=-(l+r),+l 


in  which 


j{u!k  -UJn  )  '■ 


J[{uk-‘^n)((m+{l+r)q)/0)-uJk{3/0)+‘^n(m/P)] 


—pjl{<^ri/P}{m-g)] 


,-j[i(^k-<^n){{m-{-{l+r)q)/P)-uJki9/P)-\’i^n{'m/P)] 


_e  j[i^n/P)(m  5)]  j  (72(5, m,r), 


for  k  ^  n 


-[(m  +  (1  +  r)^)  -g]  \  ^2(5^,  m,  r), 


for  k  =  n 


where 


G2(</,m,r)  = 


1,  for  <  m  +  (1  +  r)q 

0,  otherwise 
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+M/2 

-1  (--f?)-! 

+M/2 

-(r--l)g-l 

E 

E  E 

E 

E  E 

:=-M/2 

n=-Af/2 

r~-‘L  m  =—rq 

n^O 

in  which 


=  < 


oi[{‘^k-<^n){(g-k-(i-\-(-)q)/li)-'^k(9lP)+'^n{rn/P)] 


^pi[{<^k-<^n)((m+Tq)l(3)-(jj^(gl0)+tjjn(mll3)\ 


,-j[(‘^*-‘^n)((3+(l+^)g)//?)-Wfc(fl//3)+Wn(m//3)] 


p9+m  ^  .  .  . 


where 


[-j{oJk  -  a;„)] 

_g-i[(w/c-a/„)((m+rq)//?)-u;fc(3//3)+a-„(m//3)]J^ 

+  (1  +  ^)9)  -irn  +  ^9)]^  Gsig,  £,  m,  r), 

f  1,  m  +  r^  <5r  +  (l  +£)9  ) 


ioT  k  ^  n 


for  k  = 


n 


[  0,  otherwise 

f  2Re[sBs}gj]dt  =  (  2[sBSisi]dt  = 

Jo  Jo 


p^+™6o,o^o,iG4(5,z,m) 

g=Oi=-Lo  77i=“(l4-t)q+l 


q-l  _1  _(t_l)g_l 

+  2P{i-pfY,  T.  E  /-""VVi 

g—Oi=—Lo  m—  —  iq 


T  —  —  max(5f,  m  +  iq) 


(25) 


where 


Gi[g,i,m)  =  { 


^[{m +  i)q)  -  g\,  g  <  m  +  {I -\- i)q 


0, 


otherwise 


14 


[  2Re[s*gSAci]dt  =  f  2Re[sBSAci]dt  =  Ab-aci  +  Bb-aci  +  Gb-aci 
Jo  Jo 


where 


9-1  +A//2  9-1 

Ab-aci  ^  13  13  ¥^7 

5=°  k=-M/2 


in  which 


^7 


J(^k 


+P^'^"bo,oblo 


[-j^k] 


o-j^k{T-(m/l3))  _  -jaifc((max(9,77i)-m)//J) 


9-1  +M/2  _1  (-<9)-l 

Bb-aci  =  ^’(1  “  13  13  13  13 

it=-M/2  ^=-'^^™=-(l+^)9+l 
ki^O 


in  which 


(^8  =  — 


oi‘^fc((l+09//5)  _  p3'^k{(3-m)IP) 


Gz{g,(-,m) 


[-j^^k] 


G5{g,i,rn) 


where 


Gs{gJ,m)  = 


’  1,  g  <m-\-{l+i)q  ] 

0,  otherwise 


9-1  +M/2  -1  -(^-1)9-1 

Cb-aci  =  ^(1  -  13  13  13  13  *^9 

fe=-M/2  ^=-'^ 


(26) 
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in  which 


<)^9  — 


+P^-^'^ho,oh\ 


kl 


[~]^k\ 


Ji^k(T-{m/P))  _  juik{(mdLx{g,Tn+eq)—m)/(3) 


,-i'^k{T-{mlP))  _  -j<A)k((ma.x(g,m-\-Cq)—m)l0) 


I  2Re[s*j^jSAci]dt  —  f  2Re[sisiSACi]dt 
Jo  Jo 

—  J^ISI-ACI  +  BiSI-ACI  +  Cjsi-ACl 
+  D1SI-ACI  +  EisI-ACI  +  Fisi-ACI 


where 


in  which 


V’lo 


-1  +M/2  q-1 

Aisi-aci  P{1  -  pf  ^  ^  Y1  XI  ‘<^10 

z=-Lo  g=-(H-t)?+l 

it  540 


+p^^'"bo,iblo7Z^]  -  1]  Ge{g,i,Tn) 
[  3^k\  ^ 


where 


G6{g,i,m) 


1,  m 


0,  otherwise 


<  5r  +  (1  + 


-1  (-t?)-!  +M/2  -1  i-eq)-l 

Bisi-aci  =  ^(1-p)'  E  E  EE  E 

i=-Lo  s=-(l+t)g+l  k=-M/2 
ki^O 


in  which 


V^ii  —  P^^'^bo^ibk,^ 


J^k 


Ji^k[{'i-/P)rr>in{g+(l+i)q,m+{l+i)q)-m/P]  _  juik[-m/P] 


’*  ’  [-j^k]  ^ 


(27) 
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in  which 


<pn  — 

where 


in  which 

V^13 


in  which 

= 


-1  +M/2  -1  -(^-l)g-l 

Cisi-Aci  =  -P(i  -  pY  'Yh  Y  Y  Y  Y 

t=->Lo  5=  — (l+z)q+l  k——M/2  m=—iq 

k:^0 


gi‘^fc[(l//3)(s+(l+«)9)-™/^]  _  gi")t[^9//5] 


+P^^^Kbh 


1 


[-JWfcJ 


G7{g,iJ,m) 


G7{g,iA,m) 


I"  1,  m  +  iq  <  g  {I -\-  i)q 
[  0,  otherwise 


-1  -(t-l)9-l  +MI2  9-1 

Disi-aci  =  p{^-pf  Y  Y  Y  Y  V^13 

k=-M/2 
k^O 


=  p^'^”^boibko-—  _  giwfc[(l/^)max(ff+!g,m)-m//?] 


J>^k 


in9+mi  L*  \p-3<^klT-{m/0)]  _  -jwk[{l/P)ma.x{g+iq,m)-m/p] 


-1  -(«-l)9-l  +M/2  -1  (-^g)-l 

Eisi-aci  =  P{^  ~  pY  Y  Y  Y  Y  S  *^14 

i=-Lo  g=-iq  k=z-Mt2  ^=“■-'^’^=“‘(1+^)9+! 
k^O 


P^^'^bo,ibk,(- 


Ji^k 


[-Ml 


=>iwfc[(l//?)(l+^)g]  _  „ii*'fc[(l//3)(s+*g)-m//3] 


Gs{g,iJ,m) 
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where 

f  1,  g  -{■  iq  <  m  +  {I  +  l)q 

Gs{g,i,i,m)  =  < 

0,  otherwise 

-1  +M/2  -1 

Fisi-aci  =  P(1  -  /j)^  ^  X]  S  Fi5 

i=-Lo  g=-iq  k=-M/2 
k^O 

in  which 

—  p^'^^hpibk  i  ^  Qj^k[T-{m/p)]  _  ^jwk[{l/0)ma.xig+iq,Tn+iq)-m/0] 

’*  ’  jl^k  ^  J 

_|-^9+™5g  .5*  _ i _  ^-:'^k[T-{m/P)]  _  ^-jwk[{l/P)max{g+iq,m+iq)-m/0]  /2g) 

’’  ’  [-M]  ^ 

We  have  now  laid  down  all  terms  of 

rT  {'T  X* 

/  \s{t)\^dt  =  K^  s'{t)s'{tydt  =  [  \s'{t)\^dt 

Jo  Jo  Jo 

This  will  be  used  to  compute  the  deterministic  signal  detection  statistic  X  where 

X  =  7^  r  \s{t)\^dt 
Jo 

This  will  be  used  to  compute  probabilities  of  bit  error  for  the  dense  WDM  system. 

D.  THE  DECISION  VARIABLE  Y 

The  decision  variable  Y  appears  at  the  output  of  the  integrator  (see  Fig.  1). 
It  consists  of  a  signal  component  X,  which  was  presented  in  Section  C,  and  a  noise 
component  N 

Y  =  X  +  N  (29) 

where 

T 

X  I  \s{t)\^dt  (30) 

Jo 
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and 


N=  r  n{t)dt  (31) 

Jo 

We  note  that  n{t)  is  a  zero  mean  Gaussian  random  process  with  two-sided  current 
spectral  density  Nq  (A^/Hz),  and  is  a  zero  mean  Gaussian  random  variable  with 
variance  NqT.  Y  is  compared  to  a  threshold  Vr  to  determine  whether  a  “0”  or  “1” 
bit  is  presented  at  the  output. 


E.  PROBABILITY  OF  BIT  ERROR  FOR  THE  DENSE  WDM 
SYSTEM 


For  a  detection  threshold  Vr  and  an  ACI/ISI  bit  pattern  ■0?  =  {bk,(,bo,i}  for 
k  —  -M(2,  •  •  • ,  +MI2,  k  ^  0,  £  =  -A,  •  •  • ,  0,  and  i  =  -Aq,  •  •  • ,  -1;  the  conditional 
probability  of  bit  error  for  the  dense  WDM  system  utilizing  the  decision  variable  Y 
described  by  Eqs.  (28)-(30)  is  given  by  [2,  4] 


1  (Vt-X^ 

2^1  v%r  , 


(32) 


where 


Q(x)  =  ^j"  (33) 

Xq  is  the  value  of  X  given  by  Eq.  (14)  or  (30)  with  the  value  of  \s{t)\'^dt  obtained 
by  summing  Eqs.  (22)-(27)  when  6o,o  =  0-  is  obtained  in  the  same  manner  with 
bo  o  =  1.  The  average  probability  of  bit  error  Pb  is  given  by  taking  the  expected  value 
of  Pe(0p)  given  in  Eq.  (32)  over  all  possible  bit  patterns  ipp  [2]-  Let  us  define  an 
ACI/ISI  bit  pattern  set  0  =  {0p}j  P  —  1?' ' '  ,NPAT.  NPAT  is  the  total  possible 
number  of  bit  patterns  in  the  set  0  =  {0p}'  Then 


a  =  (34) 

Vl’p) 

If  we  count  up  the  number  of  independent  bits  in  0p,  there  are  M(A  4- 1)  •+•  Aq  bits. 
Using  the  assumption  that  a  0  or  1  is  equally  probable,  we  see  that  there  are  two 
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possible  ways  to  fill  each  bit  position.  Then,  by  the  generalized  multiplication  rule,  if 
we  let  NPAT  be  the  total  number  of  possible  bit  patterns  in  the  set  ip  =  {ipp} 

NPAT  =  (2),  (2)2  •  • .  (35) 

If  we  assume  that  all  bit  patterns  mtp  =  {tpp]  are  equiprobable  (an  excellent  assump¬ 
tion),  then  by  the  law  of  total  probability 

~  2A7(I+T)+rT-^'=(^2)  H  1-  2M(l+i)-i-Lo  ^^(^npat) 

and 

1  NPAT 

p=i 
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III.  NUMERICAL  RESULTS 


We  attempt  to  generate  four  probability  of  bit  error  graphs.  We  limit  our 
expression  for  X  given  by  Eq.  (14)  or  (29)  and  the  sum  of  Eqs.  (22)-(27)  by  setting 

Lq  =  1  L  —  0  <f>k  —  0 

and 

2'KkI 
=  -Y- 

where  k  is  an  integer  and  /  is  the  normalized  channel  spacing  integer  (/  >0).  When 
we  do  this,  our  21  terms  (“clusters  of  summations”)  reduce  to  ten  terms  (“clusters 
of  summations”).  In  the  interests  of  brevity,  we  refer  the  reader  to  Appendix  B  for 
the  details  on  how  we  came  up  with  the  ten  programmable  terms  of  X  for  this  case. 
Referring  to  Appendix  B,  the  probability  of  bit  error  equations  we  use  to  generate 
these  graphs  mirror  the  equations  presented  in  Chapter  II,  Section  E.  Referring  to 
Appendix  C,  one  finds  the  programming  strategy  for  the  ten  terms  and  the  final 
graphs,  the  actual  programs  for  each  of  the  ten  terms,  and  the  final  four  graphing 
programs  which  produce  one  graph  each.  Each  of  the  four  graphs  plot  probability  of 
bit  error  versus  signal-to- noise  ratio  Z  =  1^  PyjT j Nq  in  dB.  Each  graph  has  five  traces. 
One  trace  shows  either  Single  Channel  (SC)  operation  without  Fabry-Perot  (FP) 
filtering  or  SC  operation  with  FP  filtering  and  without  ISI  and  ACI  for  comparison 
with  the  other  four  dense  WDM  traces.  The  other  four  traces  show  probability  of 
bit  error  versus  Z  (dB)  for  four  selected  values  of  the  normalized  channel  spacing 
integer  /,  or  equivalently  the  number  of  adjacent  channels  M.  See  Appendix  C  for 
the  relationship  between  /  and  M.  Each  of  the  four  graphs  is  plotted  for  a  different 
value  I3T,  the  free-spectral  range-bit  period  product.  We  also  note  that  the  four 


21 


graphs  are  generated  with 


p  =  Single  Cavity  Fabry  —  Perot  Filter  Power  Reflectivity  =  0.99 


We  now  present  the  values  of  f3T  the  four  values  of  M,  and  the  corresponding 
values  of  /  used  for  each  of  the  four  graphs.  Each  separate  value  of  M,  or  equivalently 
/,  produces  one  of  the  four  dense  WDM  traces  for  use  in  comparison  with  the  fifth 
trace  showing  single  channel  operation. 

For  Fig.  2:  f3T  —  500  and 


M 


For  Fig.  3:  —  1000  and 


M  = 


For  Fig.  4:  jST  =  1500  and 


M  = 


For  Fig.  5:  =  2000  and 


M  = 


m  ^  ^  ^ 

L  1=4  1=5  1=8  /=20j 


1^  ^  1^  ^ 

L  1=5  1=6  [=9  1=20 


^  m  7^ 

L  1=7  1=9  7=12  1=20 


^  2^  IM  ^ 

L  1=8  1=9  1=12  1=20 


Looking  at  Figs.  2-5,  we  see  that  only  Fig.  2  {/3T  =  500)  is  a  complete  graph. 
For  I3T  =  500  we  were  able  to  compute  all  ten  terms  of  X  (see  Appendix  B)  for  all 
four  given  values  of  M ,  or  equivalently  /,  and  completely  compute  and  simulate  all  ISI 
and  AC  I  effects  of  our  complete  model  given  the  constraints  given  at  the  beginning  of 
this  chapter  {Lq  =  I,  L  =  0,  (jik  =  0,  ujk  =  2'iTkIfT).  For  Figs.  3-5  {I3T  =  1000,  1500, 
and  2000,  respectively),  we  were  unable  to  arrive  at  a  solution  for  Eq.  (24)  (see  Term 
#5  in  Appendix  B),  which  is  the  very  important  K'^  Jq  \sAci\^dt  term.  It  is  very  com¬ 
putationally  intensive  as  it  involves  computing  a  quadruple  summation  multiple  times 
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Figure  3:  Probability  of  bit  error  vs.  Z  (db)  for  —  1000.  Only  7  =  20  trace 
accurate  (within  1/16  dB). 
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Figure  4:  Probability  of  bit  error  vs  Z  (db)  for  =  1500.  Only  /  =  20  trace 
accurate  (within  1/16  dB). 


for  the  various  adjacent  channel  bit  “loadings”  (see  Appendix  B  for  the  details), 
as  well  as  having  to  repeat  the  same  very  lengthy  calculations  for  each  of  the  four 
given  values  of  M,  or  equivalently  I.  However,  we  will  still  be  able  to  draw  some 
conclusions  for  the  /  =  20  traces  of  these  graphs  as  the  channel  frequency  spacing 
between  adjacent  channels  given  by 

=  S  =  ?  = 

where  Rb  is  the  channel  data  rate  or  channel  bit  rate.  Note  that  A/  is  proportional 
to  I.  At  these  larger  frequency  spacings,  ACI  elTects  become  quite  negligible.  Equa¬ 
tion  (37)  explains  why  we  call  I  the  normalized  channel  spacing  integer  as  adjacent 
channels  are  separated  from  each  other  at  an  integer  multiple  I  of  the  bit  rate  Rb- 
Let  us  assume  a  system  design  of  a  filter  free- spectral- range  (3  =  3800  x  10^  Hz  = 
3.8  GHz.  Then  we  can  say  that  for  I3T  =  500  our  bit  rate  Rb  =  7.6  Gb/sec,  for 
/3T  =  1000  our  bit  rate  is  Rb  =  3.8  Gb/sec,  for  j3T  =  1500  our  bit  rate  is  Rb  ~  2.5 
Gb/sec,  and  for  0T  =  2000  our  bit  rate  is  =  1.9  Gb/sec.  It  is  noteworthy  that,  in 
the  time  domain  model,  calculating  probability  of  bit  error  as  presented  in  Chapter 
II  and  Appendix  A  is  easiest  for  the  highest  channel  data  rate  Rb,  or  equivalently  the 
smallest  bit  period  T. 

We  define  power  penalty  as  the  increase  in  signal-to-noise  ratio  required  for  the 
dense  WDM  system  to  achieve  a  10“^^  bit  error  probability  over  the  Single  Channel 
(SC)  system  which  achieves  this  goal  operating  at  Z  =  12  dB.  We  can  now  make 
some  power  penalty  statements. 

Looking  at  Fig.  2  {^T  =  500)  we  see  that  the  power  penalty  increases  as  I 
decreases,  or  equivalently  the  number  of  adjacent  channels  M  increases.  Looking  at 
Eq.  (37),  we  see  that  as  I  decreases,  the  interchannel  frequency  spacing  A /  decreases. 
Thus,  it  is  logical  that  ACI  effects  would  increase.  In  short,  the  channels  become  much 
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more  tightly  “packed”.  Looking  at  the  results  for  Term  #5  in  Appendix  C  when 
/?r  =  500  we  see  that  the  magnitude  of  the  contribution  for  /  =  20  is  only  about  two 
to  four  percent  of  the  contribution  for  /  =  4.  Returning  to  Fig.  2  {(iT  —  500),  we 
see  that  the  power  penalty  for  /  =  20  and  M  =  24  is  approximately  2  1/8  dB;  for 
7  =  8  and  M  =  60,  the  penalty  is  approximately  2  3/8  dB;  for  7  =  5  and  M  —  98,  the 
penalty  is  approximately  2  7/8  dB;  and  finally  for  7  =  4  and  M  =  124  the  penalty  is 
approximately  3  1/2  dB.  It  is  in  this  case  of  7  =  4  and  M  =  124  that  we  see  the  true 
beauty  of  dense  WDM.  If  we  assume  ^  =  3.8  GHz,  then  the  channel  bit  rate  is 
7.6  Gb/sec,  and  the  aggregate  system  bit  rate  Rb  is  124  x  7.6  Gb/s  =  942  Gb/sec, 
which  is  quite  close  to  1  Tb/sec  (1  Tb/sec  =  10^^  b/sec).  Thus,  for  slightly  more  than 
doubling  the  signal-to-noise  ratio  Z  =  R  PyjTjNo,  we  have  increased  the  aggregate 
bit  rate  of  one  single  fiber  close  to  a  terabit  per  second.  Also,  the  dense  WDM 
system  utilizes  a  noncoherent  receiver  and,  thus,  is  much  cheaper  to  implement  than  a 
coherent  optical  FDMA  network  with  expensive  synchronization  circuitry.  Looking  at 
Z,  however,  we  see  that  increasing  signal-to-noise  ratio  is  not  that  easy  to  do.  Upper 
limits  have  already  been  approached  in  making  low-noise  current  amplifiers;  thus,  the 
true  technical  problem  to  be  solved  in  the  dense  WDM  system  utilizing  single-cavity 
Fabry-Perot  filters  is  to  place  an  optical  amplifier  before  the  FP  filter  (see  Fig.  1) 
to  boost  the  received  optical  power  P.  With  this  problem  solved,  economical  dense 
WDM  systems  utilizing  single-cavity  FP  filters  with  aggregate  bit  rates  with  many 
terabits  per  second  are  possible. 

Now,  in  Figs.  3-5  we  have  already  seen  that  at  7  =  20,  the  effects  of  ACI 
become  quite  negligible  as  the  channels  are  spaced  farther  apart  in  frequency.  Thus, 
we  can  conclude  that  the  7  =  20  probability  of  bit  error  traces  are  accurate  for 
=  1000,  1500,  and  2000.  The  rest  of  the  traces  for  7  <  20  are  not.  A  manuscript  in 
preparation  by  Tri  T.  Ha  entitled,  “A  Discrete  Time  Approach  to  Sensitivity  Analysis 
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of  Direct  Detection  Optical  FDMA  Networks  with  OOK  Modulation,”  considers  the 
same  problem  as  this  thesis,  but  uses  a  discrete  time  approach.  The  probability  of 
bit  error  graphs  in  this  future  paper  are  generated  with  the  same  parameters  used  for 
this  thesis.  For  the  discrete  time  approach,  the  N  =  =  1000  and  1500  traces  for 

/  =  20  match  the  /  =  20  traces  for  this  thesis  to  within  less  than  1/16  of  a  dB,  and 
the  N  =  /3T  =  2000  trace  for  /  =  20  matches  the  one  in  this  thesis  to  within  1/5  of 
a  dB. 

Finally,  using  the  case  l3T  =  500  we  feel  that  our  time  domain  model  is  very 
accurate  as  Fig.  2  agrees  with  the  corresponding  discrete  time  graph  for  N  —  500  to 
within  1/4  of  a  dB  for  /  =  4  and  M  =  124,  agrees  to  within  1/8  of  a  dB  for  /  =  5 
and  M  =  98,  and  matches  almost  exactly  for  /  =  8  and  M  =  60,  and  /  =  20  and 
M  =  24.  These  differences  are  well  within  the  realistic  bounds  of  numerical  error  as 
our  time  domain  approach  requires  many  billions  of  calculations. 
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IV.  CONCLUSION 


We  have  presented  a  complete  model  and  equation  for  the  probability  of  bit 
error  of  an  OOK,  dense  WDM  system  employing  single-cavity  Fabry-Perot  filters 
as  channel  demultiplexers.  Our  expression  completely  models  Inter-Symbol  Interfer¬ 
ence  (ISI),  Adjacent  Channel  Interference  (ACI),  and  phase  offset  between  Channel 
0  and  Channel  k  for  all  M  of  the  adjacent  channels.  No  simplifying  mathematical 
assumptions  have  been  made  to  arrive  at  the  final  answer  contained  in  Eq.  (14),  Eq. 
(29),  and  the  sum  of  Eqs.  (22)-(27).  Others  [1,  3]  have  worked  on  this  problem; 
however,  the  solution  is  arrived  at  in  the  frequency  domain  and  often  involves  sim¬ 
plifying  mathematical  assumptions  [1,  Eqs.  (15),  (36),  (39)  and  3,  Eq.  (5)].  The 
closed  form  expression  we  have  presented  for  probability  of  bit  error  for  the  dense 
WDM  system,  although  mathematically  rigorous  and  an  extremely  complete  model, 
has  one  significant  drawback — it  is  extremely  computationally  intensive.  The  single 
complete  graph  for  probability  of  bit  error  when  =  500  took  three  months  to 
compute  using  multiple  SPARC  10  workstations  working  24  hours  per  day.  Even  for 
our  limited  case  oi  Lq  =  L  =  0,  Uk  =  2TrkIjT,  and  (j)k  =  0,  we  never  arrived  at 
solutions  for  the  other  three  values  of  =  1000,  1500,  and  2000.  The  solution  is  to 
use  a  discrete  time  approach.  Again,  we  are  confident  of  the  accuracy  of  our  model 
as  the  same  problem  was  attacked  independently  via  discrete  time  approach,  and  our 
graph  for  =  500  matched  the  discrete  case  of  =  500  to  within  the  bounds 

of  numerical  error.  Finally,  we  again  mention  that  to  practically  implement  dense 
WDM  systems  with  aggregate  bit  rates  of  many  terabits  per  second,  large  increases 
(doubling,  tripling, . . . )  in  the  SNR  of  the  dense  WDM  system  over  that  required  for  a 
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single  channel  system  are  necessary.  To  accomplish  this  an  optical  amplifier  will  have 
to  be  placed  in  front  of  the  Fabry-Perot  filter  to  greatly  increase  the  filter’s  received 
optical  power. 
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APPENDIX  A 

DERIVATION  OF  |s(<)p(it  FOR  USE  IN  X  =  ng  |s(t)pA 

Before  we  begin  the  derivation  of  Jq  \s{t)\‘^dt,  we  must  first  define  the  variables, 
equations,  and  terms  involved: 

Y  —  The  decision  variable  appearing  at  the  output  of  the  integrator  (see  Fig.  1) 
according  to 

Y  =  X  +  N 

X  —  The  deterministic  signal  portion  of  the  decision  variable 

X  =11  r  IsitWdt 
Jo 

7?.|s(f)p  —  Output  current  of  the  photodetector 

1Z  —  Responsivity  of  the  photodetector  (A/W) 

n{t)  —  Postdetection  thermal  noise  with  two-sided  current  spectral  density  Nq  (A^/Hz) 

N  —  Random  signal  component  of  the  decision  variable  Y.  N  is  &  zero-mean,  Gaus¬ 
sian  random  variable  with  variance  NqT 

N  =  f  n(t)dt 
Jo 

p  —  Power  reflectivity  of  the  single-cavity  Fabry-Perot  filter 

/3  —  Free  spectral  range  of  the  single-cavity  Fabry-Perot  filter  (Hz) 

T  —  Data  bit  period  (s) 

P  —  Received  optical  power  (W) 
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w(')  -  {  ol  !therwise  } 

Desired  Channel/ Channel  of  Interest  —  Channel  0 

Adjacent  Channel  —  Channel  k,  A:  =  — M/2,***,— 1,  I,***,  +M/2  where  M  is  an 
even  integer 

6o,o  —  The  detected  bit  of  interest  in  the  interval  0  <  t  <  T 

bo^i  —  Bit  in  Channel  0  during  the  time  interval  [zT,  (z  +  l)r]  where 

bo,i  G  {0, 1} 

bk^i  —  Bit  in  Channel  k  during  the  time  interval  [iT^  {£  +  1)T]  where 

bk,t  e  {0,6^-^*} 

and 

j  =  ^/^ 

(f)k  —  Phase  offset  between  Channel  0  and  Channel  k.  (f>k  is  assumed  to  be  a  uniformly 
distributed  random  variable  between  [0, 27r] 

(l>k  ~  i7[0,27r] 

uik  —  Radian  frequency  spacing  between  Channel  0  and  Channel  k.  Channels  are 
symmetric  around  Channel  0,  i.e., 

LCk  = 

h{t)  —  Impulse  response  of  the  Channel  0  single  cavity  Fabry-Perot  filter 

m  =  (1  -  p)  (( - 1)  (39) 


34 


The  Dirac  Delta  Function  defined  in  two  parts 


t  <  0 
t  >  0 


(40a) 


fO+ 

/  S(t)dt  — 
Jo- 


(40b) 


Complex  baseband  equivalent  data  signal  in  Channel  0 

bo{t)=  E  bo,iPT{t-iT)  (41) 

Lo 

An  integer  greater  than  zero  representing  the  number  of  bits  in  Channel  0  that 
are  trailing  the  detected  bit  6o,o- 

Complex  baseband  equivalent  data  signal  in  Channel  k 

bk{t)=  E  bk,te^^'‘*PT{t  -  iT)  (42) 


An  integer  greater  than  zero  representing  the  number  of  bits  in  Channel  k  which 
trail  the  0*’**  bit  in  Channel  k^  bkfl. 

Received  complex  baseband  signal  at  input  of  the  Channel  0  Fabry-Perot  filter 

+M/2 

^(^)  =  ^/P6o(^)+  E  ^bk{t)  (43) 

fc=-M/2 


We  now  begin  the  derivation  of  Jq  \s{t)\^dt  which  allows  the  computation  of 


=  nr\s(t)f 

Jo 


where  s{t)  is  the  output  of  the  Fabry-Perot  filter  with  an  impulse  response  of  h{t) 


[W,rw7)i^ 
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As  stated  before 


M<)  =  (!-/>)£/<  (i-|) 

+M/2 

r{t)  =  ^bo{t)+  Y1 
k=-M/2 
k:^0 


Substituting  Eqs.  (41)  and  (42)  yields 

Co  +M/2  0  N 

^  bo,iPT{t  -  iT)  +  E  E  h,ee^^^'pTit  -  iT) 

k=-M/2  ' 

A:  9^0 

We  are  interested  in  the  detected  bit  6o,o  during  the  detection  interval  0  < 
we  need  only  evaluate  s{t)  =  ro(t)  in  the  interval  0  <  i  <  T 


/oo  roo 

h{t  —  T)r(T)dT  =  /  r[t  —  T)h(T)di 

-OO  7  —  00 


s{t)  —  ro(t)  =  r{t)  *  h{t) 

=  h{t)  *  r{t) 

/  -\-Ml2 

=  h{t)*\^s/Pbo{t)  +  ^  E  ^k{t) 

k=-M/2 
k^O 

As  convolution  is  a  linear  operator  and  distributes  over  addition 

/  -\-M/2 

s{t)  =  ro{t)  =  iy/p  {h{t)  *  bait))  +  VP  E  *  hii) 

k=-M/2 
k^o 

/  +M/2 

s{t)  =  ro{t)  =  VP  ih{t)  *  bo{t)  +  E  H^)  *  hit) 

k=-M/2 
k^O 


I  (44) 

t  <  T,  so 

(45) 

(46) 

(47) 

(48) 
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Now  looking  at  the  term  within  the  parentheses 

+A//2  0 

h{t)  *bo{t)  +  ^  h{t)  *  bk{t)  =  h{t)  *  ^  bo^iprit  -  iT) 

k=-M/2  •=-£'0 

+M/2  0 

+  E  E 

fc=-M/2 


Applying  the  distributive  property  of  convolution 

+M/2  0 

h{t)*bo{t)+  E  h{t)*bk{t)=  E  h{t)  *  {bo,iPT{t  -  iT)) 

fc=-M/2  •=-'^0 

*5^0 

+M/2  0  /  \ 


+M/2  0  /  \ 

+  E  E  ‘(O »  »*y"“pT(i  -  <r) 

,  ^=-L  \  / 


k^^Mf2 
A:  5^0 


Expressing  convolution  in  terms  of  its  integral  definition 


+Af/2  0  -oo 

;i(t)*6o(0+  E  h{t)*bk{t)=  -  7-)  [6o,,px(r  -  iT)]  dr 

i=-Lo 


k=-M/2 

fc#0 

+M/2  0 


u  -oo 

+  E  E  /  dr 

/I  T  J  —  OO  *• 


k=-Ml2 


Substituting  the  expression  for  h{t  —  r)  [Eq.  (39)]  into  term  Ai,  above 

Ai  =  E  /  (^-p)f^P^^(i-^-^)l^o,iPT(T-iT)dT  (52 

Factoring  (1  —  p)  and  interchanging  the  order  of  integration  and  summation  yields 

Ai  =  (1  -  p)  E  E  I  ^  ~  f  1 

i=-Log=o  \ 
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Using  the  identity:  /_“  x{t)6{t  —  tQ)dt  —  x(to),  the  integral  evaluates  to 


px{t  —  gl l3  —  iT).  We  then  have  the  following  for  Ai 


Q  fOO 

^1  =  “  '^)  [^o,.Pt(t  -  iT)]  dr 

r  J  —  OO 


i=~Lo 


=  (!-/>)!]  U 


i=-Log=0  \  ^  / 

Recalling  the  expression  for 

+^/2  0 

^1=  E  E  /  h{t-T)\bk,ie^->^^PT{r-iT)]dT 

k=^M/2  ^ 


Substituting  the  expression  for  h{t  —  r)  given  by  Eq.  (39)  into  Eq.  (55)  above,  we 


+M/2  0 


^  ~  f  ^  poo  ^ 

Y.  Z  I  (i-rtE- 

/7  r  — OO  n 


k=-M/2 


p9Sit-T-^]  \h,(ei'^>‘^PT{r  -  ^r)  dr  (56) 


Factoring  1  —  p,  interchanging  the  order  of  integration  and  summation,  and  factoring 
bk/  and  p^  from  the  integral 

+M/2  0  OO  /-OO  /  a\- 

Bi  =  {l-p)  E  E  /  8  it  -  T  - pt{t  -  iT)dT  (57) 

k=-Ml2  ^ 


Applying  the  property  x{t)S{t  —  tQ)dt  =  x(fo) 


e^‘"*>T(r  -  iT)dT  =  | 


Substituting  Eq.  (58)  into  Eq.  (57)  yields 

+A^/2  0  OO  ,  ,  , 

Bi  =  Y  h{t-T)  -  iT)]  dr 

0—  T  — OO  ■' 

A:^0 

+  M/2  0  CO  ( 


+M/2  0  CO  / 

(1  -  /^)  E  EE  f 

l=-Lg=0  \ 


k=-M/2 
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Substituting  Eq.  (54)  and  Eq.  (59)  into  Eq.  (51)  and  then  substituting  this  result 
into  Eq.  (48)  yields 


s{t) 


ro(t)  =  y/P 


0  CO 

(i-rt  E 

i=— Lo  g’-O 


Ai 


+M/2  0  oo  f  n  \ 

»./n  e=-L9=0  \  P  / 


k=-M/2 
k:^0 


Bi 

We  can  now  factor  out  (1  —  /?)  from  our  expression  above.  Then 

0  oo  /  ^ 


s{t)  =  ro{t)  =  y/P  {I  -  p) 


T^P^bo,iPT 

Li=-Log=0  \  P  ) 


+M/2  0  oo  (  n  \ 

A:  9^0 


(60) 


(61) 


Looking  at  Eq.  (61).  We  separate  the  terms  involving  the  detected  bit  of  interest  in 
Channel  0,  6o,o,  the  trailing  bits  in  Channel  0,  6o,,'  *  7^  0,  and  the  terms  involving  bits 
in  the  other  channels,  bk^t-,  A;  0 


s{t)  =  rQ{t)  =  y/P  [I  -  p) 


K 


^  fPbo,oPT 
*b(*) 


+  51  £  P^bo,iPT  6  -  f  - 

i=^Lo  0=0  \  P  / 


-Lo  g 

SlSj(i) 
-\-Mf2  0  CO 


+M/2  0  00  /  n  ^ 

+  E  EE  <  -  I  - 

«,-L,=o  V  /5  / 


A:=-M/2 
it^O 


5AC/(*) 


(62) 
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where  sb(0  is  the  desired  signal,  sisi{t)  is  the  Intersymbol  Interference  (ISI)  signal, 
and  SAci{t)  is  the  Adjacent  Channel  Interference  (ACI)  signal.  Using  this  compact 
notation,  we  have 

s{t)  =  ro{t)  =  K  +  Sisi{t)  +  SAci{t))  (63) 

Sketching  SB{t)  for  6o,o  =  1 


We  can  also  sketch  the  intersymbol  interference  due  to  the  trailing  pulses  in  Channel 

0 


♦ 


ISI  contributions  to  the  pulse  of  interest 


Now  let  us  assume  in  one  pulse  interval  T  that 


T 

1  =  9 
0 


(64) 
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where  1/^  divides  evenly  into  the  period  T  making  q  an  integer  quantity.  We  are  only 
interested  in  the  detection  interval  0  <  i  <  T,  so  we  only  need  to  evaluate  5(<)  =  ro(<) 
in  the  interval  0  <  t  <T.  Looking  at  the  desired  signal,  which  is  the  system  response 
to  the  bit  of  interest  bo  o 


We  can  see  that,  during  the  interval  0  <  t  <  T,  the  response  is  a  sum  of  scaled, 

shifted  rectangular  pulses.  Therefore 

% 

g-1 

^B{t)  =  X)  P^K'oPT 

5=0 

We  now  direct  out  efforts  towards  developing  an  expression  for  S[si{t)  during  the 
interval  0  <  t  <  T.  We  begin  by  examining  the  ISI  term  for  t  =  —1  in  Eq.  (62), 
denoted  by  —IIS I. 

-  1  75/  =  £  p^bo,-iPT  +  for  -  T  <  f  <  oo  (67) 

g=0  \  P  ) 
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By  examing  the  sketch  of  —IIS I 


we  see  that  —1/5/  will  also  be  a  sum  of  shifted,  scaled  rectangular  pulses  for  0  <  t  < 
T.  The  -1/5/  expression  will  have  two  parts,  the  first  of  which  is 

=  Y1  p'K-\Pt  >  for  0  <  <  <  T  (68) 


But  with  the  q'-'^  shift  we  no  longer  sum  scaled  pulses  beginning  at  t  =  0  and  end¬ 
ing  at  1/^,2/^,-  •  •,  but  now  we  sum  pulses  which  progressively  scale  and  begin  at 
2/^,  •  •  •,  and  end  at  t  =  T.  This  second  part  of  the  expression  looks  much  like 
the  expression  for  5^(0  during  the  interval  0  <  t  <  T  [Eq.  (66)].  Thus,  we  may  write 
the  second  part  of  —1/5/ 


+  z2  P^^o-iPt 


9-Q 


MTjj  ■ 

/ 


for  0  <  i  <  T  (70) 


Examine  the  ISI  term  for  t  =  — 2  in  Eq.  (62),  —21  SI 


2  75/  =  ^  bo-2P^PT  ( ^  ^ 


We  sketch  —21  SI 


-2T<t<oo  (71) 


-2T  -2T  +  -  -T 

r' 


■2r+- 


-2T  + 


(9-1) 


-2T  + 


(29-1) 


-T  + 


-2T  + 


(39-1) 


By  analyzing  this  sketch,  we  can  obtain 

-  2  75/(1)  =  . 


for  0  <  i  <  T  (72) 


-2/57(2)  =  Zp^^o,-2Pt 


I _ 

9-29\U  V  ^ 


for  0  <  <  <  r  (73) 


Note:  These  rectangular  pulses  shift  into  the  interval  0  <  t  <  T  on  the  2o‘*'  shift. 


Similar  to  —I/S'/: 


we  may  write  directly 


-3  ISl  =  p^bo.^pT 

g=2q+i 


(mt)\ 

\9  -  2q) 


4g-i 

+  P^K-sPt 

g=3q 


for  0  <  t  <  r  (76) 


Then  from  Eqs.  (70),  (74),  (76),  the  expressions  for  the  -I/S'/,  -2/S'/,  and  -3ISI 


pulses,  we  are  ready  to  use  the  pattern  developed  to  write  an  expression  for  the  ISI 


between  0  <t  <T.  Recalling 


sisi{t)  =  Z)  D  P^f>o,iPT  \  , 

i=-Lo  5=0  \  P  J 

we  may  write  directly  from  the  pattern  developed 


for  -LoT  <t<oo  (77) 


..,(0=  i;  'e" 

i=-Lo  s=-(14-t)g+l  Vy  +  ll+O? 


-(i-I)g-l 

+  Z  P^^0,iPT 


g-\riq 


g  +  iq 


for  0  <  f  <  r  (78) 


We  now  have  expressions  for  5b(<)  and  s/s/(0)  during  the  interval  0  <  t  <  T.  We 
have  only  the  ACI  term  Sy[ci{i)  to  calculate 

4'M/2  0  oo  /  17  \ 

SAciit)  =  Z  Z  Z  iov  -LT<t<oo 

k=-M/2  ^ _ t 

k^O  * 


*  Note:  This  portion  of  SAci{t)  is  almost  exactly  like  the  ISI  term  with  which  we 
just  worked  [see  Eq.  (62)].  The  only  difference  is  the  complex  exponential  factor. 
Looking  at  the  0''*'  term  of  the  ACI  in  the  k'‘^  channel 


OACI  =  Z -  1]  > 


for  0  <  f  <  oo  (80) 


we  sketch  OACI 


We  see,  that  in  a  manner  similar  to  the  SB{t)  signal,  the  pulses  add  as  before,  but 
now  the  amplitude  of  the  pulses  is  modulated  by  p3ej<^k(t-{g/P))  versus  just  before. 
Then 


5=0 


for  0  <  t  <T  (81) 


The  —lACI  will  be  similar  to  the  —IISI  term  [see  Eq.  (70)].  However,  now  we  not 
only  have  p^  but  also  the  complex  exponential  gated  by  the  pulses  in  the  summation. 

-1  AC  I  = 


T _ 

9-q 


9-9 

/? 


,  for  0  <  t  <  T 


The  derivation  of  —2  AC  I  and  —ZAC  I  will  proceed  in  a  similar  manner  as  —21  SI 
and  —31  SI  [Eqs.  (74)  and  (76)].  We  see  that  two  different  type  pulses  gate  the 
product  of  p^  and  a  complex  exponential  term.  These  gated  terms  are  then  summed 
over  the  period  0  to  T.  We  then  use  our  expressions  for  OACI,  —lACI,  ...  to  write 
an  expression  for  SACi{t)  during  the  interval  0  <  f  <  T  in  a  manner  similar  to  the 
•s/s/(0  expression.  We  may  now  write  the  final  expression  for  the  ACI 


SAClii)  = 


\MI2 

p-1 

E 

E 

^=-M/2 

lg=0 

A:^0 

\ 

-1 

i-^qyi 

+  i: 

E 

^=-L 

s=-(l+^)g+l 

mt) 

9  +  {l+£)q 


9=-ig 
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PT 


T- 


9  +  k 


t  - 


'nil 


/ 


for  0  <  t  <  T 


(83) 


/ 


We  have  now  derived  S[si{t),  and  SACi{t)  during  the  interval  0  <  t  <  T.  We 

now  develop  the  expressions  needed  to  compute  |s(t)p  and  \s(t)\^dt.  Substituting 
the  expressions  for  5s(t),  sisi{t),  and  SACi{t)  for  the  interval  0  <  t  <  T  into  Eq.  (63) 
yields  s{t)  in  the  interval  Q  <t  <T 

s{t)  =  VP  —  p)[sB{t)  +  ^/5/(f)  +  SAci{i)]i  ior  0  <t  <T  (84) 

'  K  Iqq  ' 

where  ^^(t),  sisii^)i  ^•nd  SAcii^)  ^re  described  by  Eqs.  (66),  (78),  and  (83)  respec¬ 
tively.  Then 

|s(t)p  =  Kh'{t)s'{ty  (85) 

Dropping  the  (t)  notation  for  convenience  in  the  three  terms  of  s'{t)  and  noting  that 
complex  conjugation  is  a  linear  operator 

s{t)*  —  Sg  +  s}sj  -b  Saci  (86) 


So 

s'{t)s'{t)*  =  {sb  +  SJSI  +  SaCi){s*b  +  -Sj5j  +  •s)4C/) 
Multiplication  yields 

s'{t)s'{t)*  =  SbSq  +  SBS*jgj[  -t-  SbS*aci 

+SjsiSq  +  Sisis}si  +  SlSI^*ACI 


Psaci^b  +  •Sac/'S/S/  +  ^aci^aci 


Rearranging  and  simplifying  yields 

|s(i)|"  = /^V(t)5'(t)*  =  K^x 


SbS*b  +  SISIS*JSI  +  SaCIS*acI 


hsP  |s/s/P 


(87) 


(88) 
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+  ■SB'S/5/  +  -S/SZ-Sb 

■■ - V - ' 

2Re[sBs}sj]  =  2Re[sisiSB] 

+  ^B^*ACI  +  ^ACI^B 

' - V - ' 

2Re[s*BSAci]  —  2i?e[sB5^c'j] 

+  ^ISI^ACI  “I"  ^ACI^ISI 
" - ^ 

2Re[s},<;jSAci]  =  2Re[sisis\(;<i]  (89) 


We  will  integrate  each  of  these  terms  over  0  to  T  to  compute  |s(t)pc?t.  To  avoid 
confusion,  we  must  realize  that  constant  must  be  carried  through  to  the  final 
calculation  of  all  quantities.  Recalling  Eq.  (66) 

/  \ 


g-l 

=  P^^OfiPT 
5=0 


T 


T-l 

\\ 


K 


) 


We  know  6o,»'  £  (0,1),  p^  is  real,  and  pr(-)  is  also  real.  Then 


(90) 


SB  =  &B 


(91) 


and 


/  / 

q-l 

Y1  /^o.oPt 

g=o 


\ 


\V 


rp  J 


/ 


Factoring  6o,o  it  is  a  constant  to  the  summation,  we  obtain 


(92) 


/ 


|sb|  —  6o,o 


/ 


g-l 

J^P^PT 

5=0 


V 


\ 


T  -  ^ 

\\  /?> 


9 

t  -  - 


\ 


(93) 
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Performing  the  squaring  operation  yields 


/ 


I  o  2  l2 

I'SS  - 


{ 


q-l  q-1 
5=0  m=0 


V 


\\  h 


t-i- 


Pt 


rr.  rn' 
T - 

vv 


'  m 

K~l, 


(94) 


Then 


K^\sb?  =  P{l-p)X.o 


5-1  5-1 

'EE 

5=0  m=0 


W-i. 

w 


Pt 


Tfl 


'  m  ^ 

y~lj 


/J 


for  0  <  t  <T  (95) 


Now  we  want  to  compute  Jq  K^\sB\^dt  [see  Eq.  (89)]. 

r  K^\sB?dt  =  P(1-p)XoEEp’'"”x--- 

5=0  771=0 


f 


PT 


T-l 

VV 


o 

t  _ 


PT 


(  m ' 
'  T _ 

VV 


^  m  ^ 


dt 


(96) 


Looking  at  Eq.  (96),  we  have  interchanged  the  order  of  integration  and  summation. 
As  the  pulse  functions  are  the  only  functions  dependent  on  time  in  the  expression, 
we  look  at  the  product  of  the  two  unit  amplitude  pulse  functions  in  rj 
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I  I  =4'  Product  is  0 
Y/}  ==^  Product  is  1 


We  have  arbitrarily  assumed  that  m  >  but  our  logic  would  also  work  if  we  assumed 
g  >  m.  The  product  of  two  pulses  is  1  for  m//?  ^  t  <  T  and  0  otherwise.  Thus,  we 
see  that  the  integral  of  the  product  of  the  two  pulse  functions  in  rj  will  be  a  square 
area  of  height  1  x  {T-  largest  of  gl(3  or  m/^).  Thus 

[K^Usfdt  =  P(1  -  g  '")'))  ,  for  0  <  (  <  r 

^=0m=0  \  \  P  // 

(97) 

where  max(5r,m)  is  defined  by 

max (xi, 0:2):  Choose  largest  of  x\  or  X2,  which  are  both  positive. 

If  Xi  =  X2,  then  max(xi,X2)  =  xi  =  X2. 

We  now  need  to  compute  K‘^\sisiYdt.  We  see  from  Eq.  (78)  that  s/5/  is  real  and 
not  complex  since  feo,t  €  {0,1}.  Thus 


sisi  —  s)si  ^■nd  |5/s/|^  = 


(98) 


k/5/P  = 


-1  (-'7)-l 

r  E 

i=-Lo  s=-(l  +  i)?+l 


I^KiPT 


a 


(  pm  \ 

\9  +  (1  +i)q) 
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To  perform  the  squaring  we  set  up  another  version  oi  a  +  b  with  the  ^^-index  changed 
to  m  and  the  Tindex  changed  to  r  to  account  for  cross-product  terms.  We  call  this 
new  representation  c  +  d.  Thus 


1 5/5/ 1^  =  (fl  -f-  d)  — 

A  B  C  V 

Performing  the  multiplications  yields 

-1  (-w)-i  -1 

A  =  ac  =  ^  m  X) 

i=—LQ  (l+t)g+l  r=— Lo  m='-(l+r)94-l 


X 


Pt 


(3T{t) 


g  +  {l+i)q)^^  \m  +  {l  +  r)q^ 


mt) 


(100) 


B  =  ad  = 


-1  (-S'?)-!  -1  -(r-l)ij-l 

E  E  EE  p‘*’"Ki\r 

i=—Lo  ^=— (l-|-i)qi+l  t^  —  Lq  m——rq 


(  mt)  \ 

+  + 


V 


T- 


m  +  rq 

r~r~. 


I 


C  =  bc 


-1  -(t-l)g-l  -1  (-r?)-! 

E  E  E  E  X 

i=~Lo  g~—iq  r=— Lq  m=  — (l+r)g+l 


it  -  + 


Pt 


V 


9  +  n 


Pt 


mt) 


m  +  (1  +  r)q^ 


) 


V=bd  = 


-1  -(t-l)g-l  -1  -(r-l)?-l 

E  E  E  E  /■'"ioA.rX 

i=— Lo  g~—iq  t=  —  Lq  m=—rq 


Pt 


Tit 


'9  +  iq\\\ 

.  ^  J) 


9  +  k 


Pt 


Tit- 


m  +  f'qW  \ 


Tn-\-  rq 

n~, 


(101) 


(102) 


(103) 
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Integrating  over  the  interval  of  interest:  lO\sisi\^dt  • 

I<^  \sisi\‘^dt  =  Cdt  +  £  Vdt  (104) 

The  equations  for  A,  B,  C,  and  V  are  Eqs.  (100)-(103),  respectively.  We  now  analyze 
term  A  [Eq.  (100)]  by  sketching  the  multiplication  of  the  pulses 

p,  (_mt)  ] 

y  \m  +  il+r)q) 

Product  is  0 

— _  i  YyA  Product  is  1 


^(l  +  f)7  m+(l+r)9 


pT{t) 


\5  +  (1  +i)q) 

Again,  we  have  arbitrarily  assumed  one  pulse  lasts  longer  than  the  other.  Since  the 
height  of  each  scaled  pulse  is  1,  we  see  that  the  area  under  the  product  is  1  times  the 


minimum  of 


q+{l  +  i)q  Tn  +  {l-{-r)q 

— , — -  or  - - — 

p 


We  can  apply  the  same  logic  as  used  before  to  arrive  at  the  value  of  the  term 


IC  /  Adt 


Then 


rl  -1  -1  (-Tq)-l 

I  Adt  =  Pii  -  pf  E  E  E  VO 

t=:-Lo  5f  =  ~(l+t)(7+l  r  =  -Lo  m=-(l+r)q+l 


(105) 
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where 


V?o  — 


min  (gr  +  (1  +  i)q,  m  +  (1  +  r)q) 


where 


min  (^  +  (1  +  i)q^  m  +  (1  +  r)q)  is  defined  by 

min(a:i,i2)*  Choose  smallest  of  Xj  or  X2,  which  are  both  positive. 

If  Xi  =  X2,  then  min(xi,X2)  =  Xi  =  X2. 

We  now  compute  Vdt.  Looking  at  the  expression  for  V  [Eq.  (103)],  we  sketch 

the  product  of  the  pulses,  arbitrarily  assuming  m  rq  >  g  iq 


I  I  =^>  Product  is  0 
Y/A  ==r'  Product  is  1 


g  +  tq  m  +  rq 


We  see  the  product  of  pulses  only  exists  from  the  maximum  of 

g  -Yiq  m-\-rq 

(3  ~r~ 

Thus,  the  integral  of  the  pfoduct  of  the  pulses  will  be 

m&x{g  +  iq,  m  +  rq) 


Then 


-1  -1  -(r-l)g-l 


K^l  Vdt  =  E  E  E  p’*’"KAr 

®  i=—Lo  g=—iq  r=—Lo  m=:—rq 


max(5  +  iq,  m  +  rq) 


(106) 
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giving 

-1  -1  -(i-l)g-l 

e  =  E  E  E  E  /■"’"W.  x--- 

r=  — Lo  m=— (l4-7')g+l  9=~'^^ 


(108) 


If  we  look  at  the  expression  above  Eq.  (108),  we  see  that  B  is  identical  to  C  [Eq. 
(102)]  in  every  way  except  the  order  in  which  the  summations  appear.  The  order  of 
the  summations  can  be  rearranged  because  the  terms  inside  of  the  summations  are 
completely  separable  in  relation  to  the  sets  {m,r},  {2,5)}.  Thus,  we  can  rearrange 
the  order  of  the  summations  and  conclude 


B  =  C  (109) 

and  that 

f  Bdt  +  K^  [  Cdt  =  2K^  [  Bdt  =  2K^[  Cdt  (110) 

Jo  Jo  Jo  Jo 

We  choose  to  integrate  twice  the  value  of  B  [Eq.  (101)].  Interchanging  the  order  of 

integration  and  summation  yields 

y  /  -1  -1  -(r-l)g-l 

2K^J  Bdt  =  2K^  I  E  E  E 

®  \»=— Lo  5=— (l+i)g+l  r=—Lo  m=—rq 
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(Ill) 


m)  \ 

g-V{l-\-i)q) 


PT 


( 


\ 


m  + 

~}) 


m-\-  rq 

~T~. 


\ 


} 


dt 


Sketching  the  product  of  the  two  pulses,  we  see  the  product  will  exist  for  m  +  rq  < 
g  +  {I  +  i)q  and  will  be  zero  otherwise. 


Product  is  0 
Product  is  1 


We  create  a  “gating”  function  to  account  for  this 


mt)  ] 

g  +  {l  +  i)q) 


Pt 


Gi{g,i,Tn,r)  = 


Then  finally 

2K^  r  Bdt  = 

Jo 


I  ’  0, 


for  m  -f  <  gf  +  (1  -f-  i)q 
otherwise 


(112) 


/(■“  r  Bdt+K^  r  cdt 

Jo  Jo 

-1  -1  -(r-l)7-l 

E  E  E  /+“VA.'G'i(s,i. "!,>■) 

t=— Lo  5=— (l+Oq+l  r=— Lo  m=— rg 

(113) 
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Finally,  combining  the  results  of  Eqs.  (105),  (106),  and  (113)  we  have 

[  sisis}sidt  =  [  \sisi\^dt  =  [  Adt+[  Bdtf  C  dt  +  [  Vdt 

Jo  Jo  Jo  Jo  Jo  Jo 

(114) 

where 

K^J  Adt  =  P{l-pf  ■£  E  E  E  V’o 

^  i^z—Lo  ^r=  — (l-j-t)g+l  r~~Lo  m=-*(l+7*)<?+l 

in  which 

90  =  +  +  ■>."■  + (1  +  ’;)^)^ 

r  Bdt  +  K^  r  Cdt  =  2K^  r  Bdt 

Jo  Jo  Jo 

-1  (-w)-l  -1  -(r~l)g-l 

=  2P{\-I,f  E  E  E  E  p’*’"bo,bo.rG,(g,i,m,r) 

t=-Lo  5=-(l4-t)g4-l  r=~Lo  m=-rg 


and 


Vdt  =  p(i-pf  Y.  E  E  E  /'""V.-V.  X 

^  —  n  —  —  tn  r*— *_r.«  TT7~— Tn 


T- 


izz—Lq  g=—iq  r=:— Lo  m^z—Tq 

'  max(^  +  iq,m  +  rq) ' 

0  ^  , 


Now  we  must  compute  |5^c;/p  and  Jq  \sACl\^dt.  We  see  that  we  have  a  complex 
bk,e  as  bk,i  G  {0,  as  well  as  a  complex  exponential  [see  Eq.  (83)]. 

We  begin  by  noting  that  the  complex  conjugate  operator  has  the  following  properties 

[uv)*  =  u*  V* 

{u  +  v)*  =  u*  +  u* 


Taking  the  complex  conjugate  of  Eq.  (83) 


MC/ 


+A//2 

E 

k=-M/2 

k^O 


I 


\ 


q-1 


•g=o 


T  -  - 


t  - 


Id. 


fd. 
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Letting 


and 


-1 

+  E 

e=-L 


.3=-(H-^),+l 


(  mt)  \ 

Vfl'  +  (1  +  (^)q) 


9^-iq 


ii(33)' 
(■-  W))  I 


c 


c 


f 


d 


d' 


\g  +  {l+i)q) 

(  mt)  \ 

^  \m  +  (1  +  T)q] 


Pt 


Pt 


K-m)) 


T- 


m  +  rq 

rj~. 


Then,  with  a  change  of  indices,  we  can  express  as 


MC/ 


+M/2  g^l 
m=0 

n=— M/2 


n^O 


A* 


(115) 

(116) 

(117) 

(118) 

(119) 

(120) 

(121) 
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-hM/2  -1  (-rq)-l 

+  E  E  E 

n=-M/2  ’•=-^'^=-(1+09+1 


+M/2  _l  -(r-l)g-l 

+  E  E  E 

n>f/o  m—-Tq 

71=  — M/2 

n^O 


(122) 


+M/2  (j  — 1 

k=-M/2 

k:^0 


+M/2  -1  (-«g)-l 

+  E  E  E 

k=-M/2  ^=-i'9=-(i+e)q+i 


+M/2  _i  -(^-l)g-l 

+  E  E  E 


k--M/2 

k^O 


-L  9=-iq 


(123) 


Then 


k.io/P  =  sjc;4c;  =  (A  +  B  +  C)  (A-  +  B-  +  C’) 


(124) 


Now,  distribution  of  multiplication  over  addition  yields 


AA*  +  AB*  +  AC*  +  BA*  +  BB*  +  BC’ 


+CA*  +  CB*  +  CC’ 


(125) 
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Rearranging  yields 


=  AA*  +  BB*  +  CC*  +  2  -(AB*  +  BA* 


+2  [-  (AC*  +  CA*)]  +  2  (BC*  +  CB’ 
.2  .  .2 


Using  the  identity  Re[Z]  —  {Z  +  Z*)/2  we  have 


AA*  +  BB*  +  CC*  +  2i?e[AC* 


(126) 


+2Re  [AB*]  +  2i?e[BC* 


(127) 


Multiplying  the  terms  above  yields 

H-M/2  q  —  1  -{-M/2  q  —  l 

AA*  -  y,  Yl  Yl  V 

k=-M/2  n=-M/2 


(128) 


+M/2  -1  +A//2  -1  (-'■?)-l 

BB-=  E  E  E  E  E  E  7 

k^-M/2  <=-‘s=-(l+»)«+l  „=-M/2  '=-!■  »i=-(l+')«+l 


(129) 


where 


+M/2  -1  -(f-l)g-l  +M/2  -1  -(r-l),-l 

E  E  E  E  E  E  c 

l=-M/2  '=-‘  »=-'•  .=-M/2 


(130) 


where 


To  compute  2Re  [AC*],  we  know  Re[-]  is  a  linear  operator  so  it  may  be  moved  across 


the  summations.  After  multiplication  of  terms  A  and  C*  we  get 

+M/2  5-1  +M/2  -1  -(r-l)9-l 

2iie(AC1  =  2  E  E  E  E  E  « 

te-M/2  '="  «.-M/2 


(131) 
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where 


We  know  b^fi,  6* ,.,  and  are  complex.  We  will  calculate  k  in 

Eq.  (131)  using  the  identity 

MZ]  =  i(Z  +  Z-) 

Since  (uv)*  —  u*v*,  by  extension  it  can  be  shown  that 

{uvwY  =  u‘v*w* 

Also  note  that  if  S'  =  Uj  +  a2  +  cts  +  04  H —  ■  where  Gi,  02,  ^3,  G4,  • ' '  complex,  then 
i^efS”]  =  Re[a\]  +  Re[a2\  +  i?e[a3]  +  Re[aY[  +  •  •  • 

We  will  also  use  the  fact  that  if  z  is  complex  and  a  is  constant,  then 

Re[az]  =  aRe[z\ 


Using  the  facts  above,  we  can  see 

L-  r=  J{Wk-'^n)t-Wk{glP)->ruin(ml(3)\ 

p^+'^ad' 

Then 

p^+'^ad' 

with  a  and  d'  defined  by  Eqs.  (116)  and  (121),  respectively. 


(132) 


(133) 


Now,  we  will  turn  to  computing  2Re  [AB*].  This  will  have  five  summations  just  like 
the  2Re  [AC*]  term  above. 

+M/2  g-l  +Af/2  -1  (-rg)-l 

2i?e[AB*]=  E  E  E  E  12  ^  (^34) 

k=-M/2  n=-M/2  r=-Lm=-{l+T)g+l 

n^O 

By  inspection  of  Eq.  (123) 

After  letting  the  1/2  generated  when  applying  Re[-]  operator  to  the  argument  cancel 
the  2  outside  of  the  summations,  we  get 

—  [0*:,o0n,r® 

+^fc  o6„  p^'^'”ac'  (135) 

With  a  and  c'  defined  by  Eqs.  (116)  and  (119),  respectively. 

We  now  compute  the  term  2i?e[BC*].  This  term  will  have  six  summations.  Moving 
the  Re[]  operator  inside  of  the  six  summations  to  the  argument  yields  the  following 
by  inspection 

t  =  Re  [bk,tK,  ^^30) 

After  letting  the  1/2  generated  when  applying  the  Re[-]  operator  in  i  [Eq.  (136)] 
cancel  the  2  outside  of  the  summations,  we  can  then  write  down  the  expression  for 
2i?e[BC*] 

+M/2  -1  (-^g)-l  +M/2  -1  -(r-l)g-l 

2JJe(BC-l  =  E  E  E  E  E  E  ‘  (137) 

fc-M/2  (l+')«+l  ,=_„/2 

n^o 

I  =  +  blfbn,r  X  •  •  • 

g-iK^fc -Wn)t-Wk (g//?)+u/n (m/p)] j  ^+m 


with  c  and  d'  defined  by  Eqs.  (118)  and  (121). 
Now  we  must  find 


I  K^\sAci\^dt  =  [  AA*dt  + r  BB*dt  + r  CG*dt 

Jo  Jo  Jo  Jo 

+K^  ^  2Re  [AC*]  dt  +  2Re  [AB*]  dt 

+K^  2Re[BC*]dt  (138) 

where  we  recall  that 

=  P(1  _  py 


In  the  term  A  A*  [Eq.  (128)],  we  have  the  pulse  multiplication  aa'  where 


a.  =  PT 


a'  =  pT 


\ 


) 


Arbitrarily  assuming  m  >  g,vre  sketch  the  product  a  x  a'. 


Product  is  0 
Product  is  1 


The  product  is  1  and  exists  from  the  maximum  of  {mfp,g//3)  to  T.  Either  g  or  m 
can  be  the  largest.  Hence,  we  will  be  integrating 

^  =  ed(‘*’fc-Wn)i-i*//t(3//3)+Wn(m//3)]  (139) 

with  the  lower  limit  of  integration  being  {lfl3)ma,x{g,m).  The  upper  limit  of  inte¬ 
gration  is  T.  We  consider  two  cases  for  Eq.  (128). 
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uik  = 


Case  1:  k  ^  n  u)k  7^ 
When  k  ^  n,  ujk  ^  oJn  because 


(140) 


where  again  is  the  radian  frequency  spacing  between  Channel  0  and  Channel  k, 
and  Au>  is  the  uniform  radian  frequency  spacing  between  adjacent  channels. 


J.  +M/2  q-l  +M/2  q-l 

[  AA*dt  =  Pil-pf  ^  ^  X 

JO 


k=-Ml2  n=-M/2  ° 


i 


(l//3)max(5,m) 


k^O 

idt 


n^O 


(141) 


Letting 


X  =  j[{uJk-‘^n)t-UJkig/^)  +‘^n{m/l3)] 

dx  =  j[{uJk-U!n)]dt 


then 


I 


(l//?)max(5,m) 


Cdt  = 


J[(uJk-(^n)t-UJk(9/P)+<^n{m//3)] 


j{uJk  -UJn) 


^  Li[(wjc-ui„)r-Wfc(g//?)+UJ„(m//3)] 

j{uJk  -UJn)  L 


(l//3)max(g,m) 


(142) 


Case  2:  k  =  n  =>■  ojk  = 

Again  similar  logic  to  that  used  above  [Eq.  (140)]  allows  us  to  conclude  OJk  =  Ct;„.  For 
this  case,  the  complex  exponential  ^  reduces  to 


^  _  gi[0-i*'fc(s//3)+w„(m//3)] 


(143) 
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Now,  since  ujk  =  oJn-,  we  can  further  reduce  the  complex  exponential  to 


(144) 


Therefore,  for  Uk  ~ 


/  ^dt  -  /  dt  =  r  -  -  max(5,  m' 

'  {1  / P)  nia.x{g ,m)  J {1  / p)  max(g .m)  (j 


Thus,  we  may  write  the  expression  for  lO  AA*dt 


+M/2  q-1  +M/2  q-1 


le  I  AA-dt  =  P{l-p)^  E  E  E  Ea>’^”X'^i  (1«) 

*^0  - n  _ n 


k=-M/2  ^  °  n=-M/2 


where 


J[((^k-i^n)T-u>k{9/P)-^<^n{m/P)] 


_^i((‘^fc-‘^n)(l//3)max(5,m)-a>/k(5//3)+cc-n(m/)3)] 


for  k  ^ 


Now  we  turn  to  the  BB*  term  [Eq.  (129)].  We  have  the  pulse  multiplication  cc'  where 


c  =  pt 


mt) 

g  +  (1  +  ^)9 


c  =pt 


m  +  (1  +  r)q 


Arbitrarily  assuming  m  +  (1  +  r)q  >  +  (1  +  £)q,  we  sketch  the  product  c  x  c'. 
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) 

□ 


Product 

Product 


is  0 
is  1 


The  product  exists  and  is  1  from  0  to  minimum  of 

fg  +  {l+C)q  m  +  (l+r)g\ 

[  P  '  13  J 

Hence,  we  will  be  integrating  the  same  complex  exponential  (  [Eq.  (139)]  as  for  the 
AA*  term  [Eq.  (128)],  but  the  upper  limit  of  integration  will  be 

i  nun{g  +  (1  +  %,  m  +  (1  +  r)q) 

while  the  lower  limit  is  0.  We  again  consider  two  cases. 


Case  1;  k  ^  n  ujk  ^  qJn 

After  interchanging  the  order  of  integration  and  summation,  we  will  be  integrating 
[see  Eq.  (129)] 


(147) 


Then 


/■(1/0)  min(g+(l+e)g,m+(l+r)q) 

/  ^dt 

Jo 
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|(1/^)  min(£?+(l+^)<?,m-H(l+r)g) 


j{uk  -  LOn) 


j{uJk  -UJn) 


pj[(<^k-<^n)t-Uk(g/P)+UJn{m/P)] 


J[{uJk-^n){l/P)  rmn(g+(W)q,m+{l-\-r)q)-ujk{g/P)-{-uJn{'m/p)] 


_Jl~^k{9/P)+^n{m/P)] 


(148) 


Then,  factoring  the  result  yields 


r{l/P)  min(5+(l+£)g,m+(H-r)9) 


(dt  = 


j{uJk  -tOn) 


J[-‘^k(9/0)+<^n{m/0)]  gi[{wfc-u;n)(l//3)min(s+(l+^)?,m+(l+r)q))  _  (149) 


Case  2:  k  =  n  =>  =  uin 

This  is  similar  to  the  A  A*  term  for  uik  =  <^n  [see  Eq.  (144)] 


^  _  QM^^k/P){,m-g)] 


(150) 


/  (dt  =  e^l(‘^»=/^)(’"-s)l  /  dt 

Jo  Jo 

^  ^j[{^k/mm-9)]  }-mm{g  +  {l+e)q,m  +  {l+r)q)-Q  (151) 


r(l//3)min(£f-l-(l+^)(?,m+(l+r)g) 


Thus,  we  may  write  the  expression  for  Jq  BB*dt 

.j'  +M/2  —1  (— ^g)-l  4-Af/2  — l  (—rg)  — 1 

m-<n  =  p{\-f,f  E  E  E  E  E  E 


fc=-A//2 

ki^O 


(152) 
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where 


h,eK,r 


j{uJk  -^n) 


7[(u;fc-wn )(!//?)  mm(p-f-(l+^)g,m+(l+ 


09)]_l)  ^ 


for  k  ^  n 


=  ^  min(5r  +  (1  +  £)q,  m  +  (1  +  r)q)  J  , 

for  k  =  n 

\  / 
Now  we  integrate  the  CC*  term  [Eq.  (130)].  We  have  the  pulse  multiplication  dd^ 


where 


d'  =  pt 


m  +  rq 


d-„  \  \  ^  JJ  d'-u  V  V  ^ 

Looking  at  the  pulse  multiplication  graphically  where  we  arbitrarily  assume  m  +  rq  > 


g  +  (q  for  the  purpose  of  sketching  the  situation,  we  have 


9  +  ^q 


m  +  rq 

~r~ 


□  =»  Product  is  0 
\//\  Product  is  1 


Looking  at  the  equation  for  CC*  [Eq.  (130)],  we  see  that  after  interchanging  the  order 
of  integration  and  summation,  we  will  be  integrating  the  same  complex  exponential 
as  for  the  AA*  and  BB*  terms 


—  pi[(‘^fc-Wn)<-Wfc(g//3)+a;„(m//0)] 


4  =  e- 


(153) 
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where  the  lower  limit  of  integration  will  be 


1  /  . 


while  the  upper  limit  will  be  T.  Then  for 
Case  1:  k  ^  n  =>  LOk  ^ 

rT  1 

(dt  =  - - 

(l//?)max(5+^g,m+rg)  j[^k  ^  ^n) 


I 


1 


j{uJk  -  W„) 


j[(u}k-OJn)t-iUk  (5//5)  +  t^n (m/P)] 

\(l/P)max(g-\~£q,m-\-rq) 
J[i<^k-’^n)T-u>k{g/P)-{-u>n{m/p)] 


_^jl{^k-^n){l/P)  Tnax{g-\-iq,m-Jt-rq)-uJk{g/P)-\-ujn{Tn/P)] 


(154) 


Case  2:  k  =  n 

Again,  similar  to  the  AA*  term,  for  k  =  n,  the  term 


—  pjK^k-<^n)t-uik{g/P)+w„{m/0)] 


^  =  e 


(155) 


reduces  to 


_  pM^k/0)im-g)] 


(  =  e' 


(156) 


and  then 


r  ^dt  ^  r  dt 

J(l/P)ma.x(Q-\-iQ,Tn-\-rQ)  J(l 


—  gi[(t^fc//3)(™-s)] 


'(l/P)  max{g+lq,m+rq) 

r  -  ^  max(5  +iq,m  +  rq) 


(157) 


Thus  we  may  write  the  expression  for  Jq  CC*dt 


+M/2  -1  +M/2  -1  -(T-l)q-l 


a:' 


2111  ^  Zl  1111  ' 

J  cc-dt  =  p{i-py  E  E  E  E  E  E 

^  L_  g~-iq  _ ,,,,  r=-L  m--rq 


A:=-M/2 
ki:0 


n=-M/2 
n^i^O 


X  </?3 


(158) 
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</?3  -  < 


^k,tK,r 


j{u}k  -  UJ„) 


Jl(ujk  -Wn  )T-(Ajk  (5//3)+u/n  (m//3)] 


^^jl(^k-i^n)(l/^)max(g+iq,m-hrq)~ufk(9/j3)+i^n(m//3)] 


for  A:  ^  n  > 


1 


T  -  -  ma.x{g  +  iq,Tn  +  rq) 


ior  k  =  n 


Now  we  turn  to  the  2Re  [AC*]  term  [Eq.  (131)].  We  have  the  pulse  multiplication  of 
ad'  where 


a  =  pj 


V 


/, 


d'  =  pf 


N  } 


H 


m  +  rg\'\\ 


Looking  at  this  graphically  by  arbitrarily  assuming  m  +  rq  >  g,yve  sketch  the  product 


a  X  d' 


Product  is  0 
Product  is  1 


Then  we  can  see  that  the  product  ad'  will  exist  from  (1  / /3)  max(p,  m  +  rg)  to  T,  which 
are  the  lower  and  upper  limits  of  integration,  respectively. 

Case  1:  k  ^  n  =>  ^  uJn 

This  expression  has  the  two  complex  exponentials 


( 


„jUuk-w„)t-w,,(g/P)+u)„(m/p)] 


(t 


Q  -  e 


(160) 


to  integrate.  Letting 


u  =  -j  -  uj^)t  -  LOk  {gjP)  +  (m//?)] 


d\l  =  -j{LOk  -  Un) 


the  integration  yields 


'{l/P)  max(5,m+rg) 


[-j{uk-iJn)] 


'-j{uJk-UJn)] 


(1//?)  max(^, m+rg) 


‘j[{^k-^n)T-LUk{9/P)+UJn  (m/P)] 


:,-j[{uJk-^n)(l/P)nia.x{g,Tn-\-rq)-uJk(9/P)+^n{Tn/P)] 


(161) 


l{l/P)ma.x(g,m‘^rq) 


(dt  =  - - -  c 


j[(u>k-iOn)T-ujk{g/P)-^u;n{m/p)] 


j{uJk  -UJn)  ^ 

^p3[(^k-^n)(l/P)  m8Lx(g,m-\-Tq)-u>k  (9/P)+^n{m/P)] 


(162) 


Case  2:  k  =  n  =>  =  Un 

Looking  at  the  expression  for  2Re  [AC*]  [Eq.  (131)],  we  see  that  for  cu*.  =  the  two 
complex  exponentials  reduce  to  and  which  are  constants 

to  the  integration.  Thus,  we  will  integrate  the  product  of  the  pulse  functions  ad'  from 
0  to  T.  This  is  the  only  time  function  for  this  case.  Then,  looking  at  the  sketch  below 
Eq.  (158) 

rT  pT  2 

/  ad'dt—  /  dt  =  r  —  —  max(o,  m  +  ro)  (163) 

Jo  d(l/0)max(fl,m+rg)  ^ 


(163) 


Thus,  we  are  ready  to  write  the  expression  for 


r 2Re[kG*]. 
Jq 
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Then, 


-fM/a  q^l  +M/2  -1  -(r-l)g-l 


A'"  /  2Ae[AC-l*  =  />(I-rt“  E  E  E  E  E  ('<54) 

3=0  r=-L  )n=-rg 


where 


<^4  =  < 


P^^^h,ob:,r 


j{u;k  -tOn) 


k=-M/2  n=-M/2 

k^O  719^:0 


7  [(wfc  -  Wn )  T- Offc  (p //?) + u/n  (  m //?)] 


oi[(t^A:-t*^n)(l//?)  max(5,m+rg)-u;fc(5//?)+a;n(m//?)] 


I  g+mL>K  ^  [  1  L-i[(wfc-w„)r-a^fc(p//?)+a;„(m//?)] 

I  P  /:,0  n,r  [ [~j(a;fc-arn)] 


^p-j[{^k-^n){l/P)  inax(g,m+rg)-u;>.(g//?)+ct»n(m//3)] 


,  for  k  ^  n 


r  -  i  max(5',  m  +  rq) 


ioT  k  =  n 


We  now  turn  to  the  2i?e  [AB*]  term  [Eq.  (134)].  We  have  the  pulse  multiplication 
ac'  where 

/  \ 


a  =  pr 


^  rp  _  £_ 

VV  0 


9  \  \ 


c  =  pr 


/3T{t) 


^m  +  {l+r)q^ 

To  form  a  sketch,  we  arbitrarily  assume  m  +  (1  +  r)q  >  g,  yielding 


/ 
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I  I  Product  is  0 
Y//\  ==»  Product  is  1 


We  see  that  the  product  will  exist  for  m  +  (l  d-r)?  >  g  and  will  be  0  otherwise.  Thus, 
unless  m  +  (l  +  r)q  >  g,  the  integral  of  2Re  [AB*]  is  0.  We  will  later  create  a  “gating” 
function  to  account  for  this. 


Case  1:  k  ^  n  ^ 

For  this  case  where  m  +  (1  +  r)^  >  g  we  will  be ‘integrating  the  same  complex 
exponentials  as  before 


^  ^  ^jl{^k-^n)t-iOk{9/PH^n(m/l3)] 


(165) 


and 


g  —  ^-M>^k-‘^n)t-ut;{g/0)+w„(m/l3)] 


(166) 


We  have  computed  these  integrals  before  and  will  just  substitute  the  new  upper  and 
lower  limits  in  the  final  answer  for  this  term 

IT  T  •  (m  +  (l  +  r)9) 

Upper  Limit  :  ^ - — 

13 

Lower  Limit  :  4 

/? 

If  m  +  (1  +  r)^  <  g,  the  integrals  of  the  two  complex  exponentials  are  zero  as  the 
product  of  the  pulse  functions  ac'  -  0.  Again,  we  will  create  a  “gating”  function  to 
allow  for  this. 
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Case  2:  k  —  n  ==^  LOk  =  ijJn 

Looking  at  the  2Re  [AB*]  term  [Eq.  (134)].  We  see  that  the  two  complex  exponentials 
reduce  to 


These  are  constants  to  the  integration  with  respect  to  t.  Thus  as  before,  we  will  only 
be  integrating  the  product  of  the  pulse  functions  ac'  from  0  to  T.  The  integral  has 
value  if  m  +  (1  +  r)q  >  g  and  is  0  otherwise.  Then,  if  m  +  (1  +  r)q  >  g,  the  limits  of 
integration  are  the  same  as  for  the  case  Uk  /  ujn  and 

rT  /•(m+(l+r)g)//?  1 

/  ac'dt=  /  dt  =  -[{m  +  (1  +  r)q)  -  g]  (167) 

Jo  Jg/P  H 

Now  we  will  create  the  “gating”  function  to  null  the  integrals  if  the  product  ac'  is  0 
from  0  to  r. 

To  indicate  the  zeroing/nulling  of  all  integrals  if  >  m  +  (1  +  r)q,  we  will  create 
another  gating  function  called  G2{gtm,r). 

(  1,  for  <  m  +  (1  +  r)q  ) 


G2ig,m,r)  =  < 


(168) 


[  0,  otherwise 


We  are  now  ready  to  write  the  expression  for 

J.  +M/2  9-1 

2Re[kB*]dt  =  P{1- pf  Y.  E 

k=-M/2 
k^O 


2Re  [AB*]  dt.  Then 


-l-M/2  -1  (-rg)-l 

E  E  i: 


n^O 


<^5 


(169) 
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where 


( O^'^'^bunh*  _ _  \pi[{'^k-<^n){{m-\-(\+r)q)/l3)-uJk{glP)+uJn(ml(3)\ 

\r  K,U  n,T  j(uik-OJn)  L 


oi[(‘^n/^)(m-g)] 


J-  ( /)S+’TiA»  fj  _ 1 _  \f,-jl(<^k-<^n)({m+{l+r)q)/p)-Wk{g/l3)+u:„{m/0)] 


_e  ili<^n/p)(m  g)]  j  G2{g,m^r),  for  A:  ^ 


;[(m  +  (1  +r)5)  -5]  G'2(5,m,r), 


for  A;  =  n 


where 


02(9, m,r] 


1,  for  5  <  m  +  (1  +  r)^ 


0,  otherwise 


Now,  finally  we  must  write  out  the  last  (sixth)  term  of  the  ACL  We  now  need  to  find 


r  2Re  [BC* 
Jo 


Again,  we  look  at  the  expression  for  2Re  [BC*]  [Eq.  (137)]  to  see  the  pulse  multipli¬ 
cation  cd'  where 


c  =  pt 


mt) 

+  (1  +  £)q 


d'  =  pT 


T{t- 


m  +  rq 

.~T~ 

m  +  rq 


We  sketch  at  this  pulse  multiplication  by  arbitrarily  assuming  H-  (1  -f-  tjq  >  m-\-rq. 
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Product  is  0 
Product  is  1 


Case  1:  k  ^  n  =>  Uk  /  u;„ 

Here  we  will  be  integrating  the  complex  exponentials 


(170) 


and 


(171) 


The  integrals  will  have  the  following  limits 


Upper  : 
Lower  : 


g  +  (1  -ki)q 


m  +  rq 

~T~ 


The  integrals  will  have  value  for  m  +  rg  <  g  +  (l  +  £)q  and  will  be  zero  otherwise. 
We  have  computed  these  integrals  before  and  will  just  substitute  the  new  limits  in 
the  final  answer  for  this  term.  We  will  create  another  “gating”  function  to  account 
for  the  integrals  “turning  off”  when  m  +  rg>g  +  (l  +  i)q 

[1,  m  +  rq  <  g  +  {1  +  £)q  \ 

G3{9,£,m,r)  =  I  >  (172) 

(  0,  otherwise  J 

Case  2:  k  =  n  ==^  Uk  = 
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Looking  at  the  2Re  [BC*]  term  [Eq.  (137)],  we  see  again  that  the  complex  exponentials 
reduce  to 

which  are  constants  to  the  integration  with  respect  to  t.  Thus,  in  a  similar  manner 
as  before,  we  will  only  be  integrating  the  product  of  the  pulse  functions  cd'  from  0 
to  T.  The  integral  has  value  if  m  +  <  p  +  (1  +  i)q  and  is  zero  otherwise.  We  have 

already  created  the  “gating”  function  [Eq.  (172)]  to  account  for  this  turning  on  and 
off  of  the  integral.  Then,  ii  m  +  rq  <  ^'  +  (1+  £)q  the  limits  of  integration  are  the 
same  as  for  the  case  ^  c<;„  and 

/  cd  dt=  dt  -  -[{g  +  {1  +  £)q)  -  {m  +  rq)]  (173) 

^0  J  {m-\-rq)/p  jJ 

We  are  now  ready  to  write  the  expression  for 

r  2Re[BC*]dt. 


Then 

j.  +A//2 

-1 

+M/2 

-1  ~(r-~l)g-l 

/  2Re  [BC*]  dt  =  P(1  -  pf  Y. 

E  E 

E 

E  E 

J  0 

e=-L  g=-{l+e)q+l 

n=-M/2 

tz=:^L  m  =—rq 

A:  9^0 

n^O 

(174) 
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where 


a+mL  U*  1  \pj[{'^k-<^n){(g  +  {l+^)<l)/P)—‘^k(9/0)+'^n{m,/P)] 

P  ^k,lOn^r  j[ojk-uJn)  r 


_gi[(wfc-Wn)((m+r9)//3)-a/ife(a//J)+w„(m//3)] 


I  -a+mL*  L  _ _  L-j[(u/k-Wn)((p+(l+^)9)//3)-c^fc(fl//3)+w„(m/^)] 

+P  [ 

_e-i[K-‘^n){(m+ra)//3)-wfc(3//3)+t^n(m//?)]  ^  771,  r),  for  A:  7^ 


for  k  —  n 


-[{g  +  (1  +  i)q)  -  (m  +  rq)]  Gsig,  £,  m,  r) 


where  Gz{g,i,m,r)  is  defined  by  Eq.  (172).  We  have  now  completed  integrating  all 
six  terms  of  the  ACL 

Recalling  the  expression  for  \s{t)\'^  [Eq.  (89)],  we  see  that  we  still  have  to  compute 


the  following  terms 


2i?e[5B'S/5/]  =  2Re[sisiSB]i 
2Re[sBSACi]  =  ‘2Re[sBS*y^cIh 


2Re[s}siSACi]  =  2.Re[sjsis*Aci]j 

multiply  each  by  the  constant  =  P(1  —  p)^  and  integrate  each  from  0  to  T.  We  will 
use  either  form  of  the  above  three  expressions  depending  on  the  ease  of  computation. 
Recall  from  Eq.  (66)  that 


for  0  <  t  <  r  (175) 


Let  us  look  at  sbs}^^  term.  Looking  at  the  expression  for  SB(t)  during  the  interval 
0  <t  <T  [Eqs.  (66)  or  (175)],  we  see  that  sb(A)  is  real  as  bo,i  6  {0,1}.  Looking  at 
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the  expression  for  s/s/(0  during  the  interval  0  <  t  <  T  [Eq.  (78)],  we  see  also  that 
•5/s/(t)  is  real  as  again  bo^i  £  {0, 1}.  Then  Sfsi  —  s*jgj  and  sss^gj  =  sb>s/5/,  and 


2ile[sBs}5j]  =  2sbsisi 


(176) 


Performing  the  multiplication  yields 

2Re[sBS*fgj]  =  2[sb5/5/] 


9-1  -1  (-‘9)-i 

=  2i:  E  E  /■"“VoV.w 

g=Qi=-Lo  m=-(l+i)q+l 


/?)  ^(m  +  (H-!)g 


9-1  -1  -(t-l)9-l 

+2E  E  E  /■""'VoV.PT 

g=Oi=-Lo  m=-iq 


T  t- 


T  -i. 


m  +  iq 

.~T~^ 

m  +  iq 


(177) 


We  see  that  when  we  multiply  by  and  integrate,  the  integral  will  distribute  across 
addition  to  the  two  terms  e  and  f.  The  integrals  then  move  inside  of  the  summations 
to  the  pulse  function  products.  Now  we  need  to  figure  the  proper  limits  for  the 
integration  of  the  two  pulse  function  products  in  e  and  f  [Eq.  (177)].  We  note 
that  the  products  of  pulse  functions  are  the  only  time  functions  in  the  expression 
2Re[sBS*isi\i  so  they  are  the  only  integrands  in  the  expression.  We  examine  the  pulse 
multiplication  in  term  e  [Eq.  (177)].  Let 


E  =  pt 


r 


0m 

m  +  (1  +  i)q 


(178) 


Sketching  E,  we  arbitrarily  assume  m  +  (1  +  i)^  >  g. 
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I  I  =^>  Product  is  0 
Y/\  Product  is  1 


The  integral  Edt  will  have  value  for  <  m  +  (1  +  i)q  and  will  be  zero  otherwise. 
If  <  m  +  (1  + 

fT  r{m+{i+i)q)/P  J 

i  "  Lf  * " •*’"  "^  **  ■*■  “  ’I 

Since  only  the  pulse  function  product  is  being  integrated,  we  will  create  another 
gating  function,  which  also  gives  us  the  value  of  the  above  integral  when  it  turns  on. 
This  is  similar  to  Gi{g,i^m,r)  [Eq.  (112)]. 

T  f  4  [("^  +  (1  +  0?)  -  S']  >  5  <  m  +  (1  +  z)g  1 

Edt  =  G4{g,i,m)  =  I  ^  I  (180) 

{  0,  otherwise  J 

Look  at  the  pulse  multiplication  in  term  f  [Eq.  (177)].  Let 


Sketching  F,  we  arbitrarily  assume  m  +  iq  >  g. 


(181) 


79 


We  see  that 

/•^  1 

Jo  ^  ^  ~  ^  max(5,  m  +  iq)  (182) 

We  can  now  write  down  the  following  integral 

tT 

y  2Re[sBS*jgj]dt 

Substituting  integrals  Jq  Edt  and  JoFdt  for  the  pulse  function  products  in  the 
2/2e[sBsJ^j]  expression  [Eq.  (177)],  and  multiplying  by  yields  the  following  result. 

Note: 

tT 

Edt  is  substituted  for  the  pulse  function  product  in  e  in  Eq.  (177). 

[T 

F  dt  is  substituted  for  the  pulse  function  product  in  f  in  Eq.  (177). 
Recalling  that  lO  =  P(1  —  we  have 

/  2Re[sBS*i^j]dt  =  f  2[sBSisi]dt 
Jo  Jo 

=  2P{1  -  pYYJ,  m  P^'*'”'f>o,obo,iG4{g,i,m) 

g=0  i=—Lo  m=— (l+j)9+l 

9-1  -1  -(i-l)9-l  r  ,  -1 

+  2P(1  -  />)^  X!  H  -  -a  inax(5,  m  +  iq)  (183) 

g=Oi=  —  Lo  P 
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G 


q-l  +M/2  -1 

+  E  E  E  E 

^=°  k=-Ml2  ‘^=-^rn=-{\+l)q+\ 
k^Q 


H 


q-\  +M/2  -1  -(^-1)9-1 

+  E  E  EE 

k=-Ml2 

A:^0 
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Pt 


Pt 


■V- 

I 


\ 


/ 


(184) 


The  pulse  products  are  G,  H,  and  I  as  shown  above.  Then  we  note  again  that  6o,,'  is 
real.  Also,  G  {0,e^'^''}  and  are  complex.  We  then  have 


2Re[sgSACi]  =  2i?e[sBS^c/]  = 

9-1  +M/2  9-1 

k=-M/2  '”=° 

A:  9^0 


9-1  +M/2  -1  (-^9)-l 

+E  E  E  E  (/'""Vo 


3—^  k-^M/2  ^=“(1+^)'?+! 


9-1  +M/2  -1  -(^-1)9-1 

+  E  E  EE  (/"""Vo 


®=°  k=-M/2 
A:^0 


(185) 

We  see  that  computing  Jq  2Re[s*QSACi]  will  involve  integrating  each  of  the  six 


terms  inside  of  the  three  clusters  of  summations.  The  integrable  part  of  each  of  these 
six  terms  is  a  complex  exponential  of  the  form 


(186) 


or 

fj.  =  (187) 

The  pulse  products  G,  H,  and  I  will  provide  the  limits  of  integration. 
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Looking  at  G  [Eq.  (184)] 


G  =  pr 


Tit-i 


T  t- 


T  -  ^ 


(188) 


Sketching  G,  we  arbitrarily  assume  m  >  g. 


Product  is  0 
Product  is  1 


i.  m. 


Since  the  height  of  the  pulse  product  G  is  1  from  (1//3)  max(p,  m)  to  T,  we  see  that 
we  will  have  the  following  limits  of  integration  for  the  complex  exponentials 


(189) 


(190) 


which  are  the  integrands  in  the  first  cluster  of  summations  in  the  2Re\s*QSAci\  term 
[Eq.  (185)1 

Upper  Limit  :  T 


Lower  Limit  :  —  max(p,  m) 


Turning  to  H  [Eq.  (184)] 


H  =  pt 


rp 


pm 

m  +  (1  +£)q 


(191) 
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Sketching  II,  we  arbitrarily  assume  m  +  {I  £)q  >  g. 


I 


Product  is  0 
Product  is  1 


We  see  that  for  5'  <  m  +  (1  +  i)q,  the  upper  and  lower  limits  of  integration  for  the 
integrands  t?  =  and  g.  =  contained  in  the  second  cluster  of 

summations  in  the  2Re[sgSAci]  term  [Eq.  (185)]  will  be 


Upper  Limit  :  ^(m  +  (1  +  i)q) 

Lower  Limit  :  ^ 

H 


We  will  now  create  another  “gating”  function  to  null  the  integrals  of 

and  fi  =  e  Jn  the  second  cluster  of  summations  in  the  2i?e[sgs>ic/]  term 

[Eq.  (185)].  This  nulling  of  the  two  integrals  will  occur  if  S'  >  m  +  (1  -f  £)q. 

So  we  define 


f 

GsigJ^m)  =  { 


9  <Tn  +  {l+£)q  I 


(  0,  otherwise  J 

Now  we  compute  the  integrals  of  t?  =  ^  _  g-ia/*(<-(m//?)) 

upper  and  lower  limits  of  integration  71  and  ai 


(192) 


for  arbitrary 


•/ai 


— 1 . 


71 

<*1 


(193) 
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Similarly  for  /u 


r  fidt  =  r 

Jai  J  a\ 


71 

ai 


Now  we  turn  to  the  pulse  product  term  I  in  Eq  (184) 


Sketching  I,  we  arbitrarily  assume  m  +  £q  >  g. 


(194) 


(195) 


We  see  from  the  sketch  that  the  upper  and  lower  limits  of  integration  of  the  integrands 

^  _  g-iwfc(t-(m//3))  jjjg  third  cluster  of  summations  in  the 
2Re[s*gSACi]  term  (Eq.  (185)]  are 


Upper  Limit  :  T 

Lower  Limit  :  —  ma.x{g,  m  -t-  £q) 

Now  that  we  have  the  three  sets  of  limits  of  integration  for  each  of  the  clusters  of 
summations  and  the  indefinite  integrals  of  and  /x,  we  are  prepared  to  write  the 
expression  for 

(  2Re\s*gSACi\dt  =  f  2Re[sgSAci]dt 

Jo  Jo 
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Then,  recalling  that  —  P {I  —  pY 


where 

in  which 

(p7 


in  which 

¥^8  = 

where 


lRe[sgSACi]di  =  f  2Re[sBSAci]dt  =  Aba  +  Bba  +  Cba  (196) 

Jo 


9-1  +M/2  9-1 

ABA  =  P(l-l>fY,  E  Ev>7 

0  =  0  ,  772=0 


2^k 


g-jc^,(T-(m//3)) 


9-1  +M/2  -1  (-€9)-1 

Bba  =  P{1-p)'‘T.  E  E  E  V>8 

A:=-M/2 

A:^0 


J<^k 


giwfc((l+^)9/^)  _  gjwitCis-m)//?) 


05(9, 


,-J‘^fc((l+^)?//3)  _  p-i<^k{.(a-m)IP) 


m) 


^5(5^ )  tti) 


<  m  +  (1  + 
otherwise 


9-1  +M/2  -1  -(^-1)9-1 

C'ba  = -P(l  - H 

fc=-M/2 

fc#0 
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in  which 


_  n3+'^h  h,  _  iu-fc((max(s,m+^g)-m)//?) 

mg  —  n  OQ  QOk  t  .  c  e 

_L  „3+m7  /«  ^  \f,-}f^k(T-{ml0))  -iwk{{ma.x{g,m-\-lq)-m)l  P) 

The  last  term  we  have  to  compute  is  K'^  Jq  2Re[s*jsjSAci]dt.  We  have  already  shown 
that  Sjsi  is  real;  thus 

2Re[s*jsiSAci]  =  ‘2‘Re[sisiSAci]  (197) 

As  before,  we  will  multiply  s/5/  and  saci- 

Look  at  sisi  [Eq.  (78)],  we  will  let  the  first  term  be  A  and  the  second  be  B.  Looking 
at  Saci  [Eq.  (83)],  we  will  let  the  first,  second,  and  third  terms  be  C,  D,  and  E, 
respectively.  Then 

sjsisaci  —  (.^  "t"  +  Z?  +  E)  =  AC  4"  AD  +  AE  +  BC  +  BD  +  BE  (198) 


Then,  we  compute  all  six  of  the  above  terms  individually 
-1  +M/2  g-1 

=  E  E  EE  X  •  •  ■ 

l=:Lo  5=:-(l+t)g  +  l  ]z~-Ml2 


{  mt)  ] 

\g  +  il  +  i)q) 


(199) 


-1  +M/2  -1  i-eq)-I 

AD=  ^  E  EE  E 

i=— Lo  ij=—(l+t)g4*l  k=—M/2  77i=  —  (  1+^)9 +1 

ki^O 


mt) 

g  +  {l  +  i)q 


mt) 

m  +  (1  +^)q 


(200) 


(-I?)-!  +A//2 

AE=  Y1  E  E  E  E 

z=-Lo  ^=-(14-09+1  k=-M/2  m=-iq 

k^O 


„  (  'in 

""  U  +  (1  +  .)J  "" 


V 


'  m-\-  Iq 

n~, 


/ 


(201) 


-1  -(t'-l)9-l  +M/2  ,-l 

BC=Y,T.  E  E X  •  •  ■ 


i=-Lo  a=-iq  rn=0 


k^O 


Pt 


f--(-(^)) 

(-(’-7^))j 


\  /_/  m\\ 

^  '-W 


Pt 


m 


Ml 


(202) 


-1  -(t-l)9-l  +M/2  -1  (-^?)-l 

BD=  H  EE  E  X  .  ■  ■ 

in  — Lo  9=—i(}  A:=— M/2 


Pt 


(33) 

^  ^  yy 


\ 


Pt 


(  mt)  \ 

\m  +  (1+  t)qj 


(203) 


-1  -(«'-l)9-l  +M/2  -1  -(^-1)9-1 

BE=  Y.  Y.  E  E  E 

i=-Lo  s=-.</  jt=-M/2 


Pt 


-(T))] 


'■■(•1(33)' 

VV  V  /3  Jj/ 


(r-( 


(204) 
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When  we  distribute  2/?e[  ]  across  each  of  the  terms  AC  +  AD  H - and  use  Re[Z\  = 

[Z  +  Z*]^,  the  2  will  cancel  with  the  1/2,  and  we  will  be  integrating  the  same  complex 
exponential  terms  and  fi  =  as  we  did  before  for  the 

2Re[sgSAci]  term  [see  Eqs.  (186)  and  (187)].  Then,  all  we  need  to  do  is  investigate 

the  pulse  products  for  each  of  the  six  terms  AC  +  AD  H - to  figure  out  the  proper 

limits  of  integration.  Also,  we  may  need  to  create  other  “gating”  functions  to  null  the 
integrals  if  there  are  specific  index  occurences  which  make  the  specific  pulse  product 
0. 


The  AC  pulse  product  is 


PT 


(  mt)  \ 

\9  +  il+i)q) 


Arbitrarily  assuming  ^  +  (1  +  i)^  >  m,  we  sketch  the'  pulse  product 


I  I  Product  is  0 

Y/A  Product  is  1 


We  see  that  the  limits  of  integration  will  be 


Upper  Limit : 
Lower  Limit  ; 


^(5  +  (1  +  i)q) 

m 

J 
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Now  we  define  the  appropriate  “gating”  function 


Geig.i.m) 


<  5  +  (1  +  i)q 


0,  otherwise 


(205) 


This  nulls  the  two  applicable  integrals  when  the  pulses  do  not  overlap. 

The  AD  pulse  product  is 

\5  +  (l+09/  Vm  +  (l  +  ^)9y 

Arbitrarily  assuming  m  +  (1  +  >  y  +  (1+  i)q^  we  sketch  the  pulse  product. 


□  Product  is  0 

V/\  Product  is  1 


g  +  (1  +  i)q  m  +  (l  +  £)9 

P  0 


We  see  that  the  limits  of  integration  will  be 


Upper  Limit  :  -  min(5  +  (1  +  i)q,  m  +  (1  +  £)q) 

Lower  Limit  :  0 


The  AE  pulse  product  is 


^  U  +  (i  +  0?/^^ 


'  m  +  £q 
m  +  £q'' 

“T“. 
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Arbitrarily  assuming  g  +  {1  +  i)q  >  m  +  iq,  we  sketch  the  pulse  product. 


I  I  =>  Product  is  0 
Y/X  Product  is  1 


m  +  iq  9  +  {l  +  i)q 


We  see  that  the  limits  of  integration  will  be 


Upper  Limit  : 
Lower  Limit  : 


^{g {I  X- i)q) 
^(m  +  (q) 


Now  we  define  an  appropriate  gating  function 


G7{g,i,£,  m) 


1,  m  +  £q  <  g  +  {1  +  i)q 
0,  otherwise 


This  nulls  the  two  applicable  integrals  when  the  pulses  do  not  overlap. 


The  BC  pulse  product  is 


Arbitrarily  assuming  g  -{■  iq  >  rn,  we  sketch  the  pulse  product. 


(206) 
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1 


[  I  =»  Product  is  0 
y/A  =>  Product  is  1 


m  g  +  iq 

P  P 

We  see  that  the  limits  of  integration  will  be 


Upper  Limit  :  T 


Lower  Limit  :  max(5  +  iq,  m) 


The  BD  pulse  product  is 

/ 

PT 


T[t- 


'9  +  iq\\\ 


P 


V 


T- 


9  +  iq 

.  p  . 


Pt 


PT{t) 


^m  +  {l+e)q^ 


} 


Arbitrarily  assuming  +  iqr  <  m  +  (1  +  l)q^  we  sketch  the  pulse  product. 


Product  is  0 
Product  is  1 
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We  will  now  apply  the  methods  outlined  below  Eq.  (204)  and  perform  the  2i2e[-] 
operation  on  each  of  the  six  terms  of  s*[gjSAci  —  ^isi^ACi-  We  will  then  integrate  the 
resulting  i?  =  and  ^  each  of  the  six  clusters  of 

summations  and  multiply  by  the  proper  gating  function,  if  necessary,  to  compute 

f  2Re[s}gjSAci]dt  =  f  2i?e[s/5/s^c/]dt 
Jo  Jo 

Recalling  that  =  P (I  —  pY 
T  T 

[  2Re[s*jgjSAci]d'l^  =  f  ‘^RG[sisi^ACi]dt  =  ACR+ADn+AEfi+BCfi+BDfi+BEfi 

JO  JO 

where 

-1  (-«g)-l  +M/2  q-l 

ACr  =  P{1-p?  EE  E  Evio  (208) 

z=-Lo  5=-(H-t)q4-l  k^^M/2 
k:^0 

in  which 

V^io  =  P^^'^io,ii>k,o~ — 

J<^k 

+/'^"io.-6:o7-^  g-,w.[(l//3)(3+(l+.)g)-m//3]  _  1  G^{g,i,m) 

’  ’  [-J^k]  [  ^ 

where 

'1,  m  <  g  +  {\+i)q 

Gs{g,i,m)  =  < 

0,  otherwise 

> 

*  Note:  e^^'^AmlP-mip)  ^  ^jo  ^  ^ 

-1  +M/2  (-^g)-l 

ADn  =  P(l-pYJ2  E  E  E  E  V’U  (209) 

izz-Lo  p=_(l+{)g+l  m=~(l+^)g+l 

k^O 

in  which 


C^ll  =  p9+mfj  .bk,e^  ■ 


Geig,i,m) 
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,  0+mt  L*  ^  r„-i‘^fc[(l/)3)min(s+(l+«)g,m+(l+€)g)-m//3]  _  ^-iwfc[-m//3] 


-1  (-«9)-l  +A^/2  -1  -(^-l)g-l 

AE^^Pii-pfY,  E  E  E  E 

i=-Lo  g=-(l+i)q+l  k=-M/2  ”i=-fg 


(210) 


in  which 


<^*12  —  P^^'^bo^ibk/-.  6' 


,i<*'lc[(l//9){3+(l+«)9)-™//5]  _  pj<^k[k/P] 


+/9®+'"6o<^fc^ — T —  L-jw*[(i//3)(fl+(i+09)-W/J]  _  ^7(5-,^,  A  ”2) 

’  [-iwfc]  ^ 


where 


G7{g,iJ,m)  = 


1,  m  +  iq  <  g  +  {I  +  i)q 
0,  otherwise 


in  which 


-1  -(«-i)g-i  +M/2  q-l 

BCk=P(1-pY  E  E  E  Evis 

i=-Lo  g=-iq  k=-M/2  '"=° 


(211) 


o+mL  L  LjW[T-(m//?)]  _  ^ju;fc[(l//?)max(5+tg,m)-m//3] 


g-jwfc[T-(m//?)]  _  g-jw/c((l/^)max(g+i?,m)-m//?] 


-1  -(t-l)g-l  +M/2  -1  (-^g)-l 

BDR  =  p{i-pf  EE  E  E  E  <^14  (212) 

t=-I/0  3=-*9  k=-M/2  ^=-£'"1=-(1+^)9+1 


in  which 


SP14  0,.  L  J 

+p®+'"6o .-it  L-j-fc[(i//5)(i+^)?]  _  e-J-*[(i/^)(3+’?)— /^]]  GsC^,  ij,m) 

’  '  [-M]  ^ 
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where 


Gs{g,i,l,m)  = 


1,  g iq  <  m tjq 

0,  otherwise 


-1  -(i-l),-!  +M/2  -1 

BEn  =  P{l-pf  EE  E  E  E  915 


,=-Lo  3=-.? 


in  which 


(213) 


—  P^'^K.ih.e- —  e 


Jujk[T-{7n/^)]  __  ^ju;k[il/P)Tna.x(g+iq,Tni‘lq)--m/l3] 


+/+”6o,i6: 


[T-(m//?)]  _  ^-juJk[(l/P)msiM9+iq,Tn-\-£q)-m/P] 


We  have  now  computed  all  terms  of 


f  \s{t)\'^dt  =  f  s\t)s\tydt  [see  Eq.  (89)]  (214) 

Jo  Jo 


Now  we  recall  from  Eq.  (89)  that 


/  \s{t)\‘^dt  =  [  s\t)s\ty dt  =  [  \s\t)\^dt 

Jo  Jo  Jo 

=  [  \sB\^dt  +  (  \sisi\^dt  +  f  \sAci\'^dt 

Jo  Jo  Jo 

+  /  2Re[sBs}si]dt  +  f  2Re[sBSAci]dt  +  K'^  /  [s}sjSACi]dt  (215) 

Jo  Jo  Jo 

We  will  now  extract  the  answers  for  the  above  six  integrals  from  the  previous  rather 


lengthy  derivation  so  that  we  will  have  the  final  answer  for 


y\s{t)\^dt=^K^  y\s'{t)f 

Jo  Jo 


(216) 


in  one  single  place.  Now,  before  we  consolidate  the  terms  we  once  again  restate 

max(xi,X2):  Choose  largest  of  xi  or  X2,  which  are  both  positive. 

If  Xi  =  X2,  then  max(a;i,X2)  =  ti  =  X2. 


min(a;i,X2):  Choose  smallest  of  Xi  or  X2,  v  ;ich  are  both  positive. 

If  Xx  —  X2,  then  min(xi,X2)  =  xx  =  X2. 


Also  recall  that  K  =  \/P  (1  -  p)  and  K'^  =  P{1  -  pf.  Then  we  have 


|s(i)|V(  =  r  s'{t)s'{t)’dt  =  \s(t)\'‘dt 


(217) 


which  has  the  following  six  integrals  given  by  Eq.  (215) 


q-1  q-1 


le  f\SB\Ht  =  P(i  -  pfbioT.'L  p’*”  [t- 

gz=0m=0  V 


max(^,  m) 


(218) 


f  \sisi\^dt 
Jo 


we  recall  the  result  of  Eq.  (117) 

r  stsiPhii^t  =  £  M*  +  £ 

(219) 

j,  -1  -1  (-’■d)-! 

Adt  =  p{i-pfY.  H  E  E  VP 

®  i—  —  Lo  5=— (l+t)g+l  t  =—Lo  m=:— (l+r)q  +  l 

in  which 

a+m,  1  fmm{g  +  {l  +  i)q,m  +  {l+r)q)\ 

H>0  =  P  Oo,.Oo,r  I  - -  I 

r  Bdt  +  rcdt  =  2K^  Bdt  = 

Jo  Jo  Jo 

-1  -1  -(j--l)5-l 

2P{l-pf  I]  £  X]  H  fP'^""bo,ibo,rGi{g,i,m,r) 

i=-Lo  p=“(l+t)g+l  t=-Lo  m--rq 


where 


Gi(^,i,m,r) 


-[{g  +  (1  +  i)q)  -  (m  +  rq)],  ioT  m  +  rq  <  g  +  {I  +  i)q 

=  I  ^ 

0.  otherwise 


-1  -(*-1)9-1  -1  -(r-l)9-l 


A-*  /  Vdt  =  P(1  -  ri*  ^  E  E  E  P’^bo.ib 

-'0  t=-Lo  g=-ig  r=-Lo  m=-rg 


max(5'  +  iq,m  +  rq) 


We  know  from  Eq.  (138)  that 


r  \sACi?dt  =  r  AA*dt  +  r  BB*dt  +  r  CC*dt 

Jo  Jo  Jo  Jo 

r  2Re[AC*]dt  +  r2Re[AB*]dt 

Jo  Jo 

T 

+K^  [  2Re[BC*]dt 

Jo 

So  we  have  the  following  six  ACI  components 

'j'  +A^/2  g— 1  -\-M/2  g— 1 

A'^/  AA-dt  =  P{l~pf  Y.  E  E  E/*”XV1 

,  ....  5=0  m=0 


A:=-M/2  ^  n=-M/2 


k^O 


n^O 


where 


1 


j{u>k  -  OJn) 


oj[('^k-'^n)T-ujk(g/f3)+i^„{m./P)] 


=  < 


_pjl{^k-<^n){l/0)m^(g,m)-Wk{g/0)+i^n{m/P)] 


ioT  k  ^  n 


^kflbkfi  ^  nnax(5',  m)j  ,  for  A:  =  n 

I  kfc.oP 


(220) 


(221) 


K- 


+M/2 

-1  +A//2 

-1  (-»-9)-l 

E 

E  E  E 

x;  E 

■--M/2 

'=-‘»=-('+'l«+'  ..-M/2 

T-=z^h  m=:~(l+r)g-l-l 

k^O 

n^O 

(222) 

where 


bk/b: 


n,r 


'gi[-^fc(5//?)+Wn(m//?)] 

j{0Jk-UJn) 


^eil(‘^lc-‘^")(l/^)™in(5+(l+^)'7,m+(l+r)g)] 


<f2 


for  A;  7^  n 


bk/blr 


oj[i<^k  /P){Tn-g)] 


min((?  +  (1  +  ^)q,  m  +  (1  +  r)q) 


for  k  =  n  J 
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X  (P3 


+M/2  -1  -(^-l)g-l  +M/2  -1  -(r-l)g-l 


f  cc-dt  =  p{i-pf  E  E  E  E  E  E 

9=-iq  __  m  =  -Tq 


k^-M/2 


=  ^M/2 
n^O 


where 


J[{^k  —'^n)T-Wk  {g/ P)+UJn(rn/ /?)] 


h,eK,r  - \  1^' 


_  J  _gj[(wfc-a-„)(l//3)max(ff+fg,m+r(j)-wt(ff//3)+w„(m//3)]  ^  £qP  k  ^  H 


h  tK  T  -  —  max{g  +  £q,m  +  rq)  ,  for  A:  = 


(223) 


+M/2  q-1  +M/2  -1  -(r-l)g-l 


f  2Re[AC-]  =  P(1  -  pf  E  E  E  E  E 

•^0  n  r=— r  771=— 1 


k^-M/2  n=-M/2 


(224) 


where 


J[{u;k-f^n)T-ujk(9/P)-^i^n(rn/p)] 


_gi[(t*^fc“^n)(l//?)niax(£r,7n+r9)-wfc(p//?)4-ct;n(m//?)] 


-  i  [(t*^fc  -  wn  )  T- Wfc  (^//3) +a-n  (m //?)] 


_g--i[(a>fc-Wn)(l/i0)max(5,m+rg)-Wfc(5//3)+u;n(m//?)]  ^  k^U 


T  -  -  ma.x{g,m  +  rq)  , 


for  k  =  n 
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+Af/2  q-l  +M/2  -1  (-r?)-! 


/  2R4AB-]  =  P(1  -  pf  E  E  E  E  E 

•^0  _ rk  _ r  , 


k=-M/2  n=-M/2  "=-£'m=-(l+r)9+l 


(225) 


where 


o6*  _  Lj[{^k-^n)({Tn-\-(l+r)q)/P)-ujk(g/P)+ujn{m/P)] 


-j[(<^k-^n){{m+{'l+T)q)/P)-Wk  {g/0)+i^n(m/P)] 


_e-j((‘^->//5)('"-s)l|  I  G2{g,m,r), 


for  A;  7^  n 


where 


-[(m  +  (1  +  r)q)  -g]j  G2{g,  m,  r), 


G2{g,m,r) 


for  <  m  +  (1  +  r)q 


0,  otherwise 


for  k  —  n 


J  +M/2  _1  (-^g)-l  +M/2  _1  -(r-l)?-l 

K^l  2Re[BC*]  =  P{l-py  E  E  E  E  E  E 

°  k=-M/2  ^=-^5=-(l+^)9+l  r.=-M/2  "  = 


(226) 
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where 


pjl(^k-^n){{Tn'\-rq)/p)-uJk{9/0)‘\-^n{m/p)] 


I  ^f7-l-mL*  L  ^  [^— if(wfc~a;n)((5+{l+^)g)//5)-a;fc(^//?)+cJn(77i//?)] 

+'’  '■‘■'‘"■'i-iK--,)) 

_e-il{'*'*~“'n)((m+rq)//J)-i^»:{9//J)+UJn{m//?)]  ^  t),  foF  k  ^  Tl 


-[((/  +  (1  +  £)q)  -  (m  +  r?)]  j  (?3((?,  £,  m,  r), 


where  G3(ff,£,m,r)  is  defined  in  Eq.  (172). 
Now,  for  the  last  three  terms  of 


for  k  = 


which  are 


r  \s{t)\^dt = r  s' {t)s'itydt = r 

Jo  Jo  Jo 

T  r'l' 

[  2Re[sBS*isi]dt-  I  2Re[s*BSACi]dt 

Jo  -^0 


r  2Re[s}siSACi]dt 
Jo 

T  /"  J' 

[  2Re[sBs}si]dt  =  /  2[sBSisi]dt  = 

Jo  Jo 


2P{l-pyj2  12  iZ  P^'^'^bo,obo,iG4{g,i,m) 

g=0  i=—Lo  m=— (14-1)9+1 


q-l  _i  -(2-1)9-!  r 

+2P{i-pfj:-£  E.  p^'^’”&o,o6o,i  h  “  ^  max(5f,  m  +  iq) 

g:=Oi=-Lo  m^—iq  L  ^ 


(227) 
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where 


Gi{g,i,m)  = 


1 


[(m  +  (l  +  0g)-5], 


0, 


g  <  m  +  {I  -]r  i)q 

> 

otherwise 


T  T 

f  2Re[sQSAci]dt  =  f  2Re[sBSACi]dt  =  Aba  +  Bba  +  Cba 
Jo  Jo 


q-1  +M/2  q-1 

Aba  =  ^(1  -  XI  X  *^7 

k=-M/2 

k^O 


in  which 


<f7 


’  ’  j<^k  ^  ^ 

I  „3+mA  /,*  ^  \p-j<^k(T-{m/g))  _  -jwk((ina.x{g,m)~m)/g) 


where 

‘fs 


where 


g-l  +M/2  -1  (-^?)-l 

Bba  =  p{i-p)^E  E  E  E  »>* 

5-0  7n=-(l+^)g+l 

A:  9^0 


P^^^bo,obk,i 


J^k 


JiAjt:((\+t)qlP)  _  pi<Aik((g-m)/p) 


G5{g,l,m) 


’  '—J<^k\  ^  ^ 


Gs{g,l,m)  =  I 


1, 

0, 


CBA  =  P{l~pf'Z 

p=0 


g  <  m  +  {\  -\-  t)q 

► 

otherwise  ^ 

+M/2  -1  -(^-l),-l 

X  X  X  <^9 

*#0 
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^  '  jl^k  ^  ^ 

_l_  fl+m^oo6]^^— -^—  e-i‘*’*(r-('"/'®))  _  (228) 

T  /•X' 

K^  [  2Re[s}siSACi]d't  =  f  "^.RelsjsiSAcAdt 

Jo  Jo 

=  AC  ft  +  ADfiAEfi  +  BCfi  +  BDr  +  BEpi 


-1  (-w)-l  +M/2  g-1 

/ic«  =  p(i-,>)^  E  E  E  E'P'O 

t=-io  ff=-(l+*)9+l  k=-M/2 


(229) 


where 


V^io  =  P' 


•^■^'"60 ,-6*0—  -  1  G6(^,i,m) 

’  ’  J^k  I 

+/>^+”^6o  or^-^  e-jw*[(l/^)(g+(l+.»-m//3]  _  j  (06(^,  I,  TTl) 

'  [-J^k\  [  X. 


where 


G&{g,i,m)  = 


1,  m<g  +  {l-\-i)q 

0,  otherwise 


Note:  e^j<^k(mlp-mli3)  ^  gjo  ^  ^ 


-1  (-tq)-l  +M/2  -1  (-M-1 

yi£)H  =  ^(i-/>)^  EE  E  E  E 

t=:-Lo  g=-(H-i)g+l  jt=-M/2  ^-~^™=-(l+^)9+l 


in  which 


g,,,  =  e- 


(230) 


i  .a+mL  L*  ^  f  iwfcf(l//3)  min(5+{l4-09i”^+(l+^)9)“’^/^]  _  ^//^] 

+'■  I'  1 
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(231) 


-1  (-tg)-l  +M/2  -1  -(£-1)9-1 

AE^  =  p(i-pf  Y.  E  E  E  E  v’n 

i=  —  Lo  ^=  — (l+O'J+l  k  =  —  M/2  m=~~£q 


in  which 


¥^12  —  —  e 


where 


G7{g,i,i,m) 


m  +  <  £/  +  (1  +  i)q 


0,  otherwise 


in  which 


-1  -(t-l)9-l  +M/2  9-1 

BCR  =  p(i-pf  E  E  E  Ev’is 

.=-io  9=-.9  -=0 

k^O 


(232) 


nfl+mt  u  ^  Ipi^^klT-imlP)]  _  iaiA:[(l//?)max(g+i9,m)-m/^] 

P  Oo  jOfc  o  •  C 

JWfc  L 


+/9^+’”6o,,^>fcO 


g-jai/5[T-(m//3)]  _  g-ja/fc[(l//3)  max(9+J9,m)-m//3] 


-1  -(«-l)9-l  +M/2  -1  (-£9)-1 

BZlR  =  l’(l-rt"  D  E  E  E  E  Vn  (233) 

t=-Lo  fc=-Af/2  ^=“^^=-(1+^)9+! 


in  which 


<^14  =  \l^^'"h,^ibk,i—  -  e 


^+"^6o.-^Lr  ^  1  (S'sC^,  m) 

’  ’  [-3<^k\  '■ 


lU 


where 


g  +  iq  <  m  {I  +  tjq 
otherwise 


where  %  is  the  responsitivity  of  the  photodetector  (A/W).  This  will  be  used  to  com¬ 
pute  probabilities  of  bit  error  for  the  dense  WDM  system. 
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APPENDIX  B 


LIMITED  CASE  OF  THE  COMPLETE  MODEL 

(10  PROGRAMMABLE  TERMS)  AND  RESULTING 
PROBABILITY  OF  BIT  ERROR  EQUATIONS 


In  Appendix  A  we  derived  the  complete  expression  for  \s{t)\‘^dt.  Referring  to  the 
OOK  receiver  structure  in  Fig.  1,  we  note  that  the  deterministic  signal  component 
of  the  random  variable  appearing  at  the  output  of  the  integrator  is 

X  =n  r  \s{t)\^dt  (235) 

Jo 

where  TZ  is  the  responsivity  of  the  photodetector  (A/W)  and  s{t)  is  the  complex 
baseband  output  of  the  Fabry-Perot  filter.  Note  that  the  photodetector  detects  s{t) 
and  produces  an  output  current  '}l\s{t)\^.  We  will  call  the  integrator  output  the 
decision  variable  Y  where 

Y  =  X  +  N  (236) 

We  note  that  N  is  &  Gaussian  random  variable  with  mean  zero  and  variance  NqT 

N  ~  N(0,  NoT)  (237a) 


and 

N=  r  n(t)dt  (237b) 

Jo 

where  n{t)  is  the  postdetection  additive  white  Gaussian  noise  and  No  (N^IRz)  is  the 
two-sided  current  spectral  density  of  n{t). 

Now  we  will  set  up  a  fairly  limited  but  realistic  case  for  which  we  will  program 
the  appropriate  terms  derived  in  Appendix  A  and  generate  the  probability  of  bit 
error.  Letting 


Lo  =  1  L  =  0 


=  0 


UJk 


2TrkI 

T 
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LOk  =  2'KkI /T  is  a  special  case  of  the  radian  frequency  spacing  between  Channel  0  and 
Channel  k.  Recall  again  that  /  is  the  normalized  channel  spacing  integer  (/  >  0),  and 
T  is  the  data  bit  period  (s).  L  =  0  and  <^k  =  0  means  that  we  will  only  model  and  deal 
with  adjacent  channel  interference  (ACI)  from  the  bits  in  the  adjacent  channels  which 
occur  during  the  detection  window  0  <  t  <  T.  Recall  that  the  adjacent  channels  are 
symmetric  in  frequency  around  Channel  0,  the  channel  of  interest.  Lq  =  I  means  that 
we  will  model  the  effects  of  a  single  ISI  bit  trailing  the  detected  bit  of  interest  feo.o  in 
Channel  0.  =  0  means  that  the  phase  offset  between  Channel  0  and  Channel  k  is 

0,  or  that  all  adjacent  channels  are  bit  synchronous  with  with  Channel  0. 

Then,  since  (pk  =  0,  L  =  0,  and  bk/  £  {0,6^'^'=},  we  let 

^-A//2,o  =  ^-M/2+1,0  =  •  •  •  =  6-10  =  b  and  b  G  {0, 1}  (238) 

b~  is  the  left  adjacent  channel  bit  pattern  and  all  M/2  of  the  left  channel  bits  will  be 
1  or  0  simultaneously  during  detection  window  0  <  t  <T.  And  again,  since  <j)k  =  0, 
L  =  0,  and  bk,e  €  {0,  e'^*},  we  let 

bM/2fl  -  bM/2~i,o  =  ■■■  =  bi^o  =  b'^'  and  b'^  G  {0, 1}  (239) 

is  the  right  adjacent  channel  bit  pattern  and  all  M/2  of  the  right  channel  bits  will 
be  1  or  0  simultaneously  during  the  detection  window  0  <  t  <  T.  Also,  we  know 
Lo  =  1,  and  6o,i  G  {0, 1},  so  using  Eq.  (1)  or  Eq.  (41)  we  have 

6o,-i  e  {0,1}  and  6o.o  e  {0,1}  (240) 

We  define  an  ACI/ISI  bit  pattern 

=  {6-,6+,6o,-i}  (241) 

Note  that  b~,  and  6o_i  are  0  or  1  with  probability  1/2  yielding  eight  possible 
values  of  ipp. 
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Now  let  use  define  an  ACI/ISI  bit  pattern  set  if) 


V’ =  (242) 

We  can  denote  each  individual  element  in  the  set  tp  as  x/jp  where  p  =  1, 2,  3, . . . ,  8  as 
there  are  eight  possible  values  of  i/^p,  and 


0  =  {V’p} 


(243) 


where  p  =  1,  •  •  • ,  8.  Now  we  define 


Xo(V’p)  =  X{xf}p,bQfl  =  0)  or  X  evaluated  at  the  current  value  of  V’p  with  6o,o  =  0 

(244) 

Xi{ipp)  =  X{xJjp,bofl  =  1)  or  W  evaluated  at  the  current  value  of  t/)p  with  6o,o  =  1 

(245) 


It  can  be  shown  that  the  conditional  probability  of  error  for  the  dense  WDM 


system  given  the  ACI/ISI  bit  pattern  xpp  is  [2,  4] 

P(error|^,)  =  -Q  I  ysjp  “  )  +  2^  . 


(246) 


where 


There  are  eight  possible  members  of  the  set  xp 


xp  = 


Vt  is  the  detection  threshold  where 


II 

o 

0 

0 

05  =  1 

0 

0  ' 

^2=0 

0 

1 

06  =  1 

0 

1 

CO 

II 

o 

1 

0 

07  =  1 

1 

0 

II 

O 

1 

1 

08  =  1 

1 

1 . 

(247) 


(248) 


Vt  = 


Xo  max  +^1  min 


(249) 
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and 

Xq  max  =  max  {Xo{ipp))  (250) 

i'Pp} 

i.e.,  to  find  Xq  max,  compute  all  eight  values  of  Xq  [see  Eq.  (244)]  and  then  choose 
the  maximum,  and 

Ximin  =  mn(Xi(?/>p))  (251) 

i.e.,  to  find  Xi  min,  compute  all  eight  values  of  Xi  [see  Eq.  (245)]  and  then  choose 
the  minimum.  Then,  by  the  law  of  total  probability,  the  probability  of  bit  error  is 

Pb  =  P  (errorjz/^i)  P(V’i)  +  P  (errorlV’2)  -^(^2)  H - h  P  (errorj^g)  PH’s)  (252) 

Assuming  all  bit  patterns  are  equiprobable 

P(^,)  =  P(M  =  ■■■  =  P{4's)  =  i  (253) 


Then  we  can  easily  see  that 


Again,  we  note  that 


Pft  =  -  ^  P  (errorIV'p) 


x  =  n  \s(t)\‘^dt 
Jo 


(254) 


(255) 


where  \s{t)\^dt  was  computed  in  Appendix  A  and  consisted  of  21  terms  (“clusters 
of  summations”)  [see  Eqs.  (218)-(234)].  Given  our  parameters  Lq  =  I,  L  =  0,  (l)k  =  0, 
and  Uk  =  2'KkIlT,  11  of  the  21  terms  drop  out  and  we  are  left  to  compute  in  order 
ten  to  generate  the  probability  of  bit  error.  If  we  name  these  ten  terms  (“clusters  of 
summations”)  SUMj,  SUM2,  •  •  • ,  and  SUMio,  then  for  Lq  =  1,  L  =  0,  <l)k  =  0,  and 


LOk  =  2‘KkI  fT 


r  \s(t)\^dt  =  SUMi  +  SUM2  +  •  •  •  +  SUMio 
Jo 


(256) 


and  using  Eq.  (255) 


a:  =  7^(SUMl  +  SUM2  +  •  ■  ■  +  SUMio) 


(257) 


no 


Now  PT  may  be  factored  out  of  each  sum.  So  we  define 


SUM,- =  Pr(QUOTj)  i  =  (258) 

So 

X  =  7^P^(QU0Tl)  +  7^P^(QUOT2)  +  •  •  •  4-  7^P^(QUOTlo)  (259) 

and 

X  =  n  Pr[QUOT,  +  QUOT2  +  •  •  •  +  QUOTio]  (260) 

Also,  letting  QUOT,  +  QUOTj  +  •  •  •  +  QUOT^q  =  QUOTSUM,  we  also  see 

X  =  n  Pr[QUOTSUM]  (261) 

We  can  also  say  that,  in  general,  QUOT,,  QUOT2,  •  •  • ,  QUOT^q,  and  QUOTSUM  are 
dependent  functions  of  ?/>p  and  60, o-  Thus,  we  may  write  in  general 

QUOTi(^p,  Vo),  QUOT2(V’p,  60,0),  •  •  • ,  QUOT,o(0p,  60,0) 


and 

QUOTSUM(  V,  60,0) 

Now  we  can  see  that 


A:i(V>p)  =  7^P^[QUOTSUM(V,Vo  =  1)]  (262) 

Xo(V)  =  7^P^[QUOTSUM(V,Vo  =  0)]  (263) 


Recalling  that 


Vt 


Xq  max  +A'i  min 
2 


we  let 


V’p  max  =  Value  of  ipp  which  causes  maximum  value  of 

of  Xo  (or  maximum  value  ofX  with  60,0  =  0).  (264) 
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and 


V’p  min  =  Value  of  ipp  which  causes  minimum  value  of 

of  Xi  (or  minimum  value  ofX  with  6o,o  =  !)•  (265) 


Then  using  Eqs.  (249),  (262)-(265) 


Vx  = 


7^PT[QUOTSUM(0p  ^ax,6o.o  =  0)  +  QUOTSUM(i/>p  ^,6o.o  =  1)] 

2 


(266) 


Now  using  Eqs.  (246),  (262),  (263),  and  (266) 
P{eTrov/ipp)  = 


\q  (lQUOTSUM(.f,.6„,„  =  1)1 

QUOTSUM  {<p,^,bo„  =  0)  +  QUOTSUM  (V-,  6o,o  =  1) 


+ 


nPT\ 

vwr) 


QUOTSUM^  „,ax,  bofl  =  0)  +  QUOTSUM  (V’p  min,  &0,0  =  1) 


-[QUOTSUM(V>p,6o.o  = 


(267) 


Thus,  we  see  that  P(error|0p)  is  directly  related  to  the  signal  to  noise  ratio  which  we 
will  call  Z. 

^  nPT  ^  r. 

~'7w~ 

Thus,  when  we  compute  a  probability  of  the  bit  error  graph,  we  will  choose  a  suitable 
range  of  values  for  Z,  a  value  of  free  spectral  range-bit  period  product  and  several 
values  of  the  number  of  adjacent  channels  M.  For  each  of  these  values  of  M,  we  will 
compute  a  probability  of  bit  error  trace.  To  compute  a  point  on  a  trace,  choose  a  value 
of  Z,  compute  all  eight  values  of  P (error [i/’p),  sum  all  eight  values  of  P(error/^p),  and 
divide  by  8  [see  Eqs.  (248),  (254),  and  (267)].  The  point  is  then  plotted. 
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Recall  that 


then 


q  =  /3T 

We  will  use  this  and  the  fact  that  [see  Eq.  (259)] 


(269a) 


(269b) 


x  =  n  pr(QU0Ti ) + n  PT{q\]OT^)  +  ---  +  n  pr(QUOTio)  (270) 

V _ ^  ^  ^ 

Term  #1  Term  #2  Term  #10 

to  compute  the  probability  of  bit  error  graphs.  We  also  use  Eqs.  (256),  (257),  (258), 
(260),  (261),  (262),  (263),  (264),  (265),  (266),  and  (267)  in  this  endeavor. 

We  now  proceed  to  write  out  the  ten  terms  of  X,  which  we  will  use  to  compute 
the  probability  of  bit  error  graphs.  Before  we  begin  we  recall  that  max  {xi,X2) 
means  choose  the  largest  of  Xi  or  X2.  If  X2  =  X21  then  max  (xi,X2)  =  xi  =  X2- 
Also  recall  that  min(xi,X2)  means  choose  the  smallest  of  Xj  or  X2.  If  Xi  =  X2,  then 
min  (xi,  X2)  =  xi  =  X2.  The  first  term  to  be  programmed  will  be  Eq.  (218),  Appendix 
A.  Factoring  out  a  T  we  get 

SUM.  =  Fr(l  -  (1  -  (""-jf.’--))  (2U) 

Using  q  =  fST  and  multiplying  by  TZ  we  get 

Tern.  #1  =  UPTH  -  (l  -  (272) 

With  Lo  =  I  and  q  =  ^T,  the  first  term  (“cluster  of  summations”)  in  Eq.  (219), 
Appendix  A,  becomes 

j3T-l  pT-1  / 1  \ 

^(1  -pf  Y.  Y  iK-if  (  -3  ™n(5r,  m) )  (273) 

g=l  m=l  \P  / 

Factoring  out  a  T  and  multiplying  by  TZ  we  get 

iST-lpT-l  j 

Term  #2  =  71PT(1  -  pY{bo,-if  Y  Y  /■^’"^min(^,m)  (274) 

3=1  m=l 
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With  Lq  =  I  and  q  =  I3T^  the  second  term  (“cluster  of  summations”)  in  Eq. 
(219),  Appendix  A,  becomes 

PT-l 20T-1 

2P{l-pfYl  E  /■'”(V-.)“G,(s,i  =  -l,m,r  =  -l)  (275) 

5=1  m=(3T 


where 


= -1)  =  I  ^ 


[g-{m-  /?r)], 


for  m  —  <  g 


otherwise 


Factoring  a  T  and  multiplying  by  7^,  we  arrive  at  the  expression  for  Term  #3 

/5T-1  2/5T-1  P , 

Term#3  =  RPTx2(l-rf(V-,)'  E  E  ~  X  G.(s, !  = -1, m, r  = -1) 

5=1  m=pT 

(276) 

With  To  =  1  and  q  =  /ST,  the  third  term  (“cluster  of  summations”)  in  Eq. 
(219),  Appendix  A,  yields  Term  ^4,  after  factoring  a  T  and  multiplying  by  TZ. 

2/3T-1 2/3T-1  r 

Term#4  =  77Fr(l-p)=(6„,_0"  E  E  P’*"  i  -  -  IST,m  -  PT) 

g=0T  m=gT  L 

(277) 


The  next  applicable  term  is  Eq.  (221)  of  Appendix  A.  Our  given  parameters 
include  4>k  =  (j>n  —  0.  We  can  then  say  by  recalling  that  hk/  G  {0, 


^n,o  —  ^71,0  and  6^  o  —  bk^i 


(278) 


We  will  also  use  the  facts 


2TrkI 

LOk  =  — and  = 


27rn/ 

T 


(279) 


where  /  is  a  positive  integer  called  the  normalized  channel  spacing  integer.  Now,  by 
slightly  reworking  the  ujk  ^  c<;„  expression  by  substituting  the  results  given  in  Eqs. 

(278)  and  (279)  into  the  appropriate  places,  we  get 


-  hflbnfi  - -TT  e 


;[27r/(A:  -  n)] 


j[27r/(A:-n)-(27rA)/(g)//3r)  +  (27r7i/(m)//3T)] 


J[{2'kI / (iT){k-n)  max(5,m)~(27rA;/(5)//3T)  +  (27rn/(m)//?T)] 


for  A;  n  (280) 
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Note:  The  2-Kl{k  —  n)  term  in  the  first  phasor  of  Eq.  (280)  contributes  integer 
multiples  of  27r  to  the  phasor,  and  so  it  is  neglected  in  the  next  rewrite  of  the  term. 
Then,  by  factoring  the  common  factor  we  arrive  at 


(^1  = 


1  _  ej(2^///3T)(fc-n)tnax(s,m)  ^  for  A:  ^  U  (281) 


j[27r/(A:  -  n)] 

Finally,  after  we  factor  out  a  T  from  both  the  ^  u;„  and  Uk  =  u)n  arguments  to 
the  summations,  rearrange  the  summations  for  more  logical  programming,  and  use 
the  fact  that  q  =  I3T,  we  arrive  at  the  final  form  for  Term  #5. 

+M/2  +M/2  PT-1I3T-1 

Term  #5  =  7^P^(1  -  pf  Y.  E  iKo){Kfl)  E  E  ^  (282) 

k=-M/2  n=-M/2 

k^Q  n^O 


where 


=  < 


J(2irI/PT)[{nm)-{kg)] 

j[2'Kl{k  -  n)] 


,j[{2KkI/PT)(m-g)] 


1  - 


^  _  ^j(2irI/0T)(k-n)max(g,m) 

ma,x{g,m) 


I3T 


iox  k  ^  n 


for  k  =  n 


The  next  applicable  term  is  the  first  term  (“cluster  of  summations”)  in  Eq. 
(227),  Appendix  A.  Utilizing  Lq  =  1,  factoring  a  T,  multiplying  by  TZ,  and  noting 
that  q  =  I3T  yields 

0T-1PT-1  _  _i 

Term  #6  =  7^P^  x  2(1  -  p)\bo,o){bo,-i)  E  E  P'^"  '  r - ^ 

g=Q  m=l 


where 


04(9, i  =  -l,m) 
T 


=  < 


1 

Jr 


\m-g], 


iox  g  <ra 


0,  otherwise 

Looking  at  the  second  term  (“cluster  of  summations”)  in  Eq.  (227),  Appendix 

A,  and  after  utilizing  To  =  1,  factoring  a  T,  multiplying  by  and  substituting 
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q  =  /?r,  we  have 


/3T-1  2f5T-\ 


Term#7  =  Ki>Tx2(l-rt"(Vo)(V-i)  E  E  1^*"  1  -  ^  maxfc m  - /3T) 

5=0  m=0T  L 


(284) 


The  next  applicable  term  is  the  first  term  (“cluster  of  summations”)  in  Eq. 
(228),  Appendix  A.  Since  (j)k  =  0 


^k,0  ~  ^k,o 


(285) 


We  can  then  work  with  the  term  inside  of  the  summations 


(pj  —  P^'^'^KobkO _  \e^‘^k(T-(m/0))  _  ^jUk({m&x(g,m)-m)/0) 

’  ’  jtOk  ^ 


+p^+"^6o,0^fc,0^  :  -y  -  e->*((max(5.m)-m)//3)j  ^286) 


Then,  since  big  —  bkfi  and  after  some  rearranging,  we  have 


P^^^bo,obk,o 

P7  =  - ^ ^  X 

Uk 


J‘^k{T-{Tn/0))  _  Q-jujk(,T-{m/0))  _^juJk{(max{g,m)-m)/0)  g-jwfc((max(5,m)-m)//g) 


Then  factoring  a  —1  from  A  yields 


(287) 


Jwk{imax(g,m)—m)/0)  _  -juJk{(max{g,Tn)-Tn)/0) 


(288) 


Euler’s  Relationship  defines 


(289) 


Applying  this  relationship  yields 


2p^+’"6o,o^>fc,o  [  .  [  im  ■  f  /max(5^,m)-m 

V>r  = - - U.„  1^0,  Jr  -  ^  1  1  -  sm  U  I - ^ - 


(290) 


and  since  Uk  =  2'KkI jT 


After  factoring  a  T,  multiplying  by  7?.,  utilizing  q  =  and  rearranging  the  summa¬ 
tions  for  more  logical  programming,  we  arrive  at  the  expression  for  Term  ^8. 

-l-M/2  fiT-\  0T-\ 

Term  #8  =  'RPT{1  -  pf  S  Z)  *^7 

k=-Ml2 

A:^0 

where 

V?7  =  ~  w)  ~  (292) 

The  next  applicable  term  (“cluster  of  summations”)  is  Eq.  (229)  of  Appendix 
A.  Since  4)k  =  0,  we  know  6fc,o  =  K,o-  Using  Lq  =  1,  we  can  rework  the  argument  of 
the  summations.  Making  the  appropriate  parameter  substitutions  yields 

JWfc  L  J 

[e-i-d(.-m)//3]  _  i]  i  =  -1,  m)  (293) 

Further  symplifying  yields 


ipiQ= - Ge{g^i  — —l,m)  -. 

(^k  J 


2sm[ujk[(9-rn)/l3]] 


Using  ujk  =  2'jrkl fT  yields 


Tp^+”^{bo,-i){bk,o)^  .  .  IX.  \2TrkI 

q^w  =  — - Vib/  ’  = -l,"7)sin 


(294) 


(295) 


Substituting  (pio  inside  of  the  summations,  multiplying  by  TZ,  factoring  T,  using 
q  =  I3T,  and  rearranging  the  summations  for  logical  programming  yields  the  following 
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for  Term  #9. 


+  M/2  pT-iPT-1 

Term#9  =  7^P^(l-p)^  ^  E  E  ¥>io 

k=-M/2  5=^  ”^=° 
yfc^O 


(296) 


where 


amd 


/+-(6o,-i)(Vo)^  ,  .  , 

(^10  = - - - Geig^t  =  -l,m)sm 


TT^/ 


2'KkI 

TF 


{g  -  m) 


1, 

0, 


for  m  < 
otherwise 


G6{g,i  =  -l,m)  = 

Look  at  the  last  applicable  term  (“cluster  of  summations”),  Eq.  (232),  Appendix 
A.  Again  we  take  advantage  of  the  fact  that  <j)k  =  0,  which  means  o  =  Kfi-  Making 
the  appropriate  parameter  substitutions  yields 


"I  i  ^ 


{m/P)]  _  ^-ji^k[T-{Tn/p)] 


JbPk[(.'i-/0)max{g-PT,Tn)-m/P]  _  -juik[{l/P)ma\(g-pT,m)—m/P] 


11 


(297) 


Using  Euler’s  Relationship  we  get 

rp^+-(6o,-i)(Vo) 


<^13 


TrkI 


f  m\ 

2'KkI  1 

sin 

—  sin 

(max((/  pT,m)  m) 

(298) 


Substituting  <fi3  back  inside  of  the  summations,  factoring  a  T,  multiplying  by  TZ, 
utilizing  q  —  (5T,  and  rearranging  the  summations  for  logical  programming  yields 


+A//2  2PT-IPT-1 

Term#10=RPr(l-rt^  Y,  T,  T.  Viz 

i=-Ml2  "=“ 

fc^O 


(299) 


where 


V^i3  — 


irkl 


/  m\ 

2KkI ,  , 

sin 

1 - 

1 

^1  ■ 
1 _ 

—  sin 

{ma,x{g  /3T,m)  m) 
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For  the  given  parameters  Lq  =  I,  L  =  0,  (l)k  =  0,  and  Uk  =  we  have 

now  derived  all  ten  programmable  terms  of  X  where 


Term  #1  Term  #2  Term  #10 


(300) 


We  will  use  these  ten  programmable  terms  to  compute  probability  of  bit  error  for  the 


dense  WDM  system. 
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APPENDIX  C 


PROGRAMMING  STRATEGY  AND  COMPUTER 
PROGRAMS  FOR  GENERATION  OF  PROBABILITY 
OF  BIT  ERROR  GRAPHS 


In  this  appendix,  we  present  the  computer  programs  written  to  generate  the  proba¬ 
bility  of  bit  error  graphs  for  the  system  under  consideration.  We  attempt  to  generate 
four  graphs,  one  for  each  value  of  considered.  We  reiterate  that  /3  is  the  free 
spectral  range  of  the  Fabry-Perot  filter  (Hz)  and  that  T  is  the  data  bit  period  (s). 
/?r  is  called  the  free  spectral  range-bit  period  product.  The  four  values  of  jiT  used 
are 

I3T  =  [500  1000  1500  2000] 


In  each  of  the  four  graphs  there  will  be  five  traces  as  we  will  present  probability  bit 
error  for  four  values  of  the  normalized  channel  spacing  integer  /,  or  equivalently  the 
number  of  adjacent  channels  M,  along  with  a  probability  of  bit  error  trace  for  single 
channel  (SC)  operation  without  Fabry-Perot  (FP)  filtering  or  SC  operation  with  FP 
filtering  and  without  ISI  or  ACL  We  will  show  the  relationship  between  I  and  the 
number  of  adjacent  channels  M  later.  Now,  however,  we  present  the  four  values  of  / 
corresponding  to  each  value  of  /3T ; 


For  I3T  =  500 


7  =  [4  5  8  20] 


For  /3T  =  1000 


7  =  [5  6  9  20] 


For  I3T  =  1500 


7  =  [7  9  12  20] 
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For  I3T  =  2000 


/  =  [8  9  12  20] 


Without  a  lengthy  discourse,  we  present  the  mathematical  relationship  between 
the  normalized  channel  spacing  integer  /  and  the  number  of  adjacent  channels  M 

M=^-l  (301) 

For  the  special  case  fk  =  kAf  =  kl/T  we  can  easily  see  that  A/  =  I /T.  Then 

BT 

M  =  ^-l  (302) 


For  the  values  of  /  above  M  will  not  always  be  an  even  integer.  Thus,  an  algorithm 
to  consistently  arrive  at  an  even  integer  value  of  M  which  is  less  than  or  equal  to  the 
true  mathematical  value  presented  in  Eq.  (302)  was  devised.  To  get  the  number  of 
adjacent  channels,  we  perform  the  following  operation 


BT 


—  Q  R 


(303) 


where  Q  is  the  integer  quotient  of  the  division  operation  and  R  is  the  remainder.  To 
arrive  at  M  we 


1.  Subtract  1  from  Q  if  <5  is  an  odd  integer  and  R  =  0. 

2.  Subtract  2  from  Q  if  <5  is  an  even  integer  and  R  —  0. 

3.  Subtract  1  from  Q  if  Q  is  an  odd  integer  and  R  ^  0. 

4.  Subtract  2  from  if  <5  is  an  even  integer  and  R  ^  0. 

Using  these  rules,  we  obtain  at  the  four  values  of  M  for  each  value  of  BT. 
For  BT  =  500 

M  =  [m  ^  ^ 

1=4.  1=5  I=S  1=20 
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For  (iT  =  1000 

For  I3T  =  1500 

For  /3T  =  2000 


M  =  [m  1^  m 

7=5  7=6  7=9  7=20 


M  =  [^  ^  m  7^] 

7=7  7=9  7=12  7=20 


M=[^  ^  ^ 

7=8  7=9  7=12  7=20 

Now  we  use  these  four  values  of  I  and  four  values  of  M  for  each  value  of  l3T 
to  compute  the  four  graphs  of  probability  of  bit  error  with  five  traces  each.  We  will 
use  the  equations  and  methods  developed  in  Appendix  B  to  compute  these  graphs. 
Recall  that 


X  =  'R  Pr(QUOTi)  +  n  Pr(QUOT2)  +  •  •  •  +  7^  PT(QUOTio)  (304) 

' - ^  ^ '  ' - V - ' 

Term  #1  Teim  #2  Term  #10 

Realizing  this,  we  program  the  ten  terms  given  in  Appendix  B  in  the  following  way: 

TERM  #1  [Eq.  (272)]:  We  will  compute  the  QUOT^  portion  of  this  term  four  times, 
once  for  each  value  of  /3T  =  500,  1000,  1500,  and  2000.  We  will  use  each 
of  these  four  values  in  a  separate  program  to  compute  the  probability  of 
bit  error  according  to  the  equations  developed  in  Appendix  B. 

TERM  #2  [Eq.  (274)]:  We  will  compute  the  QUOTj  portion  of  TERM  #2  four  times, 
once  for  each  value  of  I3T  =  500,  1000,  1500,  and  2000.  Each  of  these 
values  will  be  used  in  a  separate  program  to  compute  the  probability  of 
bit  error  according  to  the  equations  developed  in  Appendix  B. 

TERM  #3  [Eq.  (276)]:  We  will  compute  the  QUOT3  portion  of  TERM  #3  four 
times,  once  each  value  of  =  500,  1000,  1500,  and  2000.  Each  of  these 
values  will  be  used  in  a  separate  program  to  compute  the  probability  of 
bit  error  according  to  the  equations  developed  in  Appendix  B. 
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TERM  #4  [Eq.  (277)]:  We  will  compute  the  QUOT4  portion  of  TERM  #4  four  times, 
once  for  each  value  of  —  500,  1000,  1500,  and  2000.  Each  of  these 
values  will  be  used  in  a  separate  program  to  compute  the  probability  of 
bit  error  according  to  the  equations  developed  in  Appendix  B. 

TERM  [Eq.  (282)]:  We  will  compute  the  value  of  QUOT5  a  total  of  48  times,  12 
times  for  each  value  of  =  500,  1000,  1500,  and  2000.  For  each  value 
of  /3r,  we  have  the  four  corresponding  values  of  the  normalized  channel 
spacing  integer  /,  or  equivalently,  the  number  of  adjacent  channels  M. 
We  will  compute  the  value  of  QUOT5  three  times  for  each  value  of  I. 
QUOT5  will  be  computed  once  for  all  M/2  of  the  lower  adjacent  channels 
being  packed  with  Is  and  the  upper  adjacent  channels  being  packed  with 
Os.  This  is  the  case:  (6~  —  l,^"''  =  0).  QUOT5  will  be  computed  once  for 
the  lower  adjacent  channels  being  packed  with  Os  and  the  upper  adjacent 
channels  being  packed  with  Is.  This  is  the  case:  (6“  =  0,  6"^  =  1).  Finally, 
we  will  compute  QUOT5  once  more  for  both  the  lower  and  upper  adjacent 
channels  being  packed  with  I’s.  This  is  the  case:  [b~  =  =  1).  These 

twelve  values  of  QUOTg  for  each  value  of  will  be  used  in  each  of  the 
four  separate  bit  error  programs  to  compute  each  of  the  four  multiple 
channel  probability  of  bit  error  traces  using  the  equations  and  methods 
developed  in  Appendix  B. 

TERM  ^6  [Eq.  (283)]:  We  will  compute  the  QUOTg  portion  of  TERM  ^6  four  times, 
once  each  for  =  500,  1000,  1500,  and  2000.  Each  of  the  values  will  be 
used  in  separate  to  compute  program  probability  of  bit  error  according 
to  the  equations  developed  in  Appendix  B. 
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TERM  41^1  [Eq.  (284)]:  We  will  compute  the  QUOT7  portion  of  TERM  #7  four  times, 
once  for  each  value  of  13T  —  500,  1000,  1500,  and  2000.  Each  of  these 
values  will  be  used  in  a  separate  program  to  compute  the  probability  of 
bit  error  according  to  the  equations  developed  in  Appendix  B. 

TERM  #8  [Eq.  (292)]:  We  compute  the  value  of  QUOTg  twelve  times  for  each 
value  of  I3T.  These  twelve  values  of  QUOTg  for  each  value  of  I3T  will  be 
used  in  each  of  the  four  separate  bit  error  graphing  programs  to  compute 
each  of  the  four  multiple  channel  probability  of  bit  error  traces  using  the 
equations  and  methods  developed  in  Appendix  B. 

TERM  #9  [Eq.  (296)]:  We  compute  the  value  of  QUOT9  twelve  times  for  each 
value  of  jST.  These  twelve  values  of  QUOT9  for  each  value  of  will  be 
used  in  each  of  the  four  separate  bit  error  graphing  programs  to  compute 
each  of  the  four  multiple  channel  probability  of  bit  error  traces  using  the 
equations  and  methods  developed  in  Appendix  B. 

TERM  #10  [Eq.  (299)]:  We  compute  the  value  of  QUOTjo  twelve  times  for  each 
value  of  I3T.  These  twelve  values  of  QUOTjo  for  each  value  of  /3T  will  be 
used  in  each  of  the  four  separate  bit  error  graphing  programs  to  compute 
each  of  the  four  multiple  channel  probability  of  bit  error  traces  using  the 
equations  and  methods  developed  in  Appendix  B. 

Note:  When  we  say  we  are  computing  QUOT^,  we  are  not  being  exactly  mathemat¬ 
ically  correct,  as  we  compute  each  of  these  terms  without  the  bit  values  5o,o  and/or 
bo -I  the  factor  (1  —  p)^.  These  are  accounted  for  in  the  final  programs.  For  example 

Term  #1  =  KPT{1  -  pfC  eJ  E/*"  (l  -  (‘""jr"')) 

' - -  '  . . ^ 

QUOTi 
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However,  the  computer  program  for  Term  #1  (QUOTj)  only  computes 


(306) 


four  times  for  (ST  =  500,  1000,  1500,  2000.  (1  -  pY  and  the  value  of  (5o,o)^  are 
accounted  for  in  the  final  programs  which  utilize  the  equations  and  methods  developed 
in  Appendix  B.  The  author  apologizes  for  the  slight  stretch  of  the  truth,  but  it  seemed 
necessary  to  succinctly  explain  the  general  method  of  computing  each  of  the  ten  terms, 
and  the  final  graphs. 

For  completness  we  also  present  the  probability  of  bit  error  equation  for  a  single 
channel  operation 

Pb  =  Q  (307) 

where  we  recall  [Eq.  (271),  Appendix  B]  that 


Z 


HPT 

VWT 


=  np 


(308) 


and 

Q{x)  =  ^l^  (309) 

The  computer  programs  for  each  of  the  ten  terms,  the  numerical  results,  and  the  final 
graph  are  now  presented. 
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%  THESIS  COMPUTER  WORK 

•/.  _ 

•/.  I  I 

%  I  TERM  #1  I 

y.  I _ I 

•/. 

y,  COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES  OF 
TERM  #1  AND  TO  WRITE  THESE  VALUES  TO  DIARY  FILE  FOR  LATER 
USE/MANIPULATION  IN  CALCULATING  THE  DETECTION  STATISTIC 
AND  THE  PROBABILITIES  OF  BIT  ERROR 


y. 
y. 
y. 
y. 
y. 

y.  JOHN  A.  STUDER 

y.  CPT,  U.S  ARMY 

y,  550-53-7181 


DATE  LAST  MODIFIED:  11  SEP  94 


rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

termlsum(i)=  0; 
for  g  =  l:betatau(i) 
for  m  =  l:betatau(i) 

termlsum(i)  =  termlsum(i)  +  (rh.o . '  ((g-l)  +  (m-l)) * .  .  . 

(l-(max( [g-1  m-1] ' )/betatau(i))) ) ; 

end 

end 

end 

diary  johnmanl.txt 

diary  on 

termlsum 

diary  off 

end 
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termlsum  = 


l.Oe+03  * 


7.0510 

8.5126 

9.0083 

9.2563 

y,  numbers/ arrows/words 
y,  below  added  after  matlab 

y  dumped  answers  to  file  using 
y  the  diary  commaind 

500 

1000 

1500 

2000 

-  betatau  values 

y,  The  numbers  displayed  above  are  the  computed  values  of 
y,  Term  #1  for  the  four  values  of  betatau  given  above. 

y,  These  values  will  be  used  later  in  other  progreims  to  compute  detection 
y,  statistics  and  the  probabilities  of  bit  error  for  various  signal  to  noise 
y,  ratios. 


I  TERM  #1  I 
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•/. 

•/. 


THESIS  COMPUTER  WORK 


y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


I  TERM  #2  I 


COMMENTS;  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES 
OF  TERM  #2  CORRESPONDING  TO  THE  FOUR  VALUES  OF 
betatau;  500,1000,1500,2000.  WE  WILL  USE  BACKGROUND 
PROCESSING  TO  WRITE  THESE  FOUR  VALUES  TO  AN  OUTPUT 
FILE  CALLED  teriii2.out.  WE  WILL  USE  THESE  FOUR  OUTPUT 
VALUES  FOR  TERM  #2  TO  LATER,  IN  ANOTHER  PROGRAM 
COMPUTE  THE  DETECTION  STATISTICS  AND  PROBABILITIES 
OF  BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER  DATE  LAST  MODIFIED:  11  SEP  94 

CPT,  U.S.  ARMY 

550-53-7181 


rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

term2sum(i)  =  0; 
for  g  =  1 :betatau(i)-l 

for  m  =  1:  betatau(i)-l 

term2suni(i)  =  term2sum(i)  +  ( (rho “  (g+m) ) *10111(5,111) )  ; 

end 


end 

terin2sumf  inal(i)  =  term2sum(i)/betatau(i)  ; 

end 

t  erm2  sumf inal 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commands  to  get  started;  intro,  demo,  help  help 

Commands  for  more  information;  help,  whatsnew,  info,  subscribe 


term2sumf inal  =  The  numbers/cirrows/words 

%  below  were  added  after 

959.5662  492.4271  328.3413  246.2563  %  the  matlab  background  job 


'/,  dumped  the  values  to  the 

~  ~  “  ’/,  file  term2.out 

'/,  500  1000  1500  2000  */,  All  values/text  not  done  by 

*/,  ^  ’/,  Sun  Stn  is  preceded  by  a  "7," 

7,  - betatau  values - 


29965020  flops. 


THESIS  COMPUTER  WORK 


y.y////j//////////.7////.y.y.7//.y.7.y;/.y.y. 


y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y, 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


1  TERM  #3  I 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES  OF 
TERM  #3  CORRESPONDING  TO  THE  FOUR  VALUES  OF 
betatau:  500,1000,1500,2000.  WE  WILL  USE  BACKGROUND 
PROCESSING  TO  COMPUTE  THE  VALUES  AND  WRITE  THESE  FOUR 
VALUES  TO  A  FILE  CALLED  termS.out.  WE  WILL  USE  THESE 
FOUR  VALUES  FOR  TERM  #3  TO  LATER,  IN  ANOTHER  PROGRAM 
COMPUTE  THE  DETECTION  STATISTICS  AND  PROBABILITIES  OF 
BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  12  SEP  1994 


format  short  e 
rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

term3sum(i)  =  0; 
for  g  =  1 :betatau(i)-l 

for  m  =  betatau(i) : ((2*betatau(i))-l) 
if  m-betatau(i)  <  g 

term3sum(i)  =  term3sum(i)  +  ((rho“(g+m))*(g-(m-betatau(i)))) ; 

end 


end 

end 

termSsumf inal(i)  =  (2*term3sum(i) )/betatau(i) ; 
end 

t  erm3  sumf inal 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 


term3sumf inal  =  */,  Numbers/arrows/words 

*/i  below  were  added  after 

1.2211e+01  4.2917e-02  1.8815e-04  9.2720e-07  */,  the  matlab  background 

X  job  dumped  the  values 

*/.  '  ^  ^  ~  */,  to  the  file  term3.out 

*/•  500  1000  1500  2000  */,  All  values/text  not 

’/•  ’/,  done  by  Sun  Stn.  is 

y.  y,  preceded  by  a 

y,  - betatau  values - 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y.  y. 

y.  TERM  #3  y. 

y.  y. 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


29990004  flops. 
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THESIS  COMPUTER  WORK 


•/. 

7. 


7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


1  TERM  #4 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES 
OF  TERM  #4  CORRESPONDING  TO  THE  FOUR  VALUES  OF 
betatau:  500,1000,1500,2000.  WE  WILL  USE  BACKGROUND 
PROCESSING  TO  WRITE  THESE  FOUR  VALUES  TO  AN  OUTPUT 
FILE  CALLED  term4.out.  WE  WILL  USE  THESE  FOUR  OUTPUT 
VALUES  FOR  TERM  #4  TO  LATER  ,  IN  ANOTHER  PROGRAM 
COMPUTE  THE  DETECTION  STATISTICS  AND  PROBABILITIES 
OF  BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


7. 

7. 

7. 

7. 


JOHN  A.  STUDER  DATE  LAST  MODIFIED:  13  SEP  94 

CPT,  U.S.  ARMY 

550-53-7181 


format  short  e 
rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

term4sum(i)  =  0; 

for  g  =  betatau(i) : ((2*betatau(i))-l) 
for  m  =  betatau(i) : ((2*betatau(i))-l) 

term4sum(i)  =  term4sum(i)  +  ( (rho" (g+m) )* . . . 

(l-(max(g-betatau(i) ,m-betatau(i))/betatau(i)))) ; 

end 

end 

term4sumf  inal(i)  =  term4s'um(i)  ; 

end 

term4sumf inal 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commands  to  get  started:  intro,  demo,  help  help 

Commaoids  for  more  information:  help,  whatsnew,  info,  subscribe 


term4sumf inal  = 

3.0440e-01 

1 . 5865e“05 

7.2482e-10 

3.2152e-14 

y,  Numbers/arrows/words 
y,  below  were  added  after 
y,  the  mat  lab  bkgd  job 

y. 

•' 

- 

- 

y,  dumped  the  values  to 
y,  the  file  term4.out 

y.  500 

1000 

1500 

2000 

y.  All  values/text  not 

y. 

y. 

V  - 

-  betatau 

values  - 

y,  done  by  Sun  Stn  is 
y,  preceded  by  a  "V' ■ 

•/.•/.•/.•/.•///.•/.•/.•/.•/.•/.mn 

*/•  */. 

y.  TERM  #4  y, 

y.  y. 

y!y.y.y.y.y.y.y.y.y.y.y.y.y.y.y! 


60010008  flops. 
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THESIS  COMPUTER  WORK 


TERM  #5 
betatau  =  500 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #5 
FOR  THE  VALUE  OF  betatau  =  500.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  500  are  I  =  [4  5  8  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #5  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l) .  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term5500.mat 
SINCE  THE  FILE  IS  IN  BINARY  FORMAT  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE 
READABLE  HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  22  OCT  94 


format  short  e 
rho  =  0.99; 
betatau  =  500; 

I  =  [4  5  8  20]  ; 

M  =  [124  98  60  24]  ; 

g  =  0 : (betatau-1) ; 

m  =  0 : (betatau-1) ; 

row  =  rho.“g; 

col  =  rho.“m; 

rhomatrix  =  col'*row; 

rowprime  =  g(ones(betatau,l) , :) ; 

colprime  =  rowprime'; 

WO  =  max (rowprime, colprime) ; 

C  =  WO/betatau; 

ONEMINUSMAX  =  1-C; 
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A  =  colprime-rowprime ;  %  m-g  for  later  complex  exponential 

for  i  =  1:4 

bitmatrix  =  [-M(i)/2:-l  zeros(l,M(i)/2) ;zeros(l,M(i)/2)  l:M(i)/2;... 

-M(i)/2:-l  l:M(i)/2]; 

const  =  ((2*pi*I(i))/betatau) ; 
for  bitpat  =  1:3 

term5sum(i jbitpat)  =  0; 
work  =  bitmatrix(bitpat, :) ; 
m  =  work==0; 
work(m)  =  []  ; 
k  =  work; 
n  =  k; 

for  kct  =  1: length (k) 
for  net  =  1: length (n) 
if  kct  ==  net 


B  =  j*((2*pi*k(kct)*I(i)*A)/betatau) ; 

D  =  exp(B) ; 

term5sum(i , bitpat)  =  term5sum(i ,bitpat)+ . . . 

sum(sum(rhomatrix . *D . *0NEMINUSMAX) ) ; 


else 

El  =  exp(j*const*((n(nct)*colprime)-(k(kct)*rowprime)))  ■J'l.VI, 
E2  =  l-exp(j*const*(k(kct)-n(nct))*WQ) ; 
term5sum(i , bitpat)  =  term5sum(i, bitpat)  +... 
sum(sum( (rhomatrix . *E1 . *E2) / ( j ♦2*pi*I (i) * (k(kct ) -n(nct) ) ) ) ) ; 

end 


end 


end 

end 


end 

save  termSBOO  termBsum 
exit 
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<MATLAB  (R)  > 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 
All  Rights  Reserved 
Version  4.1 
Jun  15  1993 


Coitunands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  term5500.mat 
>>  who 

Your  vauriables  are: 
terms sum 
>>  termSsum 
termBsum  = 


l.Oe+03  * 


1.5423  -  O.OOOOi 
0.9419  -  O.OOOOi 
0.3307  -  O.OOOOi 
0.0431  -  O.OOOOi 


1.5423  +  O.OOOOi 
0.9419  +  O.OOOOi 
0.3307  +  O.OOOOi 
0.0431  +  O.OOOOi 


1.3825  +  O.OOOOi  7.  I  =  4 

0.9001  -  O.OOOOi  y.  I  =  5 

0.3576  +  O.OOOOi  */.  I  =  8 

0.0580  -  O.OOOOi  */.  I  =  20 


» 

y. 

y.  CASE:  (b-=l,b+=0)  CASE:  (b-=0,b+=l) 


CASE:  (b-=l,b+=l) 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y.  TERM  #5  y. 

y,  betatau  =  500  % 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THIS  IS  THE  FILE  term5500.out  CREATED  FROM  THE 

BINARY  FILE  term5500.mat  USING  MATLAB  INTERACTIVE  COMMANDS. 
AT  T,  text  preceded  by  a  y.  WAS  ADDED  LATER  USING  A 
TEXT  EDITOR. 
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•/.  THESIS  COMPUTER  WORK 

•/. 

•/. 

•/.  _ 

•/.  I  I 

'/.  I  TERM  #6  I 

y.  I _ I 

y. 

y,  COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES  OF 
y.  TERM  #6  CORRESPONDING  TO  THE  FOUR  VALUES  OF  betatau  WHICH  ARE: 

y.  500,1000,1500,2000.  WE  WILL  USE  BACKGROUND  PROCESSING  TO 

y.  COMPUTE  THE  VALUES  AND  WRITE  THESE  FOUR  VALUES  TO  A  FILE 

y.  CALLED  terme.out.  WE  WILL  USE  THESE  FOUR  VALUES  FOR  TERM  #6 

y.  TO  LATER,  IN  ANOTHER  PROGRAM  COMPUTE  THE  DETECTION  STATISTICS 

y.  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO 

y.  NOISE  RATIOS . 

y. 

y.  JOHN  A.  STUDER  DATE  LAST  MODIFIED:  16  SEP  94 

y.  CPT,  U.S.  ARMY 
y.  550-53-7181 


format  short  e 
rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

term6sum(i)  =  0; 
for  g  =  l:betatau(i) 

for  m  =  1 :betatau(i)-l 
if  (g-1)  <  m 

term6sum(i)  =  term6sum(i)  +  ((rho''((g-l)+m))*(m-(g-l))) ; 

end 


end 

end 

termGsumf inal(i)  =  (2*term6sum(i) )/betatau(i) ; 
end 

termSsumf inal 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 


term6sumf inal  = 

1.8585e+03 

9.9411e+02 

6.6331e+02 

4.9749e+02 

y,  Numbers/arrows/words 
y,  below  were  added  after 
y,  the  mat  lab  bkgd  job 

y. 

- 

- 

y,  dumped  the  values  to 
y,  the  file  term6.out 

y.  500 

1000 

1500 

2000 

y,  All  values/text  not 

y. 

y.  - 

- betatau 

values  - 

y,  done  by  Sun  Stn  is 
y,  preceded  by  a  . 

•/.  y. 

y.  TERM  #6  y. 

y.  y. 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


33732508  flops. 
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•/.  THESIS  COMPUTER  WORK 

•/. 

•/. 

•/.  _ 

y.  I  I 

y.  I  TERM  #7  I 

y.  I _ I 

y. 

y.  COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  THE  FOUR  VALUES  OF 
y.  TERM  #7  CORRESPONDING  TO  THE  FOUR  VALUES  OF  betatau  WHICH  ARE 

y.  500,1000,1500,2000.  WE  WILL  USE  BACKGROUND  PROCESSING  TO  COMPUTE 

y,  THE  VALUES  AND  WRITE  THESE  FOUR  VALUES  TO  A  FILE  CALLED 

y.  term7.out.  WE  WILL  USE  THESE  FOUR  VALUES  FOR  TERM  #7  TO 

y.  LATER,  IN  ANOTHER  PROGRAM  COMPUTE  THE  DETECTION  STATISTICS 

y,  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO 

y.  NOISE  RATIOS. 

y. 

y.  JOHN  A.  STUDER  DATE  LAST  MODIFIED:  16  SEP  94 

y.  CPT,  U.S.  ARMY 
y.  550-532-7181 


format  short  e 
rho  =  0.99; 

betatau  =  [500  1000  1500  2000] ; 
for  i  =  1:4 

term7sum(i)  =  0; 
for  g  =  l:betatau(i) 

for  m  =  betatau(i) : ((2*betatau(i))-l) 
term7sum(i)  =  term7sum(i)  +  ... 

( (rho" ( (g-l)+m) ) *(1- (max(g-l ,m-betatau(i) ) /betatau (i) ) ) ) ; 

end 

end 

term7sumf inal(i)  =  2*term7sum(i) ; 
end 

term7sumf inal 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 


term7sumf inal 

9.2657e+01 

7.3500e-01 

5.1105e-03 

3.4503e-05 

y,  Numbers/arrows/words 
y,  below  were  added  after 
y,  the  matlab  bkgd  job 

y.  ‘ 

- 

- 

y,  dumped  the  values  to 
y,  the  file  term7.out 

y.  500 

1000 

1500 

2000 

y,  All  values/text  not 

y. 

y. 

y.  - 

- betatau 

values  - 

y,  done  by  Sun  Stn  is 
y,  preceded  by  a  . 

y.  y. 


y,  TERM  #7  y. 

y.  y. 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


67510004  flops. 
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y.  THESIS  COMPUTER  WORK 

y. 

y. 

y. 

y.  _ 

y.  I  I 

y.  I  TERM  #8  I 

y,  I  betatau  =  500  j 

y.  I _ I 

y. 

y,  COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #8 
y.  FOR  THE  VALUE  OF  betatau  =  500.  WE  WILL  COMPUTE 

y,  FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 

y.  betatau  =  500  are  I  =  [4  5  8  20]  .  WE  WILL  COMPUTE  THE  VALUE  OF 

y,  TERM  #8  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 

y,  ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 

y.  BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 

y.  WITH  O^S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 

y.  ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 

y,  (have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 

y.  IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 

y.  LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 

y,  (b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 

y.  DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term8500 .  out .  WE  WILL 

y.  USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 

y.  DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 

y,  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 

y. 

y.  JOHN  A.  STUDER  DATE  LAST  MODIFIED:  26  SEP  94 

y.  CPT,  U.S.  ARMY 
y,  550-53-7181 


format  short  e 
rho  =  0.99; 
betatau  =  500; 

I  =  [4  5  8  20]  ; 

M  =  [124  98  60  24] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term8sum500(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  l:betatau 
for  m  =  1: betatau 

term8sum500(i)  =  term8sum500(i)  +  ( (rho" ((g-l)+(m-l)) )/ . . . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/. . . 
betatau)* (max (g-1 ,m-l) - (m-1) ) )) ; 

end 

end 

if  ct  ==  M(i)/2 
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term85001eft(i)  =  term8sum500(i) ; 

end 


end 

term85001r(i)  =  term8sum500(i) ; 

term8500right (i)  =  term85001r(i)-term85001eft(i) 


end 

term85001eft 

term8500right 

term85001r 

exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commcinds  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 


termSSOOleft  = 

-3.7098e+02  -2.3852e+02  -9.2820e+01  -1.4109e+01  ’/.  CASE:  (b-=l,b+=0) 
term8500right  = 

-3.7098e+02  -2.3852e+02  -9.2820e+01  -1.4109e+01  */.  CASE:  (b-=0,b+=l) 


term85001r  = 

-7.4196e+02  -4.7703e+02  -1.8564e+02  -2.8217e+01  I  CASE:  (b-=l,b+=l) 

•/.  ' 

'/.  1  =  4  1  =  5  1  =  8  1  =  20 

y,  - 

•/. 

'/.  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


y.  TERM  #8  y. 

y,  betatau  =  500  */, 

y.y//.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  ALL  TEXT  PRECEDED  BY  A  */,  WAS  ADDED 

BY  A  TEXT  EDITOR  AFTER  THE  BACKGROUND 
JOB  CREATED  THE  FILE  term8500 . out . 


2142000318  flops. 
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THESIS  COMPUTER  WORK 


•/. 

y. 

•/. 

•/. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


TERM  #8 

betatau  =  1000 


I. 


.1 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #8 

FOR  THE  VALUE  OF  betatau  =  1000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  1000  are  I  =  [569  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #8  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l) .  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  termSlk.mat 
SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT  WE  WILL  THEN  USE 
MATLAB  INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE 
READABLE  HARDCOPY  RESULTS.  AFTER  THIS  IS  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 

/  •/  V  •/  •/  •/•/•/  •/  •/ »/  •/  •/  •/ 


DATE  LAST  MODIFIED:  27  SEP  94 


format  short  e 
rho  =  0.99; 
betatau  =  1000; 

I  =  [5  6  9  20]  ; 

M  =  [198  164  110  48] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term8sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  1: betatau 
for  m  =  l:betatau 

term8sum(i)  =  term8sum(i)  +  ((rho“((g-l)+(m-l)))/. . . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l) /betatau) ))-sin(((2*pi*k(ct)*I(i))/. . . 
betatau)* (max(g-l,m-l)-(m-l)))) ; 


end 


end 

if  ct  ==  M(i)/2 

termSlkleft (i)  =  term8sum(i); 

end 

end 

term81klr(i)  =  term8sum(i) ; 

term81kright (i)  =  term81klr(i)-term81kleft(i) ; 

end 

save  term81k  term81kleft  term81kright  term81klr 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks, 

All  Rights  Reserved 
Version  4.1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 
Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  termSlk 
»  who 

Your  variables  are: 

termSlkleft  termSlklr  termSlkright 

»  termSlkleft 
termSlkleft  = 

-464.3344  -327.7899  -148.2912  -29.7545  */.  CASE:  (b-=l,b+=0) 

>>  termSlkright 
termSlkright  = 

-464.3344  -327.7899  -148.2912  -29.7545  */.  CASE:  (b-=0,b+=l) 

>>  termSlklr 
termSlklr  = 

-928.6688  -655.5798  -296.5824  -59.5090  V.  CASE:  (b-=l,b+=l) 

» 

y.  - 

y.  1  =  5  1  =  6  1  =  9  1  =  20 

y.  _ 

y. 

y.  VALUES  OF  THE  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y.  TERM  #8  y. 
y,  betatau  =  1000  */, 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THIS  IS  THE  FILE  termSlk. out  CREATED  FROM  THE 

BINARY  FILE  termSlk. mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  I  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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•/.  THESIS  COMPUTER  WORK 

•/. 

•/. 

•/. 

•/.  _ 

y.  I  I 

y.  I  TERM  #8  I 

y,  I  betatau  =  1500  I 

y.  I _ I 

y. 

y.  COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #8 

y.  FOR  THE  VALUE  OF  betatau  =  1500.  WE  WILL  COMPUTE 

y.  FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER)  .  THE  VALUES  OF  I  FOR 

y.  betatau  =  1500  are  I  =  [7  9  12  20]  .  WE  WILL  COMPUTE  THE  VALUE  OF 

y.  TERM  #8  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 

y.  ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 

y.  BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 

y.  WITH  O'S.  THIS  IS  THE  CASE;  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 

y.  ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 

y,  (have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 

y.  IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 

y.  LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 

y.  (b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 

y.  DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term815k.out  WE  WILL 

y.  USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 

y.  DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 

y.  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 

y. 

y.  JOHN  A.  STUDER  DATE  LAST  MODIFIED:  27  SEP  94 

y.  CPT,  U.S.  ARMY 
y.  550-53-7181 


format  short  e 
rho  =  0.99; 
betatau  =  1500; 

I  =  [7  9  12  20]; 

M  =  [212  164  124  74] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term8sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  1: betatau 
for  m  =  1: betatau 

term8sum(i)  =  term8sum(i)  +  ((rho''((g-l)  +  (m-l)))/.  .  . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/. . . 
betatau)* (max (g-1 ,m-l)-(m-l)) )) ; 

end 

end 

if  ct  ==  M(i)/2 
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terms IBkleft (i)  =  term8sum(i) ; 


end 

end 

term815klr(i)  =  term8sum(i) ; 

termSlBkright (i)  =  term815klr(i)-term815kleft(i) ; 

end 

termSlSkleft 

termSlSkright 

term815klr 

exit 
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<MATLAB(R)> 

(c)  Copyright  1984-94  The  MathWorks ,  Inc. 

All  Rights  Reserved 
Version  4.2 
Mar  29  1994 

Commands  to  get  started;  intro,  demo,  help  help 

Commands  for  more  information;  help,  whatsnew,  info,  subscribe 

malloc  count  matches  malloc_debug_count 
malloc  count  matches  malloc_debug_count 
malloc  count  matches  malloc_debug_count 
malloc  count  matches  malloc_debug_count 
malloc  count  matches  malloc_debug_count 
malloc  count  matches  malloc_debug_count 

terms 15kleft  = 

-3.5253e+02  -2.1854e+02  -1.2475e+02  -4.5096e+01  ’/.  CASE;  (b-=l,b+=0) 


termSlSkright  = 

-3.5253e+02  -2.1854e+02  -1.2475e+02  -4.5096e+01  */.  CASE;  (b-=0,b+=l) 


termSlSklr  = 


-7.0505e+02  -4.3708e+02  -2.4950e+02  -9.0192e+01  ’/.  CASE;  (b-=l,b+=l) 

•/.  “ 

•/.  1  =  7  1  =  9  I  =  12  I  =  20 

1  _ 

•/. 

*/.  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


•/.  TERM  #8  •/. 

*!,  betatau  =  1500  */, 


NOTE;  THIS  IS  THE  FILE  term815k.out .  ALL  TEXT 
PRECEDED  BY  A  7.  WAS  ADDED  LATER  USING  A 
TEXT  EDITOR. 


36162000586  flops. 
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THESIS  COMPUTER  WORK 


TERM  #8 

betatau  =  2000 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #8 

FOR  THE  VALUE  OF  betatau  =  2000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  2000  are  I  =  [8  9  12  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #8  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  siimmation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I’S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term82k.mat.  WE  WILL 
THEN  SINCE  THE  FILE  IS  IN  BINARY  FORMAT  USE  MATLAB  INTERACTIVE 
COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE  HARDCOPY 
RESULTS.  WE  WILL  USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM 
TO  COMPUTE  THE  DETECTION  STATISTICS  AND  PROBABILITIES  OF 
BIT  ERROR  FOR  THE  VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  27  SEP  94 


format  short  e 
rho  =  0.99; 
betatau  =  2000; 

I  =  [8  9  12  20]; 

M  =  [248  220  164  98] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term8sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  1: betatau 
for  m  =  l:betatau 

term8sum(i)  =  term8sum(i)  +  ((rho“((g-l)+(m-l)))/. . . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*.  .  . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/ . . . 
betatau)*(max(g-l,m-l)-(m-l) ))) ; 

end 
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end 

if  ct  ==  M(i)/2 

tenn82kleft(i)  =  term8suni(i)  ; 

end 

end 

term82klr(i)  =  term8sum(i); 

term82kright(i)  =  term82klr(i)-term82kleft(i) ; 

end 

save  term82k  term82kleft  term82kright  term82klr 
exit 
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<MATLAB  (R)  > 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4 . 1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  term82k.mat 
>>  who 

Your  variables  are: 

term82kleft  term82klr  term82kright 

>>  term82kleft 
term82kleft  = 

-351.8931  -283.0055  -163.9044  -60.1542  */.  CASE:  (b-=l,b+=0) 

>>  term82kright 
term82kright  = 

-351.8931  -283.0055  -163.9044  -60.1542  */.  CASE:  (b-=0,b+=l) 

>>  term82klr 
term82klr  = 

-703.7863  -566.0110  -327.8088  -120.3084  */.  CASE:  (b-=l,b+=l) 

•/.  I  =  8  1  =  9  I  =  12  I  =  20 

•/.  _ 

•/. 

•/.  VALUES  OF  THE  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


•/.7.y.y.y.y.y.y.7.y.y.y.7.7.*/.y.*/.’/. 

y.  TERM  #8  y, 

y  betatau  =  2000  7, 

y.y.y.y.y.y.y.7.y.y.y.y.y.y.y.y.7.y. 


NOTE:  THIS  IS  THE  FILE  term82k.out  CREATED  FROM  THE 

BINARY  FILE  term82k.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  7,  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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THESIS  COMPUTER  WORK 


•/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


TERM  #9  I 

betatau  =  500  | 


COMMENTS;  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #9 
FOR  THE  VALUE  OF  betatau  =  500.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  500  are  I  =  [4  5  8  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #9  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE; 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term9500 .mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE  WE  WILL 
USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 

JOHN  A.  STUDER  DATE  LAST  MODIFIED:  28  SEP  94 

CPT,  U.S.  ARMY 

550-53-7181 


format  short  e 
rho  =  0.99; 
betatau  =  500; 

I  =  [4  5  8  20]  ; 

M  =  [124  98  60  24] ; 
for  i  =  1:4 

k  =  [-M(i)/2;-l  l:M(i)/2] ; 
term9sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  l:betatau-l 
for  m  =  1: betatau 
if  m-1  <  g 

term9sum(i)  =  term9sum(i)  +  ((rho~(g+(m-l)))/(pi*k(ct) . . . 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l)) . . . 
/betatau)) ; 
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end 


end 

end 

if  ct  ==  M(i)/2 

term95001eft (i)  =  term9snm(i) ; 

end 

end 

term95001r(i)  =  term9siim(i) ; 

term9500right (i)  =  term95001r(i)-term95001eft(i) ; 

end 

save  term9500  term95001eft  term9500right  term95001r 
exit 
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<108  sp254201 (SunOS)  /kepler_u2/studer>  matlab 

<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  term9500.mat 
»  who 

Your  variables  are: 

term95001eft  term9500right 

term95001r 

>>  term95001eft 

term95001eft  = 

125.1608  80.4704  31.3156  4.7599  '/.  CASE:  (b-=l,b+=0) 

»  term9500right 
term9500right  = 

125.1608  80.4704  31.3156  4.7599  '/.  CASE:  (b-=0,b+=l) 

»  term95001r 
term95001r  = 

250.3216  160.9408  62.6313  9.5199  '/,  CASE:  (b-=l,b+=l) 

» 

•/.  1=4  1  =  5  1  =  8  1  =  20 

•/.  _ 

•/. 

y.  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


NOTE:  THIS  IS  THE  FILE  term9500.out  CREATED  FROM  THE  BINARY 

FILE  term9500.mat  USING  MATLAB  INTERACTIVE  COMMANDS.  ALL 
TEXT  PRECEDED  BY  A  '/.  WAS  ADDED  LATER  USING  A  TEXT 
EDITOR. 
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THESIS  COMPUTER  WORK 


TERM  #9  I 

betatau  =  1000  j 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #9 

FOR  THE  VALUE  OF  betatau  =  1000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  1000  are  I  =  [5  6  9  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #9  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l.b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term91k.mat. 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE  WE  WILL 
USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  28  SEP  94 


format  short  e 
rho  =  0.99; 
betatau  =  1000; 

I  =  [5  6  9  20]  ; 

M  =  [198  164  110  48]; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term9sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  l:betatau-l 
for  m  =  1: betatau 
if  m-1  <  g 

term9sum(i)  =  • 


term9sum(i)  +  ((rho'‘(g+(m-l)))/(pi*k(ct)  . 
*I(i)) ) *sin(2*pi*k(ct)*I(i)*((g-(m-l)) 
/betatau)) ; 
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end 

end 

end 

if  ct  ==  M(i)/2 

terni91kleft(i)  =  term9sum(i) ; 

end 

end 

term91klr(i)  =  term9sum(i) ; 

term91kright(i)  =  term91klr(i)-terra91kleft(i) ; 

end 

save  term91k  term91kleft  term91kright  term91klr 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4. 1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

»  load  term91k.mat 
»  who 

Your  variables  are: 

term91kleft  term91klr  term91kright 

»  term91kleft 


term91kleft  = 


155.3047  109.6350 

49.5986 

9.9519 

y.  CASE:  (b-=l,b+=0) 

»  term91kright 

termOlkright  = 

155.3047  109.6350 

49.5986 

9.9519 

y,  CASE:  (b-=0,b+=l) 

»  term91klr 

term91klr  = 

310.6094  219.2701 

99.1971 

19.9038 

y.  CASE:  (b-=l,b+=l) 

» 

y. 

y.  1  =  5  1  =  6 

y.  ___ 

1  =  9 

I  =  20 

•/.  VALUES  OF  THE  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


y;/.y.y.r/.*/.y.y.y//.y.y.y.’/.y.7.’/. 

y.  TERM  #9  y. 
%  betatau  =  1000  */, 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THIS  IS  THE  FILE  term91k.out  CREATED  FROM  THE 

BINARY  FILE  term91k.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  */.  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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THESIS  COMPUTER  WORK 


TERM  #9 

betatau  =  1500 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #9 

FOR  THE  VALUE  OF  betatau  =  1500.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  1500  are  I  =  [7  9  12  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #9  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O’S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term915k.mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE  WE  WILL 
USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  28  SEP  94 


format  short  e 
rho  =  0.99; 
betatau  =  1500; 

I  =  [7  9  12  20] ; 

M  =  [212  164  124  74]  ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term9sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  l:betatau-l 
for  m  =  l:betatau 
if  m-1  <  g 

term9sum(i)  =  term9sum(i)  +  ((rho~ (g+(m-l) ) )/(pi*k(ct) . . . 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l)) . . . 
/betatau)) ; 
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end 


end 

end 

if  ct  ==  M(i)/2 

term915kleft  (i)  =  term9sTim(i)  ; 

end 

end 

term915klr(i)  =  term9sum(i) ; 

term915kright(i)  =  term915klr(i)-term915kleft(i) 

end 

save  term915k  term915kleft  term915kright  term915klr 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  term915k 
>>  who 

Your  variables  are: 

term915kleft  term915kright 

term915klr 

>>  term915kleft 


term915kleft  = 

117.9019  73.0900  41.7221  15.0822  */.  CASE:  (b-=l,b+=0) 

>>  term915kright 
term915kright  = 

117.9019  73.0900  41.7221  15.0822  */.  CASE:  (b-=0,b+=l) 

>>  term915klr 
term915klr  = 

235.8037  146.1800  83.4441  30.1644  */,  CASE:  (b-=l,b+=l) 

*/.  “ 

•/,  1  =  7  1  =  9  I  =  12  I  =  20 

•/.  _ 

•/. 

y.  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


7.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.7.y.y. 

7.  TERM  #9  7. 

7,  betatau  =  1500  7, 

y.7.y.7x/.y.y.y.y.y.y.y.7.y.y.7.y.y. 


NOTE:  THIS  IS  THE  FILE  term915k.out  CREATED  FROM  THE 

BINARY  FILE  term915k.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  7,  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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THESIS  COMPUTER  WORK 


TERM  #9  I 

betatau  =  2000  j 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #9 

FOR  THE  VALUE  OF  betatau  =  2000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  2000  are  I  =  [8  9  12  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #9  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I’S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  term92k.mat. 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE  WE  WILL 
USE  THESE  VALUES  LATER,  IN  ANOTHER  PROGRAM  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  28  SEP  94 


NNNNNNNNNNNNNNNNNNNNNNNNMNNNUNNNNNNNNNNNNNNNNNNNNNNNNNNNNUNNNNNNUNNNNNNNUNNNNk 


format  short  e 
rho  =  0.99; 
betatau  =  1500; 

I  =  [8  9  12  20]; 

M  =  [248  220  164  98] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
term9sum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  l:betatau-l 
for  m  =  1: betatau 
if  m-1  <  g 

term9sum(i)  =  term9sum(i)  +  ((rho“(g+(m-l)))/(pi*k(ct) . . . 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l)) . . . 
/betatau)) ; 


164 


end 

end 

end 

if  ct  ==  M(i)/2 

term92kleft (i)  =  term9sum(i) ; 

end 

end 

term92klr(i)  =  term9sum(i) ; 

term92kright(i)  =  terin92klr(i)-term92kleft(i)  ; 

end 

save  term92k  term92kleft  term92kright  term92klr 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4 . 1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  term92k.mat 
»  who 

Your  variables  are: 

term92kleft  term92klr  term92kright 

»  term92kleft 


term92kleft  = 

91.5258  73.0261  41.6758  15.0544  •/.  CASE:  (b-=l,b+=0) 

»  term92kright 
term92kright  = 

91.5258  73.0261  41.6758  15.0544  •/.  CASE:  (b-=0,b+=l) 

»  term92klr 


term92klr  = 

183.0515  146.0523  83.3516  30.1088  '/,  CASE:  (b-=l,b+=l) 


» 

•/. 

•/. 

•/. 

•/. 

•/. 


1  =  8  1  =  9  1  =  12  1  =  20 


VALUES  OF  THE  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


y;/x///xmy.y.y//.*/.y.y.y.y. 

y.  TERM  #9  y. 
y,  betatau  =  2000  '/. 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THIS  IS  THE  FILE  term92k.out  CREATED  FROM  THE 

BINARY  FILE  term92k.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  %  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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•/. 

•/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


THESIS  COMPUTER  WORK 


I  TERM  #10  I 

1  betatau  =  500  j 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #10 
FOR  THE  VALUE  OF  betatau  =  500.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  500  are  I  =  [4  5  8  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #10  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE;  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  terml0500 .mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT,  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM,  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


y. 

y,  JOHN  A.  STUDER  DATE  LAST  MODIFIED:  29  SEP  94 

y.  CPT,  U.S.  ARMY 
y.  550-53-7181 


format  short  e 
rho  =  0.99; 
betatau  =  500; 

I  =  [4  5  8  20]  ; 

M  =  [124  98  60  24] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
termlOsum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  betatau : (2*betatau) -1 
for  m  =  l:betatau 

termlOsum(i)  =  termlOsum(i)  +  ((rho~(g+(m-l)))/. . . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/ . . . 
betatau)* (max(g-betatau,m-l)-(m-l) )) ) ; 
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end 

end 

if  ct  ==  M(i)/2 

tennl05001eft(i)  =  termlOsumCi) ; 

end 

end 

terml05001r(i)  =  termlOsuin(i) ; 

termlOSOOright (i)  =  terml05001r(i)-terml05001eft(i) ; 

end 

save  termlOSOO  termlOSOOleft  termlOSOOright  termlOSOOlr 
exit 
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<MATLAB(R)> 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4.1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  termlOBOO .mat 
>>  who 

Your  variables  are: 


termlOSOOleft  termlOSOOright 

termlOSOOlr 


>>  termlOBOOleft 
termlOBOOleft  = 

-2.4375  -1.5672  -0.6099  -0.0927  '/,  CASE:  (b-=l,b+=0) 

>>  termlOBOOright 
terml0500right  = 

-2.4375  -1.5672  -0.6099  -0.0927  */.  CASE:  (b-=0,b+=l) 

>>  terml05001r 
termlOSOOlr  = 

-4.8750  -3.1343  -1.2197  -0.1854  */.  CASE:  (b-=l,b+=l) 


» 

•/.  1  =  4  1  =  5  1  =  8  1  =  20 

•/.  _ 

•/. 

•/,  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 


’/.•/.•/.7.7.y.*/.y.y.*/.y.y.7.y.y.7.y.y.y. 

y.  TERM  #10  y. 


y,  betatau  =  500  7, 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THIS  IS  THE  FILE  termlOBOO .out  CREATED  FROM  THE  BINARY  FILE 
termlOBOO. mat  USING  MATLAB  INTERACTIVE  COMMANDS.  ALL  TEXT 
PRECEDED  BY  A  7.  WAS  ADDED  LATER  USING  A  TEXT  EDITOR. 
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•/. 

•/. 

•/. 

•/. 

•/. 

•/. 

*/. 

•/. 

•/. 

•/. 

•/. 

•/. 

•/. 

•/. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


THESIS  COMPUTER  WORK 


TERM  #10  I 
betatau  =  1000  j 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #10 
FOR  THE  VALUE  OF  betatau  =  1000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  1000  are  I  =  [5  6  9  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #10  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  termlOlk.mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT,  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM,  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 

JOHN  A.  STUDER  DATE  LAST  MODIFIED:  29  SEP  94 

CPT,  U.S.  ARMY 

550-53-7181 


format  short  e 
rho  =  0.99; 
betatau  =  1000; 

I  =  [5  6  9  20]  ; 

M  =  [198  164  110  48]  ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
termlOsum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  betatau: (2*betatau)-l 
for  m  =  1: betatau 

termlOsum(i)  =  termlOsum(i)  +  ((rho'‘(g+(m-l)))/.  .  . 

(pi>t'k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*.  . . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/. . . 
betatau) * (maLx(g-betatau, m-1)- (m-1) )) )  ; 
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end 


end 

if  ct  ==  M(i)/2 

termlOlkleft (i)  =  termiOsum(i) ; 

end 

end 

termlOlklr(i)  =  termlOsum(i) ; 

terralOlkright (i)  =  termlOlklr(i) -termlOlkleft (i) ; 

end 

save  termlOlk  termlOlkleft  termlOlkright  termlOlklr 
exit 
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(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4 . 1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  termlOlk.mat 
>>  who 

Your  variables  are: 

termlOlkleft  termlOlkright 

termlOlklr 

>>  termlOlkleft 

termlOlkleft  = 

-0.0200  -0.0142  -0.0064  -0.0013  */.  CASE:  (b-=l,b+=0) 

>>  termlOlkright 
termlOlkright  = 

-0.0200  -0.0142  -0.0064  -0.0013  7.  CASE:  (b-=0,b+=l) 

>>  termlOlklr 
termlOlklr  = 


NOTE:  THIS  IS  THE  FILE  termlOlk.out  CREATED  FROM  THE 

BINARY  FILE  termlOl.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  7.  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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THESIS  COMPUTER  WORK 


TERM  #10 
betatau  =  1500 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #10 
FOR  THE  VALUE  OF  betatau  =  1500.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  1500  are  I  =  [7  9  12  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #10  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O'S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I’S.  THIS 
IS  THE  CASE:  (b-=0,b+=l).  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  terml015k.mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT,  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM,  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 
550-53-7181 


DATE  LAST  MODIFIED:  29  SEP  94 


format  short  e 
rho  =  0.99; 
betatau  =  1500; 

I  =  [7  9  12  20] ; 

M  =  [212  164  124  74] ; 
for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
termlOsum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  betatau: (2*betatau) -1 
for  m  =  l:betatau 

termlOsum(i)  =  termlOsum(i)  +  ((rho'‘(g+(m-l)))/.  .  . 

(pi*k(ct)*l(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l) /betatau) ))-sin(((2*pi*k(ct)*I(i))/ . . . 
betatau)*(max(g-betatau,m-l)-(m-l) ) ) ) ; 

end 

end 
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if  ct  ==  M(i)/2 

terml015kleft (i)  =  termlOsum(i) ; 

end 

end 

terml015klr(i)  =  termlOsum(i) ; 

terinl015kright(i)  =  tennl015klr(i)-ternil015kleft (i) ; 

end 

save  termlOlSk  termlOlBkleft  termlOlBkright  termlOlBklr 
exit 
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<MATLAB  (R)  > 

(c)  Copyright  1984-93  The  MathWorks,  Inc. 

All  Rights  Reserved 
Version  4 . 1 
Jun  15  1993 

Commands  to  get  started:  intro,  demo,  help  help 

Commands  for  more  information:  help,  whatsnew,  info,  subscribe 

>>  load  termlOlSk.mat 
>>  who 

Your  variables  are: 

termlOlSkleft  termlOlSkright 

terml015klr 

>>  termlOlBkleft 

termlOlBkleft  = 

l.Oe-04  * 

-1.0000  -0.6199  -0.3539  -0.1279  ’/.  CASE:  (b-=l,b+=0) 

>>  termlOlBkright 
termlOlBkright  = 
l.Oe-04  * 

-1.0000  -0.6199  -0.3539  -0.1279  ’/,  CASE:  (b-=0,b+=l) 

>>  termlOlSklr 
termlOlBklr  = 
l.Oe-03  * 

-0.2000  -0.1240  -0.0708  -0.0256  7.  CASE:  (b-=l,b+=l) 

» 

7.  “ 

7.  1  =  7  1  =  9  I  =  12  I  =  20 

7.  _ 

7. 

7.  VALUES  OF  NORMALIZED  CHANNEL  SPACING  INTEGER  I 
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NOTE: 


y.y//.‘/.‘/.y.y.r/.y.y.y////.*/.y.y.y.y. 

y.  TERM  #10  y. 

y,  betatau  =  1500  */, 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


THIS  IS  THE  FILE  termlOlSk. out  CREATED  FROM  THE 
BINARY  FILE  termlOlSk.mat  USING  MATLAB  INTERACTIVE 
COMMANDS.  ALL  TEXT  PRECEDED  BY  A  */.  WAS  ADDED  LATER 
USING  A  TEXT  EDITOR. 
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•/. 

1 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


TERM  #10 
betatau  =  2000 


I, 


COMMENTS:  THE  IDEA  OF  THIS  PROGRAM  IS  TO  COMPUTE  TERM  #10 
FOR  THE  VALUE  OF  betatau  =  2000.  WE  WILL  COMPUTE 
FOR  THE  FOUR  VALUES  OF  I  (AN  INTEGER) .  THE  VALUES  OF  I  FOR 
betatau  =  2000  are  I  =  [8  9  12  20] .  WE  WILL  COMPUTE  THE  VALUE  OF 
TERM  #10  THREE  TIMES  FOR  EACH  CORRESPONDING  VALUE  OF  I.  ONCE  FOR 
ALL  M/2  OF  THE  LOWER  ACI  CHANNELS  (have  negative  summation  index) 
BEING  PACKED  WITH  I'S  AND  THE  UPPER  ACI  CHANNELS  ALL  BEING  PACKED 
WITH  O’S.  THIS  IS  THE  CASE:  (b-=l,b+=0).  ONCE  FOR  ALL  THE  LOWER 
ACI  CHANNELS  BEING  PACKED  WITH  ZERO'S  AND  THE  UPPER  ACI  CHANNELS 
(have  a  positive  summation  index)  BEING  PACKED  WITH  I'S.  THIS 
IS  THE  CASE:  (b-=0,b+=l) .  FINALLY  ONCE  FOR  BOTH  THE  UPPER  AND 
LOWER  ACI  CHANNELS  BEING  PACKED  WITH  I'S.  THIS  IS  THE  CASE: 
(b-=l,b+=l).  WE  WILL  USE  BACKGROUND  PROCESSING  WHICH  WILL 
DUMP  ALL  TWELVE  VALUES  TO  A  FILE  CALLED  terml02k.mat . 

SINCE  THE  FILE  IS  IN  A  BINARY  FORMAT,  WE  WILL  USE  MATLAB 
INTERACTIVE  COMMANDS  AND  THE  TEXT  EDITOR  TO  PRODUCE  READABLE 
HARDCOPY  RESULTS.  AFTER  THIS  IS  ALL  DONE,  WE  WILL 
USE  THESE  VALUES  LATER  ,  IN  ANOTHER  PROGRAM,  TO  COMPUTE  THE 
DETECTION  STATISTICS  AND  PROBABILITIES  OF  BIT  ERROR  FOR  THE 
VARIOUS  SIGNAL  TO  NOISE  RATIOS. 


JOHN  A.  STUDER 
CPT,  U.S.  ARMY 


DATE  LAST  MODIFIED:  28  NOV  94 


y.  550-53-7181 
format  short  e 


rho  =  0.99; 
betatau  =  2000; 

I  =  [8  9  12  20]; 

M  =  [248  220  164  98]; 


for  i  =  1:4 

k  =  [-M(i)/2:-l  l:M(i)/2] ; 
termlOsum(i)  =  0; 
for  ct  =  l:M(i) 

for  g  =  betatau : (2*betatau) -1 
for  m  =  l:betatau 

termlOsum(i)  =  termlOsum(i)  +  ((rho“(g+(m-l)))/. . . 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*. . . 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/. . . 
betatau) *(max(g-betatau,m-l)-(m-l)))) ; 


end 


end 

if  ct  ==  M(i)/2 

tennl02kleft (i)  =  termlOsum(i) ; 

end 

end 

terml02klr(i)  =  termlOsum(i) ; 

terml02kright(i)  =  terml02klr(i)-terral02kleft(i) ; 

end 

save  terinl02k  terml02kleft  terml02kright  terinl02klr 
exit 
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>>  terml02kleft8 


terml02kleft8  = 

-6.5584e-07  I  CASE:  (b-=l,b+=0) 

>>  terml02kright8 
terinl02kright8  = 

-6.5584e-07  ’/,  CASE:  (b-=0,b+=l) 

»  terml02klr8 
terml02klr8  = 

-1.3117e-06  •/.  CASE:  (b-=l,b+=l) 

» 

>>  terml02kleft9 
ternil02kleft9  = 

-5.2745e-07  */,  CASE:  (b-=l,b+=0) 

>>  terml02kright9 
terml02kright9  = 

-5.2745e-07  •/.  CASE:  (b-=0,b+=l) 

>>  terml02klr9 
ternil02klr9  = 

-1.0549e-06  '/,  CASE:  (b-=l,b+=l) 

» 

>>  terml02kleftl2 
terml02kleftl2  = 

-3.0548e-07 
»  terml02krightl2 
terinl02krightl2  = 


•/.  CASE:  (b-=l,b+=0) 


vixvmva 

*/.  I  =  12  •/. 


vmivm 

•/.  I  =  9  '/. 


•/.  I  =  8  •/. 

vixmxvi. 


-3.0548e-07 


•/.  CASE:  (b-=0,b+=l) 


>>  terml02klrl2 


terml02klrl2  = 

-6.1096e-07  */,  CASE:  (b-=l,b+=l) 

» 

»  terml02kleft20 
terml02kleft20  = 

-1.1211e-07 
»  terml02kright20 
terml02kright20  = 

-1.1211e-07 
»  terml02klr20 
terml02klr20  = 

-2.2423e-07  ’/.  CASE;  (b-=l,b+=l) 

» 


•/.  CASE:  (b-=l,b+=0) 


y,  CASE:  (b-=0,b+=l) 


y.  I  =  20  y. 

y.y.y.y.y.y.y.y.y.y. 


y.  TERM  #10  y. 

y,  betatau  =  2000  */. 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE:  THERE  WERE  SEVERAL  MISHAPS  IN  COMPUTING  THIS  TERM.  POWER  OUTAGES, 

MY  ACCOUNT  BEING  SHUT  OFF, ETC.  IN  THE  INTERESTS  OF  BEING  ABLE  TO 
COMPUTE  RESULTS  IN  A  TIMELY  MANNER,  FOUR  COPIES  OF  THE  PROGRAM  FOR 
TERM  #10,  betatau  =  2000  WERE  MADE.  THE  I  AND  M  BLOCK  WERE  MODIFIED 
TO  HOLD  ONLY  ONE  VALUE  OF  I/M,  AND  THE  OUTER  LOOP  WAS  REDUCED  TO 
ONE  ITERATION.  THEN,  THE  FOUR  PROGRAMS  WERE  RUN  ON  FOUR  SEPARATE 
WORKSTATIONS.  IN  THE  INTERESTS  OF  BREVITY  WE  ONLY  SHOW  THE  SINGLE 
PROGRAM  FOR  COMPUTING  THE  TERM  AS  THE  FOUR  COPIES  ARE  ESSENTIALLY 
THE  EXACT  SAME  PROGRAM  BUT  THEY  DO  ONLY  ONE  OUTER  LOOP  EACH. 

THE  RESULTS  OF  THE  FOUR  PROGRAMS  ARE  TABULATED  HERE  USING  THE 
A  TEXT  EDITOR. 
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THESIS  COMPUTER  WORK 

GENERATION  OF  GRAPHS 
FOR 

betatau  =  500 

DATE:  19  OCTOBER  1994 

COMMENTS:  THIS  PROGRAM  CONTAINS  THE  RESULTS  OF  ALL  TEN  TERMS  OF 
APPENDIX  B  RESULTING  IN  THE  PROGRAM  PRODUCING  A 
COMPLETE,  ACCURATE  GRAPH  OF  THE  SYSTEM'S  PROBABILITY 
OF  BIT  ERROR  FOR  ALL  FOUR  VALUES  OF  I. 

I  =  [4  5  8  20] 


mNNNNNNNNNNNNMNNNNNNNNNNNNNNNNNmNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNm 


format  short  e 
rho  =  0.99; 
c  =  (l-rho)'2; 

I  =  [4  5  8  20]  ; 

M  =  [124  98  60  24] ; 

bitmatrix  =[000;001;010;011;  100;  101;  110;  111]; 

terms  =  [l.5423e+03  0.9419e+03  0.3307e+03  0.0431e+03 

1.5423e+03  0.9419e+03  0.3307e+03  0.0431e+03 

1.3825e+03  0.9001e+03  0.3576e+03  0.0580g+03]; 

term8  =  [-3.7098e+02  -2.3852e+02  -9.2820e+01  -1.4109e+01 

-3.7098e+02  -2.3852e+02  -9.2820e+01  -1.4109e+01 
-7.4196e+02  -4.7703e+02  -1.8564e+02  -2.8217e+0l]  ; 

term9  =  [125.1608  80.4704  31.3156  4.7599 
125.1608  80.4704  31.3156  4.7599 
250.3216  160.9408  62.6313  9.5199]; 

termlO  =  [-2.4375  -1.5672  -0.6099  -0.0927 
-2.4375  -1.5672  -0.6099  -0.0927 
-4.8750  -3.1343  -1.2197  -0.1854]; 

for  ct  =  1:4 
for  i  =  1:8 

for  boo  =  [1  0] 

XI  =  boo*c*7.0510e+03; 

X2  =  bitmatrix(i, 3) *0*959.5662; 

X3  =  bitmatrix(i ,3) *0*1 . 2211e+01 ; 

X4  =  bitmatrix(i ,3)*o*3 . 0440e-01 ; 
if  bitmatrix(i, 1)==1  k  bitmatrix (i, 2) ==0 
X5  =  o*term5(l,ot) ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrixCi ,2)==1 
X5  =  o*term5(2,ot) ; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2) ==1 
X5  =  o*term5C3,ot) ; 


y. 

•/. 

•/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


else 

X5  =  0; 

end 

X6  =  boo*bitmatrix(i,3)*c*l .8585e+03; 

X7  =  boo*bitmatrix(i,3)*c*9.2657e+01; 
if  boo  ==  0 
X8  =  0; 

elseif  bitmatrix(i , 1)==1  &  bitmatrix(i ,2)==0 
X8  =  c*term8(l,ct) ; 

elseif  bitmatrix(i,l)==0  &  bitmatrix(i ,2)==1 
X8  =  c*term8(2, ct) ; 

elseif  bitmatrix(i , 1)==1  &  bitmatrix(i ,2)==1 
X8  =  c*term8(3,ct) ; 
else 

X8  =  0; 

end 

if  bitmatrix(i ,3)  ==  0 
X9  =  0; 

elseif  bitmatrixCi , 1)==1  &  bitmatrix(i ,2)==0 
X9  =  c*term9(l , ct) ; 

elseif  bitmatrix(i , 1)==0  &  bitmatrix(i ,2)==1 
X9  =  c=t'term9(2,ct) ; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==1 
X9  =  c*term9C3,ct) ; 
else 

X9  =  0; 

end 

if  bitmatrixCi ,3)  ==0 
XIO  =  0; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==0 
XIO  =  c*termlOCl ,ct) ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrixCi ,2)==1 
XIO  =  c*terml0C2 ,ct) ; 

elseif  bitmatrixCi , l)==l  &  bitmatrixCi ,2)==1 
XIO  =  c*termlOC3,ct) ; 
else 

XIO  =  0; 

end 

INTERSUM  =  X1+X2+X3+X4+X5+X6+X7+X8+X9+X10; 
if  boo==l 

XlBETACct.i)  =  INTERSUM; 
else 

XOBETACct,i)  =  INTERSUM; 

end 

end 

end 

XlMINCct)  =  minCXlBETACct, ;)); 

XOMAXCct)  =  maxCXOBETACct, :)) ; 

VTCct)  =  CXOMAXCct)  +  XlMINCct))/2; 
end 

RPSQR_TDIVNO_DB  =  10:. 01: 20; 
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RPSQR.TDIVNO  =  10 . " (RPSQR_TDIVNO_DB*0 . 10) ; 

SINGCHAN  =  0.5*erfc(RPSQR_TDIVN0/8''0.5)  ; 
for  ct  =  1:4 

PE(ct, 1:1001)  =  zerosd.lOOl) ; 
for  i  =  1:8 

PE(ct, :)=PE(ct, :)+0.25*erfc((RPSQR_TDIVN0/2“0.5)*(XlBETA(ct,i)-VT(ct))) . .  . 
+0 . 25*erf c ( (RPSqR_TDIVN0/2~0 . 5) * (VT(ct) -X0BETA(ct , i) ) ) i 

end 

end 

PEFINAL  =  PE/8; 
figure (1) 

semilogy (RPSQR_TDIVNO_DB, SINGCHAN, ,RPSQR_TDIVNO_DB ,PEFINAL(1 , : ) , . . . 

RPSQR.TD IVNO.DB, PEFINAL (2, :) ,RPSQR_TDIVN0_DB,PEFINAL(3, :) , . . . 
RPSQR_TDIVN0_DB,PEFINAL(4, :)) 
xlabeK'Z  (dB)  0  ; 
ylabel('Pb') ; 
eLxis([lO  17  lO^C-lB)  1]) 
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y. 

y. 

y. 


y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


GENERATION  OF  GRAPHS 
FOR 

betatau  =  1000 


DATE:  21  OCTOBER  1994 

COMMENTS:  THIS  GRAPHING  PROGRAM  DOES  NOT  CONTAIN  THE  RESULTS  FOR 
TERM  #5,  APPENDIX  B  DUE  TO  AN  INABILITY  TO  ARRIVE  AT  A 
SOLUTION.  THERE  WERE  MEMORY  PROBLEMS  WHEN  A  MATRIX  SOLUTION 
WAS  TRIED  TO  SPEED  CALCULATION  OF  THIS  TERM,  AND  THE 
TIME  REQUIRED  TO  CALCULATE  USING  SIMPLE  LOOPS  WAS 
PROHIBITIVE.  THUS,  ON  THE  RESULTING  GRAPH  PRODUCED 
ONLY  THE  I  =  20  TRACE  IS  ACCURATE  IN  PREDICTING 
THE  SYSTEM’S  PROBABLILITY  OF  BIT  ERROR.  THIS  IS  BECAUSE 
AT  I  =  20,  THE  EFFECTS  OF  ADJACENT  CHANNEL  INTERFERENCE  (ACI) 
ARE  NEGLI GABLE. 


rho  =  0.99; 
c  =  (l-rho)"2; 

I  =  [5  6  9  20]  ; 

M  =  [198  164  110  48]; 

bitmatrix  =  [0  0  0;  0  0  1;  0  1  0;  0  1  1;  1  0  0;  1  0  1;  1  1  0;  1  1  1] ; 


term8  =  [-464.3344  -327.7899  -148.2912  -29.7545 

-464.3344  -327.7899  -148.2912  -29.7545 

-928.6688  -655.5798  -296.5824  -59.5090]; 


term9  =  [155.3047  109.6350  49.5986  9.9519 

155.3047  109.6350  49.5986  9.9519 

310.6094  219.2701  99.1971  19.9038]; 


termlO  =  [-0.0200  -0.0142  -0.0064  -0.0013 

-0.0200  -0.0142  -0.0064  -0.0013 

-0.0401  -0.0283  -0.0128  -0.0026]; 


for  ct  =  1:4 
for  i  =  1:8 

for  boo  =  [1  0] 

XI  =  boo*c*8.5126e+03; 

X2  =  bitmatrix(i,3)*c*492.4271; 

X3  =  bitmatrix(i,3)*c*4.2917e-02; 

X4  =  bitmatrix(i,3)*c*l .5865e-05; 
y.  NOTE  THAT  THE  RESULTS  FOR  TERM  #5,  APPENDIX  B  ARE 
y.  MISSING.  SEE  THE  COMMENTS  IN  THE  TITLE  SECTION 
y.  ABOVE. 

X6  =  boo*bitmatrix(i,3)*c*9 .9411e+02; 

X7  =  boo*bitmatrix(i,3)*c*7.3500e-01; 
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if  boo  ==  0 
X8  =  0; 

elseif  bitmatrix(i , 1) ==1  &  bitmatrix(i ,2)==0 
X8  =  c*term8(l , ct)  ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrix(i ,2)==1 
X8  =  c*term8(2 , ct)  ; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==1 
X8  =  c*term8C3,ct) ; 
else 

X8  =  0; 

end 

if  bitmatrixCi, 3)  ==  0 
X9  =  0; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==0 
X9  =  c*term9Cl , ct) ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrixCi ,2)==1 
X9  =  c*term9C2 , ct)  ; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==1 
X9  =  c*term9C3,ct)  ; 
else 

X9  =  0; 

end 

if  bitmatrixCi, 3)  ==0 
XIO  =  0; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==0 
XIO  =  c*termlOCl ,ct) ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrixCi ,2)==1 
XIO  =  c*termlOC2 ,ct) ; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==1 
XIO  =  c*termlOC3,ct) ; 
else 

XIO  =  0; 

end 

INTERSUM  =  X1+X2+X3+X4+X6+X7+X8+X9+X10; 
if  boo==l 

XlBETACct.i)  =  INTERSUM; 
else 

XOBETACct.i)  =  INTERSUM; 

end 

end 

end 

XlMINCct)  =  minCXlBETACct, :)); 

XOMAXCct)  =  maxCXOBETACct, :)); 

VTCct)  =  CXOMAXCct)  +  XlMINCct) )/2; 
end 

RPSQR_TDIVNO_DB  =  10:. 01: 20; 

RPSQR.TDIVNO  =  10 . “ CRPSQR_TDIVN0_DB*0 . 10) ; 

SINGCHAN  =  0.5*erfcCRPSQR_TDIVN0/8~0.5) ; 
for  ct  =  1:4 

PECct, 1:1001)  =  zerosCl,1001) ; 
for  i  =  1:8 
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PE(ct, :)=PE(ct, :)+0.25*erfc((RPSqR_TDIVN0/2~0.5)*(XlBETA(ct,i)-VT(ct))) . . 
+0 . 25*erf c ( (RPSQR_TDIVN0/2-0 . 5) * (VT(ct ) -X0BETA(ct , i) ) ) ; 

end 

end 

PEFINAL  =  PE/8; 
figure (1) 

semilogy(RPSqR_TDIVNO_DB,SINGCHAN,'— ' , RPSqR_TDIVNO_DB, PEFINAL ( 1, . . . 

RPSqR_TDIVN0_DB,PEFINAL(2, :) ,RPSqR_TDIVNQ_DB,PEFINAL(3, :) , . . . 
RPSqR_TDIVNO_DB , PEFINAL (4 , : ) ) 
xlabeK'Z  (dB)  0  ; 
ylabel(’Pb'); 
axis([l0  17  lO'C-lS)  1]) 
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THESIS  COMPUTER  WORK 


*/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y, 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


GENERATION  OF  GRAPHS 
FOR 

betatau  =  1500 

DATE:  25  OCTOBER  1994 

COMMENTS:  THIS  GRAPHING  PROGRAM  DOES  NOT  CONTAIN  THE  RESULTS  FOR 
TERM  #5,  APPENDIX  B  DUE  TO  AN  INABILITY  TO  ARRIVE  AT  A 
SOLUTION.  THERE  WERE  MEMORY  PROBLEMS  WHEN  A  MATRIX  SOLUTION 
WAS  TRIED  TO  SPEED  CALCULATION  OF  THIS  TERM,  AND  THE 
TIME  REQUIRED  TO  CALCULATE  USING  SIMPLE  LOOPS  WAS 
PROHIBITIVE.  THUS,  ON  THE  RESULTING  GRAPH  PRODUCED 
ONLY  THE  I  =  20  TRACE  IS  ACCURATE  IN  PREDICTING 
THE  SYSTEM'S  PROBABLILITY  OF  BIT  ERROR.  THIS  IS  BECAUSE 
AT  I  =  20,  THE  EFFECTS  OF  ADJACENT  CHANNEL  INTERFERENCE  (ACI) 
ARE  NEGLI GABLE. 


rho  =  0.99; 
c  =  (l-rho)“2; 

I  =  [7  9  12  20]; 

M  =  [212  164  124  74]; 

bitmatrix  =[000;001;010;011;  100;  101;  110;  111]; 


terms  =  [-3.5253g+02  -2.1854e+02  -1.2475e+02  -4.5096e+01 

-3.5253e+02  -2.1854e+02  -1.2475e+02  -4.5096e+01 

-7.0505e+02  -4.3708e+02  -2.4950e+02  -9 .0192e+0l] ; 


term9  =  [117.9019  73.0900  41.7221  15.0822 
117.9019  73.0900  41.7221  15.0822 
235.8037  146.1800  83.4441  30.1644]; 


termlO  =  [-l.OOOOe-04  -0.6199e-04  -0.3539e-04  -0.1279g-04 

-l.OOOOe-04  -0.6199e-04  -0.3539e-04  -0.1279e-04 

-0.2000e-03  -0.1240e-03  -0.0708e-03  -0 . 0256e-03] ; 


for  ct  =  1:4 
for  i  =  1:8 

for  boo  =  [1  0] 

XI  =  boo*c*9.0083e+03; 

X2  =  bitmatrixCi, 3) *0*328.3413; 

X3  =  bitmatrix(i ,3)*c*l .8815e-04; 

X4  =  bitmatrix(i,3)*c*7.2482e-10; 

y.  NOTE  THAT  THE  RESULTS  FOR  TERM  #5,  APPENDIX  B  ARE 
y.  MISSING.  SEE  THE  COMMENTS  IN  THE  TITLE  SECTION 
y.  ABOVE. 

X6  =  boo*bitmatrix(i,3)*c*6.6331e+02; 

X7  =  boo*bitmatrix(i,3)*c*5 . 1105e-03; 
if  boo  ==  0 


X8  =  0; 

elseif  bitmatrixCi , 1)==1  & 
X8  =  c*term8(l,ct) ; 
elseif  bitmatrixCi , 1)==0  & 
X8  =  c*term8(2,ct) ; 
elseif  bitmatrixCi , 1)==1  & 
X8  =  c*term8C3,ct) ; 
else 

X8  =  0; 

end 

if  bitmatrixCi, 3)  ==  0 
X9  =  0; 

elseif  bitmatrixCi , 1)==1  & 
X9  =  c*term9Cl,ct) ; 
elseif  bitmatrixCi , 1)==0  k 
X9  =  c*term9C2,ct) ; 
elseif  bitmatrixCi , 1)==1  k 
X9  =  c*term9C3,ct) ; 
else 


bitmatrixCi ,2)==0 
bitmatrixCi ,2)==1 
bitmatrixCi ,2)==1 


bitmatrixCi ,2)==0 
bitmatrixCi ,2)==1 
bitmatrixCi ,2)==1 


X9  =  0; 

end 

if  bitmatrixCi ,3)  ==0 
XIO  =  0; 

elseif  bitmatrixCi , 1)==1  &  bitmatrixCi ,2)==0 
XIO  =  c*termlOCl ,ct) ; 

elseif  bitmatrixCi , 1)==0  &  bitmatrixCi ,2)==1 
XIO  =  c*termlOC2 ,ct) ; 

elseif  bitmatrixCi , 1)==1  k  bitmatrixCi ,2)==1 
XIO  =  c*termlOC3,ct) ; 
else 

XIO  =  0; 

end 

INTERSUM  =  X1+X2+X3+X4+X6+X7+X8+X9+X10 ; 
if  boo==l 

XlBETACct.i)  =  IMTERSUM; 
else 


XOBETACct.i)  =  INTERSUM; 

end 


end 


end 

XlMINCct)  =  minCXlBETACct,:)); 
XOMAXCct)  =  maxCXOBETACct, :)); 
VTCct)  =  CXOMAXCct)  +  XlMINCct) )/2; 
end 


RPSQR_TDIVNO_DB  =  10:. 01: 20; 

RPSQR.TDIVNO  =  10 . ‘ CRPSQR_TDIVN0_DB*0 . 10) ; 

SINGCHAN  =  0.5*erfcCRPSQR_TDIVN0/8‘0.5) ; 
for  ct  =  1:4 

PECct , 1 : 1001)  =  zerosCl , 1001) ; 
for  i  =  1:8 

PECct,  :)=PECct,  :)+0.25*erfcCCRPSQR_TDIVN0/2''0.5)*CXlBETACct,i)-VTCct)))  .  . 
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+0 . 25*erf  c  ( (RPSQR_TDIVN0/2''0 . 5)  * (VT(ct)  -X0BETA(ct ,  i)  ) )  J 

end 

end 

PEFINAL  =  PE/8; 
figure (1) 

semilogy(RPSQR_TDIVNO_DB,SINGCHAN, '  — ' , RPSQR_TDIVNO_DB, PEFINAL ( 1, :) , . . . 

RPSQR_TDIVN0_DB,PEFIMAL(2, :) ,RPSqR_TDIVN0_DB,PEFINAL(3, :) , .  .  . 
RPSQR_TDIVN0_DB,PEFINAL(4, :)) 
xlabeK'Z  (dB)'); 
ylabel(’Pb') ; 
axis([10  17  10~(-15)  1]) 
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THESIS  COMPUTER  WORK 


y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 

7. 

y. 
y. 
y. 

y. 

y. 

y.; 

rho  =  0.99; 
c  =  (l-rho)''2; 

I  =  [8  9  12  20] ; 

M  =  [248  220  164  98] ; 
bitmatrix  =  [000;  001;  010;  Oil;  100;  101;  110;  111]; 

term8  =  [-351.8931  -283.0055  -163.9044  -60.1542 
-351.8931  -283.0055  -163.9044  -60.1542 
-703.7863  -566.0110  -327.8088  -120.3084]; 

term9  =  [91.5258  73.0261  41.6758  15.0544 
91.5258  73.0261  41.6758  15.0544 
183.0515  146.0523  83.3516  30.1088]; 

termlO  =  [-6.5584e-07  -5.2745e-07  -3.0548e-07  -1.1211e-07 

-6.5584e-07  -5.2745e-07  -3.0548e-07  -1.1211e-07 

-1.3117e-06  -1.0549e-06  -6.1096e-07  -2 . 2423e-07] ; 

for  ct  =  1:4 

for  i  =  1:8 

for  boo  =  [1  0] 

XI  =  boo*c*9 . 2563e+03; 

X2  =  bitmatrix(i,3)*c*246.2563; 

X3  =  bitmatrix(i,3)*c*9 . 2720e-07; 

X4  =  bitmatrix(i,3)*c*3.2152e-14; 

y.  NOTE  THAT  THE  RESULTS  FOR  TERM  #5,  APPENDIX  B  ARE 
y.  MISSING.  SEE  THE  COMMENTS  IN  THE  TITLE  SECTION 
y,  ABOVE. 

X6  =  boo*bitmatrix(i,3)*c*4.9749e+02; 

X7  =  boo*bitmatrix(i,3)*c*3.4503e-05; 


GENERATION  OF  GRAPHS 
FOR 

betatau  =  2000 

DATE:  15  NOVEMBER  1994 

COMMENTS:  THIS  GRAPHING  PROGRAM  DOES  NOT  CONTAIN  THE  RESULTS  FOR 
TERM  #5,  APPENDIX  B  DUE  TO  AN  INABILITY  TO  ARRIVE  AT  A 
SOLUTION.  THERE  WERE  MEMORY  PROBLEMS  WHEN  A  MATRIX  SOLUTION 
WAS  TRIED  TO  SPEED  CALCULATION  OF  THIS  TERM,  AND  THE 
TIME  REQUIRED  TO  CALCULATE  USING  SIMPLE  LOOPS  WAS 
PROHIBITIVE.  THUS,  ON  THE  RESULTING  GRAPH  PRODUCED 
ONLY  THE  I  =  20  TRACE  IS  ACCURATE  IN  PREDICTING 
THE  SYSTEM'S  PROBABLILITY  OF  BIT  ERROR.  THIS  IS  BECAUSE 
AT  I  =  20,  THE  EFFECTS  OF  ADJACENT  CHANNEL  INTERFERENCE  (ACI) 
ARE  NEGLIGABLE. 
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X8  =  0; 

elseif  bitmatrix(i , 1)==1  &  bitmatrix(i ,2)==0 
X8  =  c*term8(l,ct) ; 

elseif  bitmatrix(i , 1)==0  &  bitmatrix(i,2)==l 
X8  =  c*term8(2,ct) ; 

elseif  bitmatrix(i , 1)==1  &  bitmatrix(i ,2)==1 
X8  =  c*term8(3,ct) ; 
else 

X8  =  0; 

end 

if  bitmatrix(i ,3)  ==  0 
X9  =  0; 

elseif  bitmatrix(i , 1)==1  &  bitmatrix(i ,2)==0 
X9  =  c*term9(l ,ct) ; 

elseif  bitmatrix(i , 1)==0  k  bitmatrix(i ,2)==1 
X9  =  c*terni9(2,ct) ; 

elseif  bitmatrix(i,l)==l  k  bitmatrix(i,2)==l 
X9  =  c*term9(3,ct) ; 
else 

X9  =  0; 

end 

if  bitmatrix(i,3)  ==0 
XIO  =  0; 

elseif  bitmatrix(i , 1)==1  k  bitmatrix(i ,2)==0 
XIO  =  c*ternilO(l  ,ct)  ; 

elseif  bitmatrixCi , 1)==0  k  bitmatrixCi ,2)==1 
XIO  =  c*tennl0(2,ct) ; 

elseif  bitmatrixCi , 1)==1  k  bitmatrixCi ,2)==1 
XIO  =  c*termlOC3,ct) ; 
else 

XIO  =  0; 

end 

INTERSUM  =  X1+X2+X3+X4+X6+X7+X8+X9+X10; 
if  boo==l 

XlBETACct,i)  =  INTERSUM; 

g1  S6 

XOBETA(ct,i)  =  INTERSUM; 

end 

end 

end 

XlMINCct)  =  minCXlBETACct, :)) ; 

XOMAXCct)  =  maxCXOBETACct, :)) ; 

VTCct)  =  CXOMAXCct)  +  XlMINCct))/2; 
end 

RPSQR_TDIVNO_DB  =  10:. 01: 20; 

RPSQR.TDIVNO  =  10 . “ CRPSqR_TDIVN0_DB*0 . 10) ; 

SINGCHAN  =  0.5*erfcCRPSQR_TDIVN0/8~0.5) ; 
for  ct  =  1:4 

PECct , 1 : 1001)  =  zerosCl, 1001) ; 
for  i  =  1:8 


PE(ct, :)=PE(ct, :)+0.25*erfc((RPSqR_TDIVNQ/2~0.5)*(XlBETA(ct,i)-VT(ct))) . . 
+0 . 25*erf c ( (RPSqR_TDIVN0/2~0 . 5) * (VT(ct) -X0BETA(ct , i) ) ) ; 

end 

end 

PEFINAL  =  PE/8; 
figure (1) 

semilogy(RPSqR_TDIVNO_DB,SINGCHAN, ^ ,RPSqR_TDIVNQ_DB , PEFIMAL( 1 , : ) , . . . 

RPSqR_TDIVNO_DB, PEFINAL (2, :) ,RPSqR_TDIVN0_DB,PEFINAL(3, :) , . . . 
RPSqR_TDIVNO_DB , PEFINAL (4 , : ) ) 
xlabeK'Z  (dB)'); 
ylabel('Pb') ; 
axis([l0  17  lO'C-lB)  1]) 
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